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Measurement  of  the  water  vapor  cycle  of  the  land  surface  is  crucial  in  improving 

cvapolranspiration  in  the  hydrologic  cycle,  its  magnitude,  spatial  and  seasonal 
distributions,  and  its  relation  to  environmental  variables  remains  relatively  unknown.  In 
the  state  of  Florida,  evapotranspiration  (ET)  is  second  only  to  precipitation  in  magnitude 
within  the  hydrologic  budget. 

This  research  began  with  the  examination  of  the  spectral  responses  of  surface 
objects  in  Florida  by  the  use  of  a hand-held  radiometer.  The  relationship  between  the  field 


I the  coarse  scale  satellite  data 1 


■ extracted  from  the  NOAA 


satellite  AVHRR  data.  The  concept  of  a Normalized  Difference  Vegetation  Index 
(NDVI)  and  Diurnal  Surface  Temperature  Variation  (DSTV)  triangle  was  introduced, 
based  on  the  pattern  of  ground  objects  located  in  the  NDVI-DSTV  scattergram.  The 
NDVI-DSTV  triangle  relationship  was  used  to  calculate  the  vegetation  and  moisture 
coefficient  (VMC)  for  fractional  contribution  from  various  vegetation  cover  and  soil 
moisture  backgrounds.  A geographic  information  system  (GIS)  was  utilized  to  manage 
the  satellite  and  ground-based  data.  Regional  ET  was  then  estimated  from  the  data  layers 
in  the  GIS. 

600  m for  infrared  and  thermal  bands,  respectively.  The  concept  of  NDVI-DSTV  triangle 
was  successfully  verified  with  both  ground  and  satellite  observations  and  the  relationships 
of  NDVI-DSTV  determined  the  fractional  combination  from  agriculture,  wet  and  dry  soil 
surfaces.  Results  of  the  calculated  VMC  clearly  displayed  the  fractional  contributions  of 
crop  coefficients  from  vegetation,  wet  and  dry  soils.  Results  of  ET  estimation  show  that 
this  algorithm  did  provide  an  efficient  method  for  studying  the  spatial  variation  in  ET  and 
making  useful  estimates  of  ET  in  south  Florida. 


CHAPTER  1 
INTRODUCTION 


Evapotranspiration  (ET)  is  the  process  by  which  water  is  transferred  from  the 
earth's  surface  into  the  atmosphere.  It  includes  evaporation  from  soil  and  water  surfaces 
and  transpiration  from  vegetation.  On  land  surface,  about  60%  of  the  precipitation 
received  is  expended  through  ET  (Bnitsacrt,  1986).  Along  with  rainfall  and  runoff,  it 
controls  the  availability  of  water  at  the  earth's  surface  and  directly  affects  human 
habitation  and  activities.  Quantification  of  ET  in  the  hydrologic  budget,  therefore,  is 
essential  when  planning  for  the  development  of  water  resources  of  an  area,  estimating 
water-supply  potential,  and  understanding  the  ecological  effects  of  development.  In 
Florida,  basinwidc  studies  have  demonstrated  that  ET  is  second  only  to  precipitation  in 
magnitude  within  the  hydrologic  budget  (Jones  ct  al.,  1984).  Evapotranspiration  can  vary 
considerably  among  basins  that  exhibit  different  types  of  vegetative  cover  or  different 
proportions  of  open  water  surfaces.  Demands  for  water  supply  in  Florida  continue  to 


ain(Bidlakeetal.,  19%). 


studies  on  a regional  scale.  Hydrological  modeling  is  often  done  at  a scale  which 
involves  the  analysis  of  ET  for  a region.  Doran  (1993)  reported  that  the  need  for 
measurement  of  ET  over  large  areas  has  been  a growing  concern  as  efforts  to  use 
environmental  models  at  global  scales  are  intensified.  ET  patterns  at  a regional  scale 
integrate  factors  such  as  change  in  land  use,  rainfall  distribution  and  deforestation.  The 

studies  of  global  change  (Curran  and  Foody,  1994). 

(Penman,  1948;  Tanner  and  Pelton,  I960;  Blaney  and  Criddle,  1966;  Penman  et  «1„  1967; 
Christiansen,  1968;  Brun  et  ill.,  1972).  When  applied  to  regional  evapotranspiration 
estimation,  most  of  these  techniques,  however,  have  been  of  limited  usefulness  for  three 

measure  over  an  extended  area.  Secondly,  they  provide  point  estimation.  That  is,  these 
estimates  arc  onlv  applicable  over  a homogeneous  area  surrounding  the  location  where 
the  parameters  were  measured.  The  results  may  not  reliably  apply  to  large  heterogeneous 
areas.  Thirdly,  many  traditional  ET  methods  provide  potential  ET  or  reference  ET,  i.e., 
the  ET  with  unlimited  water  supply,  which  is  not  the  general  case.  They  have  to  depend 
on  coefficients  related  to  crop  types  or  surface  resistance  to  estimate  actual  ET.  in 
irrigation  scheduling,  the  concept  of  crop  coefficients  (Kc)  has  been  widely  used  to  relate 


efficients  in  regional  scale 


has  several  difficulties.  For  example,  coefficients  for  natural  vegetation  cover,  such  as 
hardwood  forest,  cattail  marsh,  or  cypress  swamp,  as  opposed  to  field  crops,  arc  virtually 

Even  with  the  all  crop  coefficients,  it  is  still  difficult  to  determine  the  ratio  between 
reference  ET  and  actual  ET  for  an  area  with  various  vegetation  growing  at  different 
stages.  Despite  the  instrumental  and  theoretical  progress  in  ET  estimation,  as  Hope  et  al. 
(1986)  pointed  out,  the  problem  of  estimating  regional  ET,  as  opposed  to  point  ET,  still 


Since  the  beginnings  of  Earth  remote  sensing  from  satellites,  remotely  sensed 
observations  of  the  Earth's  surface  have  provided  a promising  source  of  data  for 
examining  the  characteristics  over  land  surfaces.  Hydrologists  and  agricultural 
meteorologists  have  been  searching  for  techniques  to  calculate  regional 
evapotranspiration  from  data  provided  by  satellites.  Moran  et  al.  (1989)  suggested  that 
the  only  feasible  means  of  mapping  the  spatial  distribution  of  ET  on  a regional  scale  is  to 

The  fundamental  of  satellite  remote  sensing  is  based  on  the  knowledge  of  the 
spectral  responses  from  objects  on  the  ground.  Spectral  measurements  using  spectral 
radiometer  arc  of  crucial  importance  in  remote  sensing,  both  at  the  level  of  primary 


research  and  in  operational  applications.  For  vegetation  property  research  or  to  aid  in  the 
interpretation  of  satellite  data,  relationships  must  be  developed  between  satellite  sensed 
data  and  field  measured  spectral  data.  However,  relationships  between  field  spectral 
radiometer  and  satellite  remote  sensing  are  still  a topic  to  be  explored,  mainly  because  of 
the  spectra  mixing  problem  and  the  observational  scale  difference  (i.e„  fine  resolution, 

Field  radiometer  measurements,  while  being  close  to  the  object,  are  still  a mixture  of  the 
spectral  responses  over  the  field  of  view  of  the  instrument.  In  satellite  remote  sensing, 
the  digital  number  of  each  image  picture  clement  (pixel)  is  an  aggregation  of  all 
responses  from  objects  within  the  resolution  regardless  of  the  resolution  of  the  satellite 

optimal  spatial  resolution  for  detecting  the  characteristics  on  the  ground  for  the  subject 
area  we  wish  to  measure?  In  this  study,  the  regional  ET  is  of  interest  and  the  optimal 
spatial  resolution  of  the  relevant  optical  and  thermal  property  from  each  major  land 

Surface  temperature  and  vegetation  index  are  two  of  the  most  important  surface 
parameters  that  can  be  obtained  from  the  thermal  and  optical  portion  of  satellite  remote 
sensing,  respectively.  However,  the  relationship  of  surface  temperature  and  vegetation 


sivcly  studied.  These  I 


s.  Price  (1990)  reported  that  there  is  no  adv 


simulation  model  for  ET  estimation  if  the  surface  parameters  required  arc  not  available 
for  the  region.  The  difficulties  prompt  the  need  to  develop  a simple  algorithm  for  regional 
ET  estimation. 


The  objectives  of  this  study  are  threefold:  1)  to  study  the  field  spcctroradiomctry 

find  an  optimum  spatial  resolution  of  remotely  sensed  data  that  can  be  used  in  regional 

temperature  and  use  the  satellite-derived  NDVI  and  DSTV  for  regional  scale  ET 
estimation  in  Florida. 


CHAPTER  2 

REVIEW  OF  LITERATURE 


The  study  of  regional  ET  using  satellite  remote  sensing  included  at  least  three 
major  disciplines:  1 ) evapotranspiration  theory;  2)  remote  sensing  theory  and  methods; 

concepts  of  the  ET  process  and  traditional  methods  of  ET  estimation.  ET  characteristics 
of  Florida  were  also  reviewed.  Remote  sensing  theory  from  field  radiometry  to  satellite 
image  processing  techniques  was  reviewed.  Spectral  characteristics  of  vegetation  and 
soil,  and  the  theory  of  optical  and  thermal  remote  sensing  were  then  studied.  The  primary 
parameters  that  can  be  derived  from  remote  sensing,  vegetation  index  and  surface 
temperature,  were  then  introduced  as  a link  between  the  ET  estimation  methods  and 
remote  sensing  techniques.  Early  works  of  remote  sensing  methods  for  ET  estimation 
were  also  studied.  Finally,  the  scale  considetation  in  remote  sensing  was  also  reviewed. 


from  photosynthctically  active  vegetation.  Evapotranspiration  from  land  surfaces  is  the 
result  of  several  compound  processes  including  radiation  exchanges,  vapor  transport,  and 


2.1 


2.1.1  Physics  of  the  ET  Process 


biological  growth  operating  within  a system  ir 
There  are  three  main  factors  that  influence  the 
energy  to  provide  the  latent  heat  of  vaporizatir 
surface,  and  the  ability  to  transport  the  vapor  t 


ivolving  the  atmosphere,  plants,  and  soil, 
process  of  ET.  They  are  the  supply  of 


2.1.2  I 


Since  the  1960's,  a significant  amount  of  research  effort  has  been  devoted  to  the 

into  four  broad  categories:  energy  balance  approach,  aerodynamic  approach,  combination 
approach  and  instrument  measurement  (eddy  correlation  method). 


Evaporation  and  transpiration  require  a large  amount  of  energy  at  the  earth- 
atmosphere  interface.  Solar  radiation  usually  supplies  80  to  1 00  percent  of  this  energy 
and  is  often  the  factor  limiting  ET.  Hence,  energy  budget  is  one  of  the  earliest  methods 
used  to  estimate  evaporation.  The  energy  balance  equation,  in  general  form,  can  be 

R.  - kE  - H - G - P * Ah  = ^ (2.1) 


where  R„  - specific  flux  of  net  incoming  radiation, 

AE  = latent  heat  of  evapotranspration. 


: surface  to  the  atmosphere. 


G - specific  flux  of  heal  conducted  into  the  earth, 

P = energy  flux  stored  as  photochemical  energy  in  the  process  of 

A,  = energy  adveetion  into  the  layer  as  specific  flux, 

OM.  = the  rate  of  energy  storage  per  unit  area  in  the  layer. 

The  exact  form  of  several  of  the  terms  depends  on  the  type  of  condition  and  time 
for  which  the  energy  balance  is  described.  In  many  practical  applications,  especially  over 
land  surfaces,  the  terms  P,  Ah,  and  dM/dl  are  of  little  consequence,  so  that  Equation  2.1 
assumes  a much  simpler  form.  The  energy  balance  equation,  after  simplification,  can  be 


If  the  net  incoming  radiation,  sensible  heat  flux,  and  ground  heat  flux  can  be 
measured  or  estimated,  the  evapotranspiration  can  be  solved  as 


E = (R„-H-G)/i 


(23) 


balance  equation  for  evapotranspiration. 


Evapotranspiration  rate  can  also  be  determined  by  calculating  the  ability  to 
transport  water  vapor  away  from  the  land  or  vegetation  surfaces.  The  transport  rate  is 


r the  surface  and  I 


he  surfaces  using  Ihe  equation  developed  by  Thomt 


2. 1.2.3  £ 

energy  supply  and  by  the  energy  budget  method  when  vapor  transport  is  not  a limiting 
factor.  However,  under  normal  conditions,  both  factors  are  limiting,  so  a method 
considering  both  factors  is  required  (Chow  el  al.,  1 988).  Since  the  atmospheric  transport 
mechanisms  of  sensible  heat  are  similar  to  those  of  water  vapor,  Bowen  (1926)  assumed 
that  the  sensible  heat  flux  and  the  latent  heat  flux  are  proportional,  and  the  proportionality 
constant  being  called  the  Bowen  ratio. 


P = — 


(2.5) 


When  net  radiation  R„  and  ground  heat  flux  G are  known,  the  evaporation 


(2.6) 


10 

(«„-G) 

i(l.p) 


The  value  of  P can  be  calculated  from  gradients  of  air  temperature  and  vapor 


K„  = heat  diffusivity. 


By  substituting  0.622(e/p)  for  q,  t 


(2.10) 


where  Y is  often  called  the  psychromctr 


0.622  X A'„ 


and  the  ratio  K*  and  K»  of  the  heat  and  vapor  diffusivities  is  commonly  taken  to  be 
(Priestley  and  Taylor,  1972). 

To  determine  Bowen  ratio  p,  measurements  of  temperature  and  vapor  press 
required  at  two  different  heights.  A vety  important  step  in  Penman's  analysis  (Pern 
1 948)  is  that  he  defined  a gradient  A of  the  saturated  water  vapor  pressure  curve. 


where  e,"=e*(TJ  is  the  corresponding  saturation  vapor  pressure  at  air  temperature  T„  and 
e,'=e'(TJ  is  the  vapor  pressure  at  the  wet  surface.  Since  e,=e,.  for  a saturated  surface,  the 


anion  of  Equation  2.13  in  Equation  2,6  yields 


(2.13) 


the  gradient  A in  Equation  2.12 1 
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after  rainfall  or  dew  formation  (Brulseart  1982).  Thus,  the  specific  humidity  at  the 
surface  of  the  foliage  is  somewhat  smaller  than  the  saturation  value  at  the  corresponding 
temperature  and  Penman's  equation  is  no  longer  applicable.  This  lack  of  generality  in 
Penman's  equation  was  satisfactorily  overcome  by  Monteith  (1973)  who  combined  energy 

for  vegetation  of  any  type  in  any  state  of  water  stress.  In  Monteith's  approach,  the  surface 

foliage  elements  and  a bulk  stomatal  resistance  was  introduced  to  characterize  the  transfer 
between  the  stomatal  cavities  and  the  leaf  surface.  When  the  vegetation  is  not  actually 
wet.  the  surface  vapor  pressure  c,  is  not  equal  to  e,'.  However,  they  can  be  related  by 


(«,-o  ■ — («;-o 


(2.18) 


ratio,  and  substituting  into  Eqt 


(2.20) 


change  of  saturation  vapor  pressure  with  temperature, 

specific  heat  of  air  at  constant  pressure, 

surface  resistance  of  vegetation. 

This  equation,  commonly  called  the  Penman-Monteith  equation  (Montcith,  1973), 
not  only  considers  thermodynamic  and  aerodynamic  aspects,  but  also  includes  the 

vapor  transfer.  The  aerodynamic  resistance  describes  the  efTect  of  the  physical  roughness 
of  the  evaporating  surface  on  the  transfer  of  energy  and  mass  to  atmosphere.  The  surface 
resistance,  which  is  primarily  related  to  the  bulk  stomatal  resistance  (Thom,  197S), 
describes  the  biological  control  over  the  rate  of  transpiration.  Brutsaert  (1982)  also 
indicated  that  when  this  concept  is  used  in  the  context  of  ET  over  dty  soil  it  may  also 
account  for  the  resistance  through  the  upper  layer  of  the  soil. 

2.1 .2.4  Eddy..eorrelq|iofl.melht).d 

The  eddy  correlation  is  a direct  measurement  of  mass  transfer  that  requires 
complex  n rumen  a n and  well-trained  personnel  to  obtain  accurate  results  (Tanner  et 
al„  1985).  The  mean  vertical  (lux  of  water  vapor,  based  on  the  eddy  flux  theory,  is 


(2.21) 


wind  at  height  z;  and 

mean  vapor  pressure  at  height  i. 

The  instruments  required  for  the  eddy  correlation  method  are  a vertical 

feasible  for  field  use  in  routine  measurement  with  adequate  fetch  (Rosenberg  eta].,  1983; 
Weaver  et  al„  1986).  Fetch  is  adequate  if  there  is  sufficient  upwind  homogeneity  in  the 
vegetative  cover  so  that  the  surface  layer  equilibrium  with  the  vegetative  cover  from 
ground  surface  to  at  least  the  height  of  the  instrument.  This  upwind  distance  generally  is 
considered  to  be  about  100  times  the  instrumentation  height  (Stannard,  1993). 

Potential  ET  is  defined  as  the  quantity  of  water  evaporated  from  an  extensive  free 
water  surface  as  dictated  by  existing  atmospheric  conditions  (Shuttlcworth,  1993).  A 
concept  similar  to  potential  ET  is  reference  ET.  Reference  ET  is  defined  as  ET  from  a 
well  watered  reference  crop  mainurined  at  a given  standard  height.  Therefore,  potential 


ET  can  be  defined  exclusively  Iron 
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ological  conditions,  while  reference  ET  is 
dependent  on  the  reference  crop  selected.  Typically  grass  is  used  as  the  reference  crop. 
However,  other  crops,  such  as  alfalfa,  have  also  been  used  (Jensen  ct  al„  1 990). 

Numerous  studies  have  been  devoted  to  the  estimation  of  ET  from  meteorological 
observations,  Jensen  et  al.  (1 990),  Brutsaert  (1 982)  and  Rosenberg  et  al.  (1983).  These 

The  ET  from  a particular  crop  at  a specified  growth  stage  is  often  related  to  the 
ET  of  a reference  crop  through  the  use  of  a crop  coefficient,  Kc  (Jensen,  1 968).  The  crop 

and  Pruitt  (1977)  compiled  a list  of  crop  coefficients  for  a variety  of  crops.  Wright 
(1982)  also  derived  crop  coefficients  for  various  crops  al  Kimberly  Idaho,  from  two 

crop.  Therefore,  when  selecting  and  using  published  crop  coefficients,  the  correct 
reference  crop  ET  must  be  used  to  obtain  a reliable  estimation.  Since  they  are  empirically 
developed,  crop  coefficients  are  also  location  dependent  Hargreaves  (1994)  suggested 
that  crop  coefficients  should  be  applied  under  climatic  conditions  similar  to  the  condition 
where  they  were  developed. 

2.1.5  EyaBMcmspinilion  Studies  in  Florida 

In  Florida,  the  unevenly  distributed  rainfall  and  ET  has  been  well  recognized 
(Fcmald  and  Patton.  1984;  Smajstrla  cl  al.,  1984).  Statewide  rainfall  averages  150  billion 
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gallons  per  day  and  daily  average  ET  is  estimated  to  be  about  1 10  billion  gallons  with  the 
rate  being  the  highest  in  southern  mainland  Florida  and  lowest  in  the  Keys  and  northwest 
Florida  (Femald  and  Patton,  1984).  A study  of  potential  ET  characteristics  in  Florida 
conducted  by  Smajstrla  et  al.  (1984)  concluded  that  the  annual  potential  ET  is  1406  nun 
at  West  Palm  Beach  yet  only  1 175  mm  at  Tallahassee.  Seasonal  variability  was  found  to 
be  10%  in  the  summer  and  30%  in  winter  months. 

Evapotranspiration  from  local  and  water  basin  studies  within  Florida  also  showed 
high  spatial  and  temporal  variability.  Shih  et  al.  (1983)  conducted  a basinwidc  monthly 
ET  study  using  Penman,  pan  evaporation,  Thomthwaite,  Blaney-Criddle,  and  water 
budget  methods  in  the  Everglades  Agricultural  Area,  They  concluded  that  the  Penman 
method  provided  the  closest  prediction  while  pan  evaporation  and  the  Blaney-Criddle 
method  could  be  used  for  limited  data  availability.  ET  data  requirements  in  south 
Florida  were  studied  by  Shih  (1981, 1985)  and  the  results  showed  that  solar  radiation  and 
air  temperature  are  the  most  important  factors  among  numerous  other  climatological  data 
(pan  evaporation,  rainfall,  high  and  low  wind  velocity,  minimum  and  maximum  soil 
temperature,  period  moisture  remains  on  foliage,  relative  humidity).  According  to  the 
results  (Shih,  1981),  the  ET  was  influenced  68%  from  solar  radiation,  25%  from  mean  air 
temperature,  5%  from  wind  speed,  and  2%  from  relative  humidity.  This  indicates  that  if 
sufficient  climatological  data  are  not  available  for  using  a sophisticated  model,  an 
alternative  method  should  be  chosen  based  on  the  solar  radiation  and  temperature. 
Knowles  ( 1 996)  estimated  ET  for  the  Rainbow  and  Silver  Springs  basins  using  a regional 
water-budget  approach.  He  concluded  that  ET  accounted  for  about  73.3%  (37.9  inches) 


of  the  annual  rainfall  of  5 1 .7  inches  over  a 30-year  period  and  die  wet  season  (June 
Ihrough  September)  ET  rate  was  about  twice  the  ET  rate  during  the  dry  season  (October 
through  May). 

ET  of  crops  and  vegetation  in  Florida  were  also  studied  by  several  researchers. 

and  sugarcane  leaf  area  index  (LAI)  in  the  Everglade  Agricultural  Area.  They  found  that 

difference  in  evaporation  caused  by  elevation  occurred  at  the  period  when  the  LAI  was 
less  than  one.  Boman  (1994)  studied  ETorFlatwoods  citrus  trees  with  four  water  table 

significandy  affect  ET  rates.  A study  of  ET  from  areas  of  native  vegetadon  in  west- 
central  Florida  was  conducted  by  Bidtakc  ct  al.  (1996).  Using  energy-balance  Bowen 
ratio  (EBBR)  and  eddy  correlation  methods,  he  estimated  1,010  mm  per  year  for  the  dry 
prairie  type,  990  mm  per  year  for  the  marsh  type,  1060  mm  per  year  for  the  pine  flatwood 
type,  and  970  nun  per  year  for  the  cypress  swamp  type.  Another  study  of  ET  from 
vegetadon  (natalgrass,  dog  fennel,  dwarf  horsewccd,  and  ragweed)  was  conducted  by 
Sumner  ( 1 996)  in  a deforested  area  of  Lake  Wales  Ridge,  Florida.  He  estimated  daily 
ET  ranging  from  0.2  mm  per  day  in  late  December  1993  to  about  5 mm  per  day  in  mid- 
July  1994.  He  also  suggested  that  the  simple  Priestley-Taylor  model  would  be  preferable 

Taylor  model  is  a good  estimator  of  ET,  is  computationally  less  complex,  and  requires 
fewer  parameters. 
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Field  spcctroradiometry  is  a technique  of  fundamental  importance  in  remote 
sensing,  both  at  the  level  of  primary  research  and  in  operational  applications.  For 
fundamental  vegetation  property  research  or  to  aid  in  the  interpretation  of  satellite  dat 
relationships  must  be  developed  between  remotely  sensed  spectral  data  and  agronomii 
variables  such  as  leaf  area  index,  biomass,  and  fractional  ground  cover. 


2.2.1.  l Spectrum  and  scwlral  Kscona  swtes 

Wavelengths  of  electromagnetic  energy  arriving  at  the  earth's  surface  are  either 
reflected,  absorbed  or  transmitted.  The  proportions  of  the  energy  reflected,  absorbed  or 
transmitted  will  be  dependent  on  the  type  and  condition  of  features  on  the  earth  surface. 
The  proportion  of  reflected  energy  will  also  vary  between  wavelengths,  so  that  features 
similar  in  certain  spectral  bands  can  often  be  clearly  differentiated  in  other  spectral  bands. 

Different  objects  generally  have  spectral  reflectance  curves  of  different  shapes  and  this 
forms  a basis  for  identifying  the  object  from  remotely  sensed  data.  The  shape  of  the 
spectral  reflectance  curve  of  an  object  is  often  called  its  spectral  signature.  Spectral 
reflectance  curves  of  typical  vegetation,  soil  and  water  in  the  visible  to  middle  infrared 
bands  of  the  electromagnetic  spectrum  arc  shown  in  Figure  2.1.  It  is  obvious  that  these 


2.2. 1.2 


pigmented,  water-filled  cells  and  air  spaces.  Each  of  the  three  features  (pigmentation, 
physiological  structure,  and  water  content)  have  an  effect  on  the  reflectance,  absorption, 
and  transmittance  properties  of  a green  leaf.  Leaf  physiology  determines  infrared 

primary  pigments  are  contained  in  higher  plants,  namely,  chlorophyll  a,  chlorophyll  b. 
carotene  and  xanthophyll,  all  of  which  absorb  visible  light  for  photosynthesis. 
Chlorophyll  a absorbs  at  wavelengths  of  0.43  pm  and  0.66  pm  and  chlorophyll  b at 
wavelengths  0.46  pm  and  0.63  pm.  The  carotenoid  pigments,  carotene  and  xanthophyll. 
both  absorb  blue  to  green  light  at  a number  of  wavelengths  (Wittingham,  1974;  Curran, 

1 983).  It  is  the  combined  effects  of  the  leaf  pigments  and  the  physiological  structure  that 
produce  the  characteristic  reflectance  property  for  vegetation.  The  spectral  character  of 
light  reflected  by  a normal  plant  leaf  is  controlled  by  two  groups  of  cells  (Figure  2.2). 

The  long,  narrow  cells  below  the  upper  epidermis  arc  the  palisade  cells.  They  contain 
many  chloroplasts  with  chlorophyll  pigments  that  absorb  most  blue  and  red  light  for 
photosynthesis  and  efficiently  reflect  green  light.  The  chlorophyll  absorption  features  arc 
centered  at  0.43  pm  and  0.68  pm,  while  the  area  between  the  two  features  falls  within  the 
portion  of  the  visible  spectrum  which  results  in  the  green  color  of  plants.  Leaf  reflectance 
is  reduced  also  as  a result  of  absorption  by  three  major  water  absorption  bonds  that  occur 
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dial  occur  near  wavelengths  of  0.96  pm  and  1 . 1 |im  (Cumin.  1983).  Table  2. 1 shows  the 
relationship  between  multispcctral  reflectance  and  vegetation  amount  for  five  wavebands 
as  defined  by  Tucker  and  Maxwell  (1976)  for  a grass  canopy.  These  relationships 
provide  indications  for  the  detection  of  vegetation  amount  by  remote  sensor  in  different 
wavelengths. 

The  vegetation  spectral  reflectance  curve  also  provides  a good  indication  for 
monitoring  vegetation  condition.  The  steep  rise  in  reflectance  of  vegetation  beyond  0.8 
pm  is  called  the  red  edge  of  the  chlorophyll  band  (Horlcr  ot  al.,  1983)  and  the  slope  of  the 
red  edge  has  been  related  to  chlorophyll  concentrations  in  the  leaves  (Horleret  al„  1983). 

(Woolcy,  1971;  Sinclair  etal.,  1973),  however,  the  breakdown  of  plant  pigments  causes  a 
rise  in  blue  and  red  reflectance.  This  results  in  a positive  relationship  between  all 

is  damaged  and  the  high  albedo  between  0.7  and  1 .3  pm  is  reduced,  which  provides  a 
diagnostic  remote  sensing  technique  for  assessing  the  condition  of  plants.  The  abrupt  and 

number  of  techniques  for  assessing  the  amount  of  vegetation  present  in  a remotely  sensed 
image  (Hobler  etal.,  1983). 


23 


Figure  2.2  Generalized  diagram  of  a leaFs  structure  and  its  reflectance  characteristics 
at  visible  and  near  IR  wavelengths  (adapted  from  Avery  and  Berlin,  1985). 


Strong  chlorophyll 
and  carotenoid 


500-600  tun 
Orange/red 


* adapted  from  Tucker  and  Maxwell,  1976. 

Spectral  reflectance  curves  of  vegetation  can  also  be  used  for  early  stress 
detection.  When  a healthy  plant  becomes  strained  by  some  type  of  stress  in  its 
environment  changes  occur  in  its  spectral  curves.  These  stresses  include  moisture 
deprivation,  prolonged  inundation  by  flood  water,  nutrient  deficiency,  plant  diseases,  and 
concentration  of  soil  salts  (Murtha,  1978).  Figure  2.3  illustrates  the  changes  in  the  typical 
reflectance  properties  of  a plant  leaf  as  it  progresses  from  a healthy  state  through  different 
stages  of  damage. 


25 


2.2.1.3 


Condit  (1970)  has  classified  160  soil 


types  according  to  the  shape  of  their  reflectance  curves  within  the  0.4  to  1 .0  pm  region  of 
the  spectrum  (Figure  2.4).  Type  1 curves  have  rather  low  reflectances  with  slightly 


1.0  pm.  Type  2 curves  are  characterized  by  a generally  decreasing  slope  to  about  0.6  pm 
followed  by  a slight  dip  from  0.6  to  0.7  pm  with  a continued  decreasing  slope  beyond 
0.75  pm.  This  results  in  a typical  convex  shape  from  the  visible  to  beyond  1 .0  pm.  Type 
2 soils  are  better  drained  and  lower  in  organic  matter  than  type  I soils.  Type  3 curves 
have  a slightly  decreasing  steep  slope  to  about  0.6  pm  followed  by  a slight  dip  from  0.62 
to  0.74  pm  with  slope  decreasing  to  near  zero  or  becoming  negative  from  0.76  to  0.88 
pm.  Beyond  0.88  pm  to  1.0  pm  the  slope  increases  with  wavelength.  Type  3 soils  have 
moderately  high  iron  content.  Condit  was  able  to  reproduce  these  curves  with  a high 
degree  of  accuracy  from  measurements  at  live  narrow  bandwidths  (0.02  pm)  centered  at 
0.40. 0.54, 0.64, 0.74,  and  0.92  pm.  Stoner  and  Baumgardner  ( 1 980)  established  two 
more  types  of  soil  reflectance  curves  similar  to  type  3 by  extending  the  data  out  to  1.3 
pm.  The  type  4 reflectance  behavior  from  0.88  to  1 .3  pm  was  caused  by  high  iron 
content  and  organic  material.  In  type  5,  the  negative  slope  from  0.75  to  1 .3  pm  resulted 
from  very  high  iron  and  low  organic  concentrations,  fills  was  the  only  type  that  did  not 


i form  from  0.32  to  about 
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Wavelength  (pm) 


Figure  2.4 


Typical  soil  reflectance  curves  for  five  major  types  of  curves.  Types  1-3 
were  proposed  by  Condil  (1970)  and  types  4 and  5 by  Stoner  and 
Baumgardner  (1980). 
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Although  reflectance  in  all  spectral  regions  is  negatively  correlated  with  organic, 
the  region  around  0.57  pm  is  particularly  useful  for  monitoring  organic  matter  in  bare 
soils  since  it  is  free  of  other  major  disturbances.  Stoner  and  Baumgardner  considered 
measurements  at  0.7, 0.9,  and  1 .0  pm  to  be  essential  for  thorough  classification  of 
background  soil  reflectance.  Absorptions  at  0.7  and  0.9  pm  are  produced  by  ferric  iron 
compounds,  white  that  at  1 .0  pm  is  caused  by  ferrous  iron  compounds. 

The  0.4  to  1.0  pm  region  is  not  useful  for  monitoring  soil  moisture  content 
(Reginalo  et  al.,  1977)  although  the  entire  reflectance  curve  is  generally  suppressed  with 
increased  moisture.  The  region  centered  al  2.2  pm  displays  the  highest  correlation  with 
soil  moisture  and  is  similarly  important  with  vegetation. 

Spectral  reflectance  of  water  is  distinctive  from  vegetation,  soils  and  most  natural 
targets  on  the  ground.  The  most  significant  characteristics  of  water  reflectance  curve  is 
the  energy  absorption  at  near-infrared  wavelengths  (0.7  to  1 .4  pm).  Therefore,  locating 

near  infrared  wavelengths  (Lillesand  and  Kiefer,  1987).  However,  many  factors,  such  as 
material  suspended  in  the  water,  vegetation  existed  in  the  water,  turbidity,  soil  or  material 
at  bottom  of  the  water  body,  and  sun  glint,  can  affect  the  reflectance  of  water. 

2.2.1. 4 Field  snectroradinmeler  measurement 

Natural  targets  are  usually  illuminated  by  the  whole  hemisphere  of  the  sky  and 
thus,  receive  direct  solar  flux  and  scattered  sky  light  Interactions  al  the  surface  result  in  a 


eflected,  either  directly  from  the  surface,  < 
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after  multiple  interactions  within  the  surface  if  the  material  is  translucent  to  the  incoming 
radiation.  Natural  targets  are  generally  not  perfectly  diffuse  (Lambertian)  reflectors  and 
thus,  the  intensity  of  the  reflected  flux  varies  with  the  angle  with  which  it  leaves  the 
surface.  Consequently,  the  radiation  environment  comprises  two  hemispherical 
distributions  of  electromagnetic  radiation,  one  incoming  and  one  outgoing  and,  it  is  the 

Radiation  geometry  of  the  field  environment  is  shown  in  Figure  2.5.  Note  that  the 

two  angles,  first  the  angle  from  the  vertical  (the  zenith  angle,  6)  and  secondly,  the  angle 

energy  from  the  Sun  and  energy  reflected  to  the  sensor  can  be  thought  of  as  being 
confined  to  two  slender  elongated  cones,  each  subtending  a small  angle  at  the  target 
surface,  termed  solid  angles,  and  measured  in  steradian  (sr).  If  these  solid  angles  are 
sufficiently  small,  the  reflectance  of  the  target  can  be  defined  as: 


(2.22) 


where  dL  is  the  reflected  radiance  per  unit  solid  angle  and  dE  is  the  irradiance  per  unit 
solid  angle,  and  the  subscripts  i and  r denote  incident  and  reflected  rays  respectively. 
Both  the  radiance  and  the  irradiance  vary  in  zenith  and  azimuth;  hence,  to  specify 
completely  the  reflectance  field  at  the  target  the  reflectance  must  be  measured  at  all 


ting  in  the  ‘bi- 
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surface  is  nol  possible  and  an  alternative  is  found  through  the  standardization  of  reflected 
radiance  by  the  use  of  a panel  specified  to  be  perfectly  diffuse,  completely  reflecting  and 
viewed  under  the  same  irradiation  conditions  and  in  the  some  geometry  as  the  target.  The 

2.2.2  Vegetation  Indices 

A spectra]  vegetation  index  is  obtained  by  ratioing,  differencing,  combining  or 
transfonr  ig  spectral  data  to  represent  plant  canopy  characteristics.  Its  purpose  is  to 
differentiate  vegetation  from  the  soil  background  and  to  provide  a numerical  value  that 
can  be  related  to  the  various  plant  parameters.  The  vegetation  index  has  been  found  to  be 
very  closely  related  to  the  green  leaf  area  index  (LAI)  (Holben  et  al-  1980;  Hatfield  et  al., 
1985),  green  biomass  (Tucker  et  al„  1981;  Hardisky  et  al„  1984;  Anderson  et  al„  1993), 
the  amount  of  absorbed  photosynthetically  active  radiation  (Gallo  et  al.,  1985; 

Choudhuty,  1987),  fractional  ground  cover  (Otmsby  et  al.,  1987),  and  vegetation  amount 
(Price,  1992). 

The  development  of  vegetation  indices  is  guided  by  three  general  objectives 
including  (1)  to  enhance  the  relevant  vegetation  feature,  (2)  to  standardize  the 
representation  of  vegetation  spectral  response  (useful  in  region  to  region  or  year  to  year 


I from  several  data  sets)  (Malin 


a,  1988). 
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yews.  It  is  not  possible  in  this  study  to  cover  all  of  the  vegetation 
— ever  created.  However,  Perry  and  Lautenschlager(1984)  analyzed  the 
matical  relationships  of  some  four  dozen  vegetation  indices  derived  for  Landsat 
spectral  scanner  (MSS)  data  and  concluded  that  most  of  them  are  fiinctionally 


rhe  normalized  difference  vegetation  index  (NDVI)  is  defined  as 
ce  between  the  near  infrared  and  red  reflectances  to  their  sum. 


IR-R 


(223) 

by  the  sum  of  radiances  and 


The  index  is  called  normalized  because  it  is  divided 
normalizes  somewhat  for  differences  in  solar  spectral  radiances.  Vegetation  areas 
nerate  high  values  because  of  their  relatively  high  infrared  reflectance  and  low 
reflectance.  In  highly  vegetated  scene  areas,  the  NDVI  typically  ranges  from  0. 

. in  proportion  to  the  density  and  greenness  of  the  plant  canopy.  Clouds,  water,  at 
which  have  larger  visible  reflectance  than  near  infrared  reflectances,  will  yield 


The  DVI  is  defined  as  the  difference  between  the  near  infrared  and  red 
reflectance. 

on  ~ IR  - R (2.24) 

It  is  similar  to  NDVI  except  without  normalization  (by  dividing  the  sum). 
However,  NDVI  is  preferred  to  the  simple  DVI  for  vegetation  monitoring  because  the 
NDVI  compensates  for  changing  illumination  conditions,  surface  slope,  aspect,  and  other 


*»  ‘ J (2.25) 

can  provide  a good  measure  for  the  vegetation  area;  however,  it  may  result  in  irrationally 
large  values  for  low  IR  reflectance  pixels  such  as  water  bodies.  RVI  is  related  to  NDVI 


The  TVI  has  the  same 


t NDVI  except  it  added  0.5  to  the  NDVI  to 


m -- 


(2.26) 
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(Deeringetal.,  1975). 

2.2.2.S  EgnmilwiilMYaflMignimteii  IPVQ 

In  agricultural  scenes,  reflections  from  individual  plants  or  individual  rows  of 
plants  are  closely  intermingled  with  the  bare  soil  between  plants  and  between  rows  of 
plants  so  that  reflectances  are  mixed,  even  at  the  finest  resolutions.  Mixing  of  soil  and 

Richardson  and  Wicgand  (1977)  demonstrated  that  a plot  of  near  infrared  versus  red  band 

radiance  decreases  and  the  near  infrared  increases.  A vegetation  point  would  lie  away 
from  the  soil  line  with  the  perpendicular  distance  from  the  point  to  the  soil  line.  The 
magnitude  of  the  displacement  with  respect  to  the  soil  brightness  line,  which  is  called  the 
perpendicular  vegetation  index,  has  been  shown  to  be  related  to  the  amount  of  growing 
vegetation. 

m = pR.-IRf  - (2.27) 

where  1R,  and  1R,  are  infrared  reflectance  for  soil  and  vegetation  respectively  and  R,  and 


2.22.6 


The  "tasseled  cap"  transformation  of  Kauth  and  Thomas  (1976)  is  a linear 
transformation  based  on  the  fact  that  the  correlations  between  the  visible  and  the  near 
infrared  bands  of  Landsat  MSS  data  will  project  soil  and  vegetation  information  in 
agricultural  regions  into  a single  plane  in  multispectral  data  space.  The  transformation 

components  or  features  which  arc  directly  related  to  physical  scene  characteristics.  The 
first  component  is  similar  to  the  soil  line  in  the  Perpendicular  Vegetation  Index,  called  the 
soil  brightness  index.  The  second,  the  orthogonal  distance  from  the  line  to  a vegetation 
point,  is  the  greenness.  The  third  component  (yellowness)  is  orthogonal  to  both  soil 
brightness  and  greenness,  and  the  fourth  component  (called  "nonsuch"  because  no 
features  were  evident)  is  orthogonal  to  the  first  three.  Table  22  lists  the  coefficients  and 
combinations  of  the  four  transformed  components  for  Landsat  MSS  data. 

The  transformation  coordinates  of  yellowness  and  nonsuch  have  been  shown  to 
indicate  changes  in  atmospheric  haze  conditions  and.  therefore,  may  be  useful  for  relative 
calibration  of  the  physical  slate  of  the  atmosphere.  Jackson  ct  al.  (1 983)  further  showed 
that  yellowness  was  sensitive  to  haze  condition  and  nonsuch  was  sensitive  to  water  vapor 
absorption.  However,  disadvantages  exist  for  the  Tnssclcd  Cap  transformation.  Find, 
the  transformation  is  complex  and  time  consuming  which  is  a barrier  for  real  time  image 
processing.  Secondly,  the  transformation  combinations  in  Table  2.2  are  appropriate  only 
to  Landsat  2 MSS.  For  different  sensors,  and  probably  different  field  conditions  if  they 


i the  field ' 
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Table  2.2  Equations  for  Tasseled  Cap  transformation  combinations. 
Components  Combination 


Soil  Brightness 
Greenness 
Yellow  stuff 


0.332MSS4  + 0.603MSS5  + 0.676MSS6  + 0.263MSS7 
-0.283MSS4  - 0.660MSS5  + 0.577MSS6  + 0.388MSS7 
-0.899MSS4  + 0.428MSS5  + 0.076MSS6  - 0.04IMSS7 
-0.0I6MSS4  - 0.131MSS5  - 0.452MSS6  + 0.882MSS7 


' coefficient  based  i 


i Landsat  2 MSS  data. 


Huote  ct  al.  (1983)  found  that  at  intermediate  levels  of  vegetation  cover. 

a soil-reflected  spectral  signal  that  strongly  resembles  vegetation  spectral  signatures.  This 

the  detection  of  the  actual  vegetation  condition.  After  analyzing  a series  of  vegetation 
isolines  in  near  infrared-red  wavelength  space,  Huctc  (1988)  found  that  a shifting  of  the 


and  he  proposed  a soil  adjusted  1 


(2.28) 
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vegetation  densities  L=0.S,  and  higher  density  L=0.25.  Huete  also  stated  that  a single 


2-2.3  Ifcrmal.RenWte Sensing 

The  earth  surface  can  be  thought  to  be  a physical  system  which  is  periodically 
excited  by  the  sun  and  which  responds  in  two  ways.  First,  the  earth  rc-emits,  by 
reflection,  a portion  of  the  energy  whioh  it  receives  from  the  sun,  mainly  in  the  visible 
and  near  infrared  part  of  the  electromagnetic  spectrum.  This  part  of  the  radiation,  in 
principle,  is  the  representation  of  the  surface  characteristics  of  the  earth  which  are  uscfiil 

discriminating  soil  and  rock  types,  etc.  Secondly,  it  receives  part  of  the  solar  energy  and 
emits  directly  in  the  thermal  infrared  portion  of  the  electromagnetic  spectrum.  This 

be  used  to  infer  the  energy  exchanges  in  the  earth-atmospheric  boundary. 

2.2.3. 1 Physical  basis 

All  objects  with  temperatures  above  absolute  zero  emit  electromagnetic  (EM) 
energy.  The  magnitude  and  spectral  distribution  of  this  energy  can  be  defined  by  Plank's 
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Plank's  constant,  6.626- 1 O'"  1 si  c is  the  speed  of  light,  2.99792  10'  ms1;  2 is  the 

temperature  in  degrees  Kelvin.  This  equation  describes  the  spectral  emission  of  an  object 

equation  needs  to  be  corrected  by  multiplying  by  the  emissivity  e for  different  materials. 
The  total  radiant  emission  from  a surface  at  a given  temperature  can  then  be  determined 
by  integrating  Mx  with  respect  to  wavelength  for  all  wavelengths,  t.e. 


If-/.**./* 


(2.30) 


5.6697-10'1  Wm‘KJ  This  fourth-power  relationship  between  temperature  and  radiant 
emission  is  known  as  the  Stefan-Boltzmann  equation.  As  electromagnetic  radiation 
propagates  through  the  atmosphere,  it  will  weakly  interact  with  various  gases  or  water 
vapor  it  encounters  in  two  ways:  absorption  and  scattering.  In  absorption,  a fraction  of 
the  energy  passing  through  a volume  element  of  the  atmosphere  is  absorbed  by 

n scattering,  a fraction  of  the  radiant  energy  passing 
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through  the  atmosphere  is  redirected  by  the  atmospheric  particles.  In  either  case,  the 
energy  is  lost  from  its  original  direction  and  the  combined  effect  of  scattering  and 
absorption  is  called  atmospheric  attenuation.  Therefore,  it  is  natural  to  perform 
correction  of  the  atmospheric  effect  to  improve  the  quality  of  remotely  sensed  data. 
Atmospheric  correction  algorithms  basically  require  two  steps.  First,  the  optical 
characteristics  of  the  atmosphere  and  its  effect  can  be  modeled  by  the  radiative-transfer 
algorithms  or.  measured  by  using  special  features  of  the  target  and  the  atmosphere  (e.g., 
water  bodies  in  clear  cloudless  condition).  Secondly,  the  images  are  adjusted  by 

Considering  both  radiant  sources  in  the  ground  and  atmosphere,  the  energy  going 
through  a unit  surface  from  a unit  space  angle  per  unit  time  and  unit  spectral  interval  can 
be  written  in  the  form  of  the  radiative  transfer  equation  ( Chandrasekhar,  1 960) 
cU. 

(2.31) 

where  /,  is  the  intensity  of  radiation  at  wavelength  A passing  through  an  absorbing  and 
emitting  layer,  s the  path  length,  Ji  spectra]  source,  pHIO  the  density  of  water  vapor,  and  kx 
the  absorption  coefficient  optical  depth.  This  equation  states  that  the  incoming  spectral 

the  spectral  source  J,  can  be  reduced  to  the  blackbody  emission  of  the  layer  given  by  the 


Plank's  fiinctio 


: T.  Alter  intcgraling  Equation 
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on  ^ which  merely  depends  on  temperature 
2.3 1 along  the  complete  path  between  the  surface  and  the  top  of  the  atmosphere  with  total 
optical  depth  t,\  we  obtain  the  spectral  radiance , emerging  at  the  top  of  the  atmosphere 


h ■ tAIW-O 


/8»(7lt4))e*p(-T, 


* (>-«»)  “P(-*i> 


(2.32) 


where  Tj  is  the  optical  thickness  of  the  atmosphere  along  path  s 

h ‘ * (2.33) 

The  first  term  on  the  right-hand  side  of  Equation  2.32  represents  the  contribution 
to  the  surface  temperature  T,  via  Plank's  function  B/TJ.  Surface  emissivity  £, 
compensates  the  blackbody  emission  for  the  natural  surfaces.  In  addition,  spectral 
transmission  of  the  atmosphere  e*  further  attenuates  the  surface  emission.  The  second 
term  represents  the  contribution  to  the  received  radiance  resulting  from  the  upwelling 
radiance  emitted  by  the  atmosphere  towards  the  sensor,  and  the  third  term  is  the 
downward  atmospheric  radiance  emitted  by  the  atmosphere  that  is  reflected  (reflectance 
1-eJ  by  the  land  surface  and  attenuated  along  the  upward  path  to  the  sensor.  This 


41 

satellite-borne  sensors  are  designed  to  work  in  one  of  the  atmospheric  windows  for  which 
molecular  absorption  is  small  and  the  atmosphere  is  highly  transmissible. 

to  provide  estimates  of  sea  surface  temperature  (La  Villctte  and  Chabot,  1968;  Curtis  and 
Rao,  1969).  With  the  passage  of  time,  and  the  production  of  improved  satellite  sensors,  a 
large  amount  of  study  has  been  done  to  determine  sea  surface  temperature  from  satellite 
radiances.  For  instance,  to  estimate  sea  surface  temperature,  the  split-window  method 
utilizes  the  spectral  differential  absorption  of  the  atmosphere  in  two  different  satellite 
channels  (McClain  et  al„  1985;  Singh  ct  al„  1985;  Ho  et  al-  1986),  or  uses  the 
operational  atmospheric  sounding  system  (Scott  and  Chedin,  1981),  or  the  angular 
differential  absorption  in  one  channel  (Chedin  et  al„  1982).  The  estimation  of  sea  surface 

0.5"K  (e.g.  Barton,  1992;  Ohring.  1989).  However,  the  results  can  not  be  directly 
applied  to  land  surface  mainly  because  of  the  following  obstacles  (Price,  1 984;  Becker 
and  Li,  1990);  1)  emissivity  of  the  land  surface  is  difficult  to  define  and  is  generally 
different  from  one.  Furthermore,  land  surfaces  arc  usually  inhomogeneous,  that  is,  a large 
range  of  land  uses  may  exist  within  the  resolution  of  one  satellite  pixel  which  makes  the 
determination  of  the  effective  emissivity  difficult  (Cascllcs  ct  al„  1988);  2)  Because  the 
ground  surface  temperature  within  one  satellite  pixel  is  diverse,  it  is  not  feasible  to  seek 
accurate  ground  truth  verification  for  satellite  derived  surface  temperature.  Bare  soil,  for 
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instance,  may  vary  by  several  degrees  over  distances  of  tens  of  meters  (Vauciin  et  a). 
1982)  and  temperature  variations  between  cropped  and  fallow  areas  can  be  greatly 
diverse.  Despite  the  difficulties  mentioned,  scientists  have  been  unceasingly  searching 
for  algorithms  to  extract  surface  temperature  from  satellite  data. 

Surface  temperature  estimated  from  single  satellite  channel  is  based  on  the 
solution  of  the  radiative  transfer  equation  (Equation  2.31)  and  the  inverse  of  Plank's 
function  (Equation  2.29).  For  cloud-free  locations,  the  thermal  spectral  atmospheric 
windows  are  relatively  transparent  to  radiation  except  for  the  effects  of  atmospheric  water 
vapor  (Price,  1983).  The  integral  on  Equation  2.3 1 suggests  that  atmospheric  emission  is 
independent  of  surface  radiation.  Denoting  the  atmospheric  emission  and  surface 
reflected  atmospheric  emission  terms  by  b.  and  using  symbol  a to  represent  atmospheric 

R^-aR^  + b (2.34) 

The  equation  represents  an  aggregated  outcome  of  the  atmospheric  effects  to  the 
solution  of  radiation-transfer  equation  and  expresses  that  the  satellite  observed  radiation 
is  linearly  related  to  ground  emission. 

For  typical  surface  temperature  ranges  (e.g.,  270  to  330  degree  Kelvin),  Slater 
( 1 980)  pointed  out  that  the  Planck  equation  and  its  inverse  are  nearly  linear  at  thermal 
infrared  wavelengths.  Consequently,  the  relationship  between  the  surface  radiant 
temperature  T,  and  satellite  observed  brightness  temperature  T, , can  be  expressed  as 


iT,  + b' 


(2.35) 


The  surface  radiant  temperature,  which  is  the  atmosphere  attenuation  corrected 
black  body  temperature  at  ground  surface,  is  still  subject  to  the  influence  of  emissivity. 
For  most  natural  surfaces  in  the  10-12  pm  thermal  spectral  region,  emissivity  values 
generally  have  range  greater  than  0.95  (Salisbury  and  D’ Aria,  1 992).  Also,  at  these 
wavelengths,  Slater  (1980)  found  that  the  radiance,  R,  is  proportional  to  the  temperature 
raised  to  the  4.5th  power,  R « T*s.  Thus,  from  the  radiance  temperature  T,,  one  may 
obtain  the  surface  kinetic  (physical)  temperature  T,  using  the  relation 

no (T/>  = o(r/s  • (2.36) 

T.  - (2.37) 

Price  (1983)  further  linearized  and  combined  it  with  Equation  2.35  to 

r.  ■ ( <2 .38) 

Singh  (1985)  argued  that  the  parameters  in  the  Stefan-Boltzmann-typc  relation  are 
temperature  dependent  and  the  4.5th  power  equation  will  lead  to  signiftcant  errors  in  low 
emissivities.  He  suggested  that  the  surface  kinetic  temperature,  T„  may  derive  directly 
from  the  integration  of  Plank's  equation 


(2.39) 


The  difference  between  equations  2.38  ond  2.40,  however,  is  less  than  one  degree 
K for  most  natural  surfaces  with  cmissivity  greater  than  0.95.  For  global  or  regional 
surface  temperature  computation  the  simplicity  of  the  Stefan-Boitzmann  type  equation 
(Equation  2.38)  is  preferred.  However,  for  accurate  point  evaluation  particularly  for  the 
city  or  bore  surfaces,  the  Singh's  equation  (Equation  2.40)  or  the  split  window  technique 
is  recommended. 

2.2.3.4  Split  window  technique 

Surface  temperature  estimated  from  two  adjacent  satellite  thermal  channels, 
generally  called  split  window  technique,  was  initially  proposed  by  Anding  and  Kauth 
(1 970)  for  sea  surface  temperature  (SST)  estimation.  Split  window  technique  is  based  on 
that  radiances  measured  at  two  different  wavelengths  with  a resulting  difference  in 
absorption  and  can  be  used  to  determine  the  atmospheric  attenuation  of  the  surface 
radiance  so  that  surface  temperature  can  be  derived  (Prabhakara  ct  al.,  1974;  McMillin, 
1975).  The  general  form  of  the  split  window  surface  temperature  T,  is  obtained  by  a 
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infrared  channels. 

T,  = a,  + a,T,  + 0,7)  (2.41) 

Various  methods  (Singh  1985;  McClain  ct  al„  1983;HoetaL,  1986;  Barton  ctai., 
1 989)  have  been  proposed  to  obtain  the  coefficients  a0,  a,,  a;.  Several  split  window 
formula  which  give  the  surface  radiant  temperature  T,  as  a linear  combination  of 
temperatures  given  by  the  AVHRR's  bands  4 and  5,  T,  and  T„  are  tabulated  in  Table  2.3. 


Table  2.3  Split  window  formula  for  surface  radiant  temperature  estimation. 


T,  = T,+Z6(T,-TJ-2.2 
T,  - T4+3. 1 7(T4-T5)-0.4S 
T,  = T,+0.33(T,-T,) 

T,  = T,+2.68(T,-T,)-0.45 

T,  = T,+3.58(T,-T,)+0.87(T, -33.56) 


Dcschamps  and  Phulpin  (1980) 
McClain  etal.  (1983) 

Price  (1984) 

Li  and  McDonnell  (1988) 

Vidal  (1991) 


The  split  window  technique  is  now  used  operationally  over  oceans  with  a claimed 
accuracy  of  0.2°C  (Kerr  et  al.,  1 992).  This  accuracy  is  possible  principally  because  there 


is  vety  little  difference  between  sea  surface  temperature  and  air  temperature  near  the 


t is  the  effect  of 


ntioned  that  the  split  window  techniqu 
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surface  cmissivity,  which  is  generally  lower  thar 
all  cases.  Lagouarde  and  Kerr  (1985)  have  men 
used  directly  over  land  surface  can  result  in  errors  that  reach  6°K.  Becker  (1987)  also 
pointed  out  that  if  cmissivitics  € , and  e,  in  the  two  channels  being  considered  are 
assumed  to  be  unity,  the  error  AT  generated  on  land  surface  temperature,  by  using  the 
split  window  method,  is  significant  and  is  of  the  order  of; 

A T = 50^-^ -300^'^*  (2.42) 

where 

land  surface  temperature  calculation.  That  is,  the  classical  split  window  method  will 
provide  good  results  over  water,  slightly  less  over  fully  vegetated  areas,  and  generally 
poor  results  over  dry  bare  soil. 

Theoretical  investigation  of  the  split  window  technique  over  land  surfaces  has 
been  extensively  examined  in  the  last  few  years.  Becker  and  Li  (1990),  after  linearizing 

spectral  cmissivitics  but  not  on  atmospheric  conditions.  The  equation  and  coefficients 
they  derived  from  the  actual  emissivity  values  are 


(2.43) 


A,  = 1.274, 


P = 1 +0. 1 561 6(1  -€)/e-0.482A€/€I, 

M = 6.26+3.98(l-€)/e+38.33A6/eI, 

where  £-<e,+e,y2  and  A6=(e,-e  J/2,  (e,  and  e,  corresponding  lo  surface  eraissivities  of 
channel  4 and  channel  5 ofNOAA/AVHRR). 

Similarly,  Sobrino  et  al.  (1991),  by  linearization  of  the  Planck's  equation  and 

atmospheric  water  vapor,  viewing  angle,  and  channel  surface  emissivities.  They 
concluded  that  in  Older  to  get  an  error  below  0.4  K on  land  surface  temperature,  it  is 

These  derived  split  window  methods,  although  theoretically  possible,  are  difficult 

Data  availability,  surface  characteristics  of  applied  area,  and  the  accuracy  required  in  the 
research  are  the  key  factors  determining  the  method  in  surface  temperature  estimation 
from  satellite  data.  Single  channel  method  provides  the  most  direct  approach  while  split 
window  technique  effectively  circumvents  the  atmospheric  influence. 

Since  the  beginning  of  Earth  remote  sensing  from  satellites,  remotely  sensed 
observations  of  the  Earth's  surface  provide  a promising  source  of  data  for  synoptically 
examining  characteristics  over  land  surfaces.  Hydrologists  and  agricultural 
meteorologists  have  been  searching  for  techniques  to  calculate  ET  from  data 


i provided  by 


satellites.  The  studies  involving  remote  sensing  techniques  nnd  ET  estimation  can  be 
divided  into  four  major  categories:  I)  empirical  methods;  2)  simplified  energy  balance 
methods;  3)  numerical  boundary  layer  methods;  4)  surface  parameterization  methods  and 
5)  spatial  context  methods. 

2.3.1  EmBirisal  Meihods 

Empirical  relationship  between  the  quantities  that  can  be  measured  from  a satellite 
and  evapotranspiration  have  been  recognized  early.  Idso  et.  al.  (1975)  found  a linear 
relationship  between  evaporation  and  net  thermal  radiation.  Mencnti  (1979)  obtained 
evapotranspiration  as  a bilinear  function  of  the  remotely  estimated  temperature  and 
surface  albedo.  Seguin  and  Itier  (1983)  have  shown  that,  at  a given  location,  them  exists 
a good  correlation  between  the  midday  surface  temperature  and  daily  evapotranspiration. 
Kerr  et  al.  (1989)  found  a close  relationship  between  the  Normalized  Difference 
Vegetation  Index  (NDVI),  derived  from  NOAA  High  Resolution  Picture  Transmission 
(HRPT)  data,  and  the  actual  evapotranspiration  with  a time  lapse  of  20  days. 

simplified  energy  balance  equation  for  ET  estimation.  Hcliman  et  al.  (1 976)  used 
airborne  sensors  derived  surface  temperature  in  combination  with  ground-measured  solar 

daily  ET  based  on  an  energy  balance  equation.  Jackson  et  al.  (1977)  simplified  the 
energy  balance  equation  in  which  the  daily  evapotranspiration  value  is  given  as  a function 


temperature  both  measured  near  midday  (T,  - TJ.  Their  equation  is  expressed  as: 


ET  = Rj  -B(T,  - TJ  + A (2.44) 

the  place  of  application.  This  equation  has  become  the  subject  of  many  subsequent 
studies  (Carlson  and  Buffiim,  1989;  Caselles,  et  al..  1992;  Lagouardc,  1991;  Nieumenhuis 
el  al.,  1985;  Rambal  et  al.,  1985;  Scguin  and  Iticr,  1983;  Vidal  and  Perrier,  1989).  These 
studies  either  theoretically  justified  the  formula  or  proposed  various  means  to  determine 
the  coefficients  A and  B.  The  advantage  of  this  equation  (Equation  2.44)  is  its  simplicity, 

temperature).  Gurney  and  Hall  (1983)  used  the  Heat  Capacity  Mapping  Mission  satellite 
to  estimate  ET  over  the  Alaskan  sub-Arctic.  Thunnisscn  and  Nieuwenhuis  (1990)  used 

wind  velocity,  crop  type  and  height,  to  estimate  potential  ET  of  selected  crops.  They  also 
developed  a method  to  convert  instantaneous  crop  temperature  into  24-hour  ET  values. 
Vidal  and  Perrier  ( 1 990),  developed  a method  to  map  the  surface  temperature  from 
thermal  IR  data  using  data  from  the  NOAA  TIROS  satellite  and  to  map  the  ET  of 
sugarcane  from  surface  temperature.  Some  other  studies  similarly  utilized  oncc-daily 

evapotranspiration  based  on  surface  energy  balance  (Gurney  and  Camillo,  1984;  Hatfield 
etal.,  1983;  1984;  Rcginato  et  al 


1-,  1985;  Serafini,  1987). 


2.3.3 


Numerical  algorithm  for  deriving  realistic  values  of  evaporation  and  fluxes  at  the 
scale  of  1 00  km  by  100  km,  usually  for  climatological  purposes,  are  required  to  provide 
input  data  for  mesoscale  general  circulation  models  (GCMs).  Sophisticated  land  surface 


and  vegetation  models 

and  numerical  methods  also  have  been  employed  in  solving 

surface  energy  budget  ( 

rquation.  Wetzel  and  Chang  (1988)  utilized  one  dimensional 

boundary  layer/vegetat 

ion/soil  model  and  a explicit  numerical  method  for  the 

computation  of  grid-ce: 

ll-average  evapotranspiration  over  inhomogeneous  areas.  Taconet 

etal.  (1986)  applied  the 

■ same  boundary  layer/vegetation/soil  model  with  satellite  derived 

surface  temperatures  ar 

id  introduced  a representative  canopy  resistance  to  infer  daily 

evapotranspiration  disti 

(button.  Important  experiments  combining  field  measurements 

with  remotely  sensed  dj 

ala,  such  as  HAPEX-MOBILHY  in  France  (Andre  et  al„  1986) 

and  the  First  International  Satellite  Land  Surface  Climatology  Project  (ISLSCP)  Field 
Experiment  (FIFE)  (Sellers  el  al„  1 988)  in  USA,  has  been  designed  to  map  ET  and 


surface  parameters,  Pai 

a (1990)  used  a medium-range  forecast  model  to  simulate  the 

surface  fluxes  movemci 

at  in  the  atmospheric  boundary  layer.  Numerical  weather 

prediction  models  have 

been  used  for  estimating  surface  sensible  and  latent  heat  flaxes 

(Shukla  and  Mintz,  1982;  YehetaL,  1984). 

Several  surface  parameters  related  to  the  calculation  of  ET  can  be  determined  by 


R.=Sr-  aS, 


(2.45) 


* 6„or‘  - e,oT* 

and  T,  the  air  and  surface  temperature  (K),  respectively,  and  o is  the  Stefan-Boltzmann 
constant  (5.6697*  1 0"11  W'm'^K*4),  Among  these  parameters,  surface  albedo,  emissivity, 
and  surface  temperature  can  be  estimated  from  remote  sensing  techniques.  In  addition  to 
the  net  radiation,  soil  heat  flux  is  another  parameter  that  is  related  to  ET  and  can  be 

2.3.4. 1 Surface  albedo 

Pinker  ( 1 985)  and  Pinker  and  Ewing  ( 1 986),  considering  sun-target-sensor 
geometry  and  the  integration  of  spectral  reflectance  with  respect  to  wavelength,  satellite 
zenith  angle  and  solar  zenith  angle,  derived  an  equation  of  surface  albedo  ftom  spectral 
albedos.  Wydick  et  al.  (1987)  reported  a linear  relationship  between  the  broad  and 
narrow  band  albedo  which  can  be  approximated  by 

a = -0.70 + O.36o,  + 0.7302,  (2.46) 

where  subscripts  I and  2 denote  visible  (channel  I)  and  near-infrared  (channel  2) 
observed  albedos,  respectively. 


3.4.2  Surface  cmissivity 


Surface  cmissivity,  e,,  is  a facior  lhai  describes  how  efficiently  an  object  radiates 
energy  compared  to  a blackbody  at  the  same  temperature.  The  cmissivity  of  terrestrial 
surfaces  may  vary  significantly  as  a result  of  differences  in  soil  texture,  mineral 

Griend,  1991).  The  cmissivity  of  natural  surfaces  has  been  traditionally  determined  by 
means  of  a emissivity  box  (Bueltner  and  Kent,  1 965)  and  by  active  methods  (Becker  et  nl. 
1981).  However,  it  is  practically  impossible  to  acquire  the  overall  effective  emissivity  of 
a satellite  pixel  from  individual  measurements.  Therefore,  various  methods  (Artis  and 
Carnahan,  1982;  Dong  and  Li,  1993;  Van  dc  Griend  and  Owe,  1993)  have  been  proposed 

suggested  by  Van  de  Griend  and  Owe  (1993)  gives  the  surface  cmissivity  as  a function  of 
NDV! 

e ■ 1 .0094  + 0.047  In  (NDVI)  (2.47) 

A theoretical  based  equation  derived  from  the  Planck  radiation  formula  by  Artis 
and  Carnahan  (1982)  gives  surface  emissivity  as  a function  of  two  adjacent  thermal-band 


a(re,-r„)  I 


re  a is  constant  (a  = 1.438  * 10'2  m K),  ar 


T„,  and  T„,  the  n 
rngth  of  A,  and  i 


(2.48) 

emperature 
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2.3.4.3  Soil  heal  flux 


I flux,  G,  is  given  by  Ihc : 


ar, 

dz 


(2.49) 


where  T,  is  the  temperature  of  the  soil  and  2,  is  the  thennal  conductivity  of  the  soil.  G is 

measurements  of  the  soil  heat  flux  cannot  be  extrapolated  to  large  areas.  In  the  energy 
budget  equation,  the  soil  heat  flux  utilizes  a very  small  quantity  of  available  energy 
during  daytime  and  this  is  almost  balanced  with  nighttime  radiation  loss.  An  alternative 
approach  to  calculate  soil  heat  flux  from  remote  sensed  data  is  to  make  it  proportional  to 
another  term  in  the  energy  balance  equation.  Indeed,  many  researchers  (Choudhury  ct  al, 
1986;  Reginato  et  al.,  1985;  Idso  ct  al..  1975),  have  indicated  that  the  magnitude  of  G is 
highly  related  to  the  net  solar  radiation,  R„.  Idso  et  al.  (1975)  found  that  the  ratio  of  soil 
heat  flux  to  net  solar  radiation  is  about  0.5  for  dry  soils  and  0.3  for  wet  condition.  For 
vegetation  surfaces  under  full  cover,  Monleith  (1973)  suggested  values  for  the  ratio 
would  most  likely  be  between  0.05  and  0.1.  Reginato  et  al.  (1985)  reported  that  the  soil 
heat  flux  was  approximately  1 0%  of  net  radiation  for  bare  soil  and  sparse  canopies  and 
less  than  10%  under  a well  developed  plant  canopy.  The  empirical  relationship  between 
G/R,  and  crop  height,  h,  obtained  by  Reginato  et  al.  (1985)  was 


G = (0.1  -0. 042h)  Rb 


(2.50) 


Choudhury  et  al.  (1986)  expressed  the 
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index  (LAI)  yielding  a correlation  of  0.9.  For  regional  energy  balance  studies,  typically 
little  information  about  crop  height  and  leaf  area  is  available.  Thus,  the  forementioned 
two  methods  cannot  directly  apply  to  the  calculation  of  G by  remote  sensing  methods. 

replacement  for  plant  phytomass,  leaf  area  index  and  percent  cover  (Holben  et  al„  1980; 
Ormsby  et  al.,  1987;  Price,  1992).  In  fact,  Kustus  et  al.  (1990)  obtained  a linear 
relationship  between  G/Rn  and  the  normalized  difference  vegetation  index  (NDVI)  from 
observations  of  different  vegetation  types  and  soil  surface  conditions.  The  ground  heat 

G - (0.325  - 0.208  NDVI)  R,  (2.51) 


2.3.5  Temperature  and  Vegetation  Index  Method 

Many  studies  on  radiometric  surface  temperature  have  focused  on  the  widely 
observed  negative  correlation  between  surface  temperature  and  remotely  sensed 


irements  of  actively  transpiring  vegetation  such  as  NDVI  (Hope  et  al.,  1992; 
ni  and  Running.  1989;  and  Moran  ctal.,  1994).  Shuttleworth  and  Wallace  (1985) 
:d  the  Penman-Monteith  equation  (Montcith,  1981)  to  account  for  energy 


Gillcs  et  al.  (1 995)  used  a relation  between  NDVI  and  surface  temperature  derived 
from  multispcclral  aircraft  measurements  to  define  surface  fluxes.  Over  a large  area,  a 

distribution  of  soil  moisture  and  vegetative  cover.  Schmuggcc  and  Becker  (1991)  also 


observed  ihe  NDV1  and  surface  I 


I al.  (1995) 


showed  dial  points  of  both  low  NDVI  and  surface  temperature  correspond  to  areas  of  high 

Price  (1990)  developed  a method  for  relating  contextual  information  (the  slope  of 
the  vegetation  index  - surface  temperature  line  and  the  slope  of  the  wet  soil  - dry  soil  line) 
in  AVHRR  data  to  large  area  evapotranspiration.  Carlson  et  al.  (1990)  found  that  spatial 

variation  of  soil  water  content  modulated  by  fractional  vegetation  cover.  Based  on 
theoretical  and  experimental  evidences,  Moran  ct  al..  (1994)  proposed  a concept  termed 
the  vegetation  index/temperature  trapezoid,  which  combines  vegetation  indices  with 
composite  surface  temperature  measurements  to  allow  application  of  Crop  Water  Stress 
Index  (CWSl)  theory  to  partially  vegetated  Helds  without  knowledge  of  foliage 


2.4  Scale  Consideration  in  Remote  Sensing 


Successful  use  of  satellite  data  in  remote  sensing  studies  is  often  dependent  on  the 
correct  selection  of  spatial,  temporal,  and  spectral  scales  of  the  sensor.  At  a high  spatial 
resolution  of  several  ten's  of  meters,  a satellite  can  provide  a spatially  comprehensive 
view  of  land  surfaces.  However,  the  time  frequency  of  the  data  may  not  be  adequate  to 
monitor  hydrology  and  phenology  variations.  This  is  the  case  with  the  LANDSAT 
Thematic  Mapper  (TM)  and  the  SPOT  High  Resolution  Visible  (HRV)  radiometer,  which 
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have  repetitive  cycles  of  16  and  26  days,  respectively,  The  low  temporal  coverage  limits 
the  practical  utility  of  these  satellites  to  many  types  of  applications  such  as,  monitoring 
vegetation  dynamics,  measuring  drought  and  flood  condition,  and  detecting  regional 
environmental  change.  On  the  other  hand,  satellites  with  very  high  temporal  frequency 
but  low  spatial  resolution  sensor,  such  as  GOES  Visible  and  Infrared  Spin  Scan 
Radiometer  (VISSR)  which  have  a 30  minutes  observatory  frequency  and  an  average 
resolution  of  6 km  by  8 km,  may  not  be  adequate  for  heterogeneous  land  surface 
characteristics  and  vegetation  distributions.  The  spectral  scale  identifies  in  which  range 
the  electro-magnetic  radiation  is  able  to  be  detected  by  the  sensor.  This  is  also  an 
important  consideration  in  the  selection  of  satellite  data.  The  spectral  scale  determines 

vegetation  studies  usually  need  several  visible  and  infrared  bands  for  land  use 
classification  and  vegetation  indices  calculation.  Land  surface  temperature  and 
meteorological  studies  require  thermal  band.  The  spatial,  temporal  and  spectral  scales 
need  to  be  determined  in  remote  sensing  applications.  Table  2.4  lists  the  satellite  image 

2.4.2  Optimum  Spatial  Resolution  for  Remote  Sensed  Data 

for  large  scale  land  surface  studies.  Systems,  such  as  Landsat  Thematic  Mapper,  SPOT, 
NOAA  AVHRR,  and  new  satellites  in  planning,  furnish  users  of  digital  satellite  imagery 


age  in  Landsal,  SPOT,  GOES t 


Table  2.4 


Spatial,  temporal  and  spectral  co 


16  days 
26  days 


VISSR  6 kn 
AVHRR  1.11 


VIS,  TIR 
VIS,N1R,T1R 


more  difficult.  Several  studies  (Laity  and  Hoffer,  1981;  Cushnic,  1987;  Townshend  and 
Justice,  1988)  pointed  out  that  imagery  of  fine  spatial  resolution,  which  vary  as  a function 

necessarily  improve  per-pixel  accuracy,.  Several  factors  influence  the  choice  of  an 
appropriate  scale  for  a particular  application.  These  include  the  information  desired  about 

spatial  structure  of  the  scene  itself  (Woodcook  and  Strahler,  1 987). 


To  select  an  appropriate  data  scale,  the  spatial  structure  of  images  is  a determining 
factor  The  spatial  structure  of  images  is  expected  to  be  primarily  related  to  the 
relationship  between  the  size  of  the  objects  in  the  scene  and  spatial  resolution. 

Moellering  and  Toblcr  (1972)  used  a geographic  variance  method  for  hierarchical 
geographical  regions  to  quantitatively  analyze  the  scale  effects  of  geographic  phenomena. 
They  argued  that  geographical  processes  operate  at  different  scales  and  that  one  can 
determine  "where  the  action  is”  by  calculating  the  statistical  variances  of  the  data  at 
different  levels  of  geographic  units.  Woodcock  and  Strahlcr  (1987)  developed  a measure 
of  local  image  variance  to  aid  in  selecting  an  appropriate  image  scale  for  a variety  of  land 
cover  types.  Their  local  variance  is  measured  as  the  mean  value  of  the  standard  deviation 

computed  and  the  mean  of  these  values  over  the  entire  image  is  taken  as  an  indication  of 

study  the  heterogeneity  under  a variety  of  burning  treatments  on  a natural  prairie  area. 
They  concluded  that  their  min-max  algorithm  values  reflect  the  degree  of  textural  contrast 
throughout  the  image;  the  higher  the  resulting  textural  number,  the  greater  the  degree  of 

methods  and  the  texture  analysis  method  arc  based  on  the  assumption  that  the  variability 
of  the  geographic  data  changes  with  scale  and  resolution,  and  that  the  scale  at  which  the 


: most  of  the  process  open 
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2.4.4  Sawlliigi.jlaiignisi&aaenai  et  smdy 

The  data  required  for  evapolranspiration  study  must  fulfill  three  major  conditions: 
First,  they  must  provide  spatial  quality  that  is  sufficient  for  good  spatial  resolution  to 
identify  the  land  use/cover  on  the  ground.  Second,  the  information  mast  also  meet  the 
temporal  quality  that  satisfy  the  need  for  frequent  observation.  Thiid,  they  should  have 
spectral  bands  in  the  visible  and  thermal  range  in  order  to  provide  vegetation  and 
temperature  properties  simultaneously.  Regional  scale  evapotranspiration  studies,  due  to 
their  highly  temporal  variation  and  low  spatial  requirement,  are  often  carried  out  using 

AVHRR  of  the  NOAA  satellite. 


CHAPTER  3 

METHODS  AND  MATERIALS 
3.1  Field 


A spectroradiometer  was  used  lo  collect  spectral  responses  from  ground  objects  in 
the  field.  A collection  of  the  response  curves  of  Florida  vegetation  and  soil  were 
measured  during  1996  to  1997.  Sawgrass  and  cattail  spectral  responses  tfom  4 different 
dates  were  used  lo  1)  determine  the  most  efficient  spectral  bands  for  discriminating 
between  the  two  similar  wetland  vegetation;  and  2)  determine  the  most  desired  season  to 
distinguish  between  these  two  vegetation.  The  relationship  between  data  fiorn  field 
spectroradiometer  and  from  NOAA  satellite  was  also  established. 


In  the  field,  radiometric  reflectance  were  measured  by  a Geophysical  and 
Environmental  Research  Corporation  (GER)  Spectroradiometer  model  GER1500  (GER, 

1 996).  It  is  a portable  spectroradiometer  covering  the  ultraviolet  (UV),  visible,  and  near- 
infrared  (NIR)  wavelengths  from  350  nm  to  1050  nm.  It  uses  a diffraction  grating  with  a 

5 1 2 spectral  bands  (GER,  1 996).  With  a spectral  resolution  of  1 .7  pm,  the  GER  1 500 
spectroradiometer  collects  512  discrete  samples  from  0.35  pm  to  1.05  pm  in  0.5  seconds. 
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Up  lo  483  spectral  readings  < 


i loading  and  analysis  using 


a personal  computer  with  a standard  RS232  serial  port  and  GER  licensed  operating 
software.  A Spectralon  diffuse  while  standard  panel  (SRT-99-50),  calibrated  by 
Labspherc  (Labsphcrc,  Inc.,  North  Sutton.  NH),  was  used  for  reflectance  calibration.  The 
reflectance  of  a target  can  be  calculated  as 

<■-'> 

where  = spectral  reflectance  of  a target  (%); 

«a<c  = target  radiance  (W  cm‘!  nm'  sf 1 * 10-'°); 

Rad„f  - reference  radiance  (W  cm"!  nm"1  sr"1  * IQ"10). 

The  reference  radiance  is  a measurement  of  the  Spectralon  panel  under  the  same 
light  and  geometric  condition  as  the  target  measurement.  Spectralon  diffuse  reflectance 
material  is  a commercially  available  polytetrafluorocthylcne  (PTFE)  compound  which 
maintains  strong  spatial  and  spectral  uniformity  and  stability  over  time  (Brucgge  et  al., 
1995).  It  provides  a calibrated  reflectance  values  of  95%  to  99%  and  is  spectrally  flat 
over  the  UV-VIS-NIR  spectrum. 

Spectral  reflectance  from  vegetation  and  soil  was  collected  with  the  GERI 500 
spectroradiometer.  Each  sample  was  measured  three  to  live  limes  from  different  viewing 

each  wavelengh  was  calculated  over  the  measured  target  reflectances  and  90% 
confidence  intervals  of  the  reflectance  curve  were  calculated  as 


(3-2) 
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where  x = sample  mean; 

s = sample  standard  deviation; 

la  = the  t statistic  with  (n- 1 ) degree  of  freedom; 


= 0.10  for  90%  confidence  level 
- number  of  samples. 


aid  be  identified. 


In  order  to  locate  the  bands  for  better  discriminating  between  sawgrass  and  cal 
the  spectrum  were  divided  into  five  narrow  bands,  namely  blue  (425-475  nm).  green 
(525-575  nm),  red  (625-675  nm),  near  infrared  (725-800  nm),  and  first  infiared  recess 
(925-975  nm).  The  mean  spectral  difference  at  these  five  bands  was  computed  for  the 
original  reflectances  and  the  normalized  reflectances.  A test  statistic  (T.S.)  similar  to 
paired  1-test  was  used  to  determine  the  better  bands  as  follows: 


(3.3) 
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measured  spectrum. 

3.1.4  Comparison  of  Spectroradiometer  and  NO  A A AVHRR 

Spectral  reflectance  of  the  radiometer  in  the  red  and  near  infrared  bands  were 
extracted  to  compare  with  the  band  I (red)  and  band  2 (near  infrared)  reflectance  of  the 
NOAA  AVHRR  data.  In  order  to  match  the  red-band  range  of  the  NOAA  AVHRR,  red 
band  reflectance  of  the  spectroradiometer  was  obtained  by  averaging  the  reflectance  range 
from  0.58  to  0.68  pm.  And  similarily,  near  infrared  reflectance  was  obtained  by 
averaging  the  range  from  0.73  to  1 .05  pm  to  match  the  near  infrared  band  of  the  NOAA 
AVHRR.  Reflectances  of  sawgrass,  cattail  and  open  water  on  the  aforementioned  dates 
were  used  in  the  comparison.  Satellite  reflectances  were  acquired  from  the  calibrated 
bands  1 and  2 of  NOAA  AVHRR  images  at  the  Fort  Drum  sawgrass  and  cattail  marsh. 
Three  vegetation  indices,  IR/R  ratio,  NDVI,  and  SAVI,  were  also  generated  from  the  red 
and  near  infrared  reflectance. 

3.2  Optimal  Spatial  Resolution 

importantly,  a significant  amount  of  wetlands  which  creates  a unique  geographical  feature 
different  than  most  other  states.  The  study  of  the  optimal  spatial  resolution  in  south 


Florida  was  therefore  focused  on  these  three  major  land  use/cover  types.  Agriculture 
land  use  was  selected  fiom  the  sugarcane  and  vegetable  fields  in  the  Everglades 
Agricultural  Area.  Urban  land  cover  was  selected  from  the  commercial  and  residential 

National  Park.  In  addition,  an  area  mixed  with  all  three  types  of  land  cover  plus  part  of 
Lake  Okeechobee  and  pasture  land  west  of  Everglades  Agricultural  Area,  labeled 
mixture,  was  also  investigated  in  this  study.  Figure  3. 1 illustrates  the  location  of  the  four 
selected  land  use/cover  areas.  Letter  A represents  agricultural  area,  W,  wetlands;  U, 

3.2.1  Image  Resolution  Degradation 

A geo-rectified  Landsat  Thematic  Mapper  (TM)  image  taken  on  April  1 , 1994 
supplied  image  data  used  to  study  optimal  spatial  resolution.  The  spectral  bands  and 
ground  resolution  of  Landsat  TM  are  listed  in  Table  3.1. 

Table  3. 1 Spectral  band,  wavelength,  spectral  location,  and  ground  resolution  of 
Landsat  Thematic  Mapper  (Lillcsand  and  Kiefer,  1987) 

Band  Wavelength  Spectral  Ground 

firm) location  Resolution,  (m) 

1 0.45-0.52  Blue  30 

2 0.52-0.60  Green  30 

3 0.63-0.69  Red  30 

4 0.76-0.90  Near  Infrared  30 

5 1.55-1.75  Mid  Infrared  30 

6 10.4-12.5  Thermal  Infrared  120 

7 2.08-2.35  Mid  Infrared  30 


Location  of  the  four  selected  land  usc/covers:  letter  A stands  for 
agricultural  area  (sugarcane  and  vegetable  fields  in  Everglades 
Agricultural  Area),  W for  wetlands  (Everglades  National  Park  and  Water 
Conservation  Area),  U for  urban  (residential  area  in  west  Fort  Lauderdale), 
and  dashed  line  for  the  mixture  area. 


The  near  infrared  band  and  Ihe  thermal  band  of  the  TM  were  used  for  the  optical 
and  the  thermal  portion  of  the  spectrum  to  determine  optimal  spatial  resolution.  The  near 
infrared  band  was  selected  due  to  its  capability  to  highly  discriminate  between  vegetation 

The  sub-images  of  the  four  areas  of  interest:  urban,  agriculture,  wetland  and 
mixture,  were  cropped  from  the  original  30  m resolution  Landsat  TM  south  Florida 
images.  Each  sub-image  was  then  degraded  to  coarser  resolutions  by  averaging 
neighboring  pixels  with  a incremental  scale  of 2, 3, 4, 5, 6,  8, 10, 20, 30, 40, 50,  and  100 
to  create  image  resolutions  of  60, 90, 120, 150, 180, 240, 300, 600, 900, 1200, 1500  and 


3*2.2  Mean  Local  Coefficient  pf  Variance 

The  spatial  structure  of  images  is  primarily  related  to  the  relationship  between  the 
size  of  the  objects  in  the  scene  and  spatial  resolution  and.  when  selecting  an  appropriate 
scale  of  data,  the  spatial  structure  of  images  is  a determining  factor.  Woodcock  and 
Strahlcr  (1987)  developed  a measure  of  local  image  variance  to  aid  in  selecting  an 
appropriate  image  scale  for  a variety  of  land  cover  types.  However,  the  local  variance  for 
each  land  use  type  can  not  be  compared  to  each  other  due  to  the  high  variability  of  the 
calculated  variance.  Furthermore,  the  variance  calculated  from  one  satellite  can  not  be 
compared  with  data  from  other  satellites  or  the  chronologic  data  of  the  same  satellite  due 


67 

effects.  In  this  study,  a new  image  quantity.  Mean  Local  Coefficient  of  Variance 
(MLCV),  was  introduced  and  rated  to  study  the  appropriate  spatial  resolution.  The 
MLCV  is  measured  as  the  mean  value  of  the  coefficient  of  variance  of  a 3«3  pixel 
window  moving  across  the  image  (Figure  3,2).  The  coefficient  of  variance  (CV)  of  the 
nine  pixel  values  in  the  3x3  array  is  computed  and  the  mean  of  these  CV  over  the  entire 
image  is  treated  as  an  indication  of  the  local  variability  of  the  image.  The  reason  for  such 
an  indication  is  that  if  the  spatial  resolution  is  considerably  finer  than  the  objects  in  the 
scene,  most  of  the  measurements  in  the  image  will  be  highly  correlated  with  their 
neighbors  and  a measure  of  local  CV  will  be  low.  If  the  resolution  is  reaching  an 
approximate  size  of  the  objects  on  the  ground,  the  likelihood  of  neighboring  pixels  being 
similar  decreases  and  the  local  CV  increases.  Graphs  of  local  coefficient  of  variance  as  a 
function  of  spatial  resolution  can  then  be  used  to  measure  spatial  structure  of  images. 

3.3  NQAA  AVHRR  Data  Calibration  and  Image  Processing 
3.3.1  NQAA  T1ROS-N  Satellites 

Since  1 960,  the  National  Aeronautics  and  Space  Administration  (NASA)  and 
National  Oceanic  and  Atmospheric  Administration  (NOAA)  have  launched  and  operated 
a series  of  Polar  Orbiting  Environmental  Satellites  (POES).  This  series  of  satellites 
commenced  with  TIROS-N,  continued  with  NOAA-6  to  NOAA-14,  and  is  the  longest- 
lived  and  most  influential  series  of  Earth  observing  satellites  ever  launched  (Hastings  and 
Emery,  1992).  The  NOAA  satellites  were  initially  designed  to  observe  Earth's  weather 
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Hie  diagram  of  MLCV.  The  MLCV  is  measured  as  the  mean  value 
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from  the  Earth's  surface,  scientists  have  found  that  the  NOAA  satellites  have  a wide  range 
of  applications,  not  only  in  meteorology  and  climatology,  but  also  in  land  surface  studies. 

The  orbits  of  the  NOAA  satellites  arc  near-polar,  sun-synchronous,  and  with  an 
orbital  period  of  approximately  1 02  minutes,  produce  14.1  orbits  per  day.  Since  the 
number  of  orbits  per  day  is  not  an  integer,  the  exact  sub-orbital  tracks  do  not  repeat  on  a 
daily  basis,  however,  the  local  solar  time  of  the  satellite's  passage  is  essentially  the  same 
for  a given  latitude.  The  primary  environmental  sensor  for  the  NOAA  satellite  series  is 
the  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  which  is  a five  channel 
scanning  radiometer  sensing  in  the  visible,  near  infrared,  and  thermal  infrared  window 
regions.  The  channel  width  and  features  arc  tabulated  in  Table  3.2. 

Table  3.2  NOAA/AVHRR  spectral  ranges. 

Channel  Spectral  ranges  Description 


1 0.58-0.68  pm 

2 0.725-1.05  pm 

3 3.55-3.92  pm 


4 10.3-1 1.3  pm 

5 1 1.5-12.5  pm 


visible  (red) 

thermal  infrared 
thermal  infrared 


AVHRR  channels  I and  2 were  designed  to  discern  clouds,  land-water 


boundarie 
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terrestrial  vegetation  employing  the  computation  of  the  vegetation  indices.  Channels  3, 4, 
and  5 were  designed  to  be  used  to  measure  the  temperature  of  clouds  and  the  sea  surface 
and  for  night  time  cloud  mapping.  Current  applications  have  far  exceeded  these  original 
objectives. 

The  AVHRR  is  a scanning  radiometer,  utilizing  a rotating  scan  mirror  at  a rate  of 
360  rpm.  The  instantaneous  Field  of  View  (1FOV)  for  all  channels  is  approximately  1 .4 

altitude  of 833  km.  The  instrument  has  been  designed  in  such  a way  that  the  [FOV  of  the 
channels  can  be  made  coincident  within  ±0. 1 milliradians.  This  instantaneous  field  of 

lines  to  be  contiguous  at  the  sub-point  (Schwulb,  1978).  Figures  3.3  through  3.7  provide 
the  spectral  response  curves  for  the  NOAA  14  satellite.  Note  that  the  AVHRR  sensors 
have  near  uniform  distributions  over  their  sensitive  ranges.  The  analog  data  output  from 
the  sensors  is  digitized  to  a 10-bit  precision  digital  number  on  board  the  satellite  at  a rate 
of  39.936  samples  per  second  per  channel.  A total  of  2048  samples  is  obtained  per 


rustics  of  NOAA-AVHRR  arc  listed 
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Figure  3.4  Spectral  response  curve  for  NOAA-14  AVHRR  channel  2 (adapted  from 

Kidwell,  1991). 
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Figure  3.5  Spectral  respona 
Kidwell,  1991). 
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Figure  3.6 


Spectral  rcsponsi 
Kidwell,  1991). 


I AVHRR  channel  4 (adapted  I 
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NORA- 14  AVHRR  channel  5. 


Figure  3.7  Spectral  response  curve  for  NOAA- 14  AVHRR  channel  5 (adapted  from 

Kidwell,  1991). 


Table  3.3  NOAA/AVHRR  orbit,  ground  i 


characteristics. 


Orbit  height 
Orbit  inclination 

Scan  angle 
Ground  coverage 


S33  km 
98.8  degrees 
101.38  min 

55.4  degrees  on  cither  side  of  nadir 
2800  km 


Image  width 


Calibration 


1 . 1 km  (nadir),  2.4  km  (max.  off-angle  along  track) 

6.9  km  (max.  off-angle  across  track) 

2048  pixels 

1 0-bit  binary  (Digital  Number  0 to  1023) 
space  (zero  albedo)  for  Channels  1 and  2 

blackbody  surface  inside  spacecraft  (290  K)  and  space  (0 
K)  for  Channels  3, 4,  and  5 


Based  on  their  cloudlessncss,  five  day-and-night  pairs  of  NOAA  AVHRR  images 
were  used  in  this  study.  The  images  were  obtained  from  the  NOAA  satellite  active 


; (SAA)  using  i 


File  Transfer  Protocol  (FTP)  I 


sfer  the  images  through  the 
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Tabic  3.4  Dale,  overpassing  lime  and  orbil  number  of  the  NOAA  AVHRR  data. 


24-1996  NOAA  14 


19:10:46 

07:58:36 


3.3.2  Processing  nfNOAA  AVHRR  Dnl.-i 

Information  collccicd  via  satellite  sensors  generates  large  amounts  of  data  and  is 
transmitted  in  complex  formats.  Raw  AVHRR  data  come  with  auxiliary  data  such  as 
calibration  coefficients,  solar  zenith  angles,  earth  locations,  and  images  in  packed  1 0-bit 
radiometric  resolution.  The  raw  images  contain  both  radiometric  and  geometric  errors 
which  have  to  be  removed  in  order  to  quantitatively  analyze  and  overlay  other  data  sets, 
Howt.c.,  con. cn„o, .ally  available  software  packages  were  inadequate  for  the  processing 
of  AVHRR  data  at  the  time  the  images  were  acquired.  Research  was  conducted  to 
address  the  AVHRR  preprocessing  problem  and  a set  of  programs  written  in  C language 
were  developed  to  facilitate  the  data  and  image  preprocessing  procedures.  The  computer 
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programs  include  procedures  lo  extract  raw  image  and  auxiliary  data,  correcl  radiometric 
errors  using  the  calibration  coefficients,  correct  geometric  errors  using  the  earth  locations, 
output  solar  zenith  angles  and  telemetry  data.  This  program  also  included  subroutines  to 
rotate  the  upside-down  image  from  descending  ground  tracks  and  subroutines  to  save 
images  in  raw,  ELAS  (Graham  et  al„  1985)  and  ERDAS  (ERDAS,  1995)  file  formats. 
The  description  of  subroutines  and  program  codes  are  listed  in  Appendix  3.  A schematic 
diagram  of  the  preprocessing  of  AVHRR  HRPT  data  is  illustrated  in  Figure  3.8. 

3.3.2.  l image  and  daw  extraction 

The  AVHRR/HRPT  tape  obtained  from  NOAA/NESS  contained  about  1 1 
minutes  of  data  covering  a ground  swath  of  4500  km  (along  track)  times  2700  km  (cross 
track).  The  recorded  satellite  pass  consisted  of  about  4000  scanlines  and  each  scanline 

two  physical  records  per  scan-line.  The  data  was  digitized  to  10-bit  precision  and 
ordered  in  band  interleaved  by  element  format  (scan  clement  1 band  1,2, 3, 4. 5,  scan 
clement  2 band  1 , 2, 3, 4, 5,  etc).  The  solar  zenith  angle  and  earth  location  data  (latitude 
and  longitude)  ate  sampled  every  40  points  starting  at  point  25  of  each  scan-line  and  were 
embedded  in  the  data  record.  The  format  of  the  HRPT  physical  record  is  shown  in  Table 
3.5. 


:igure  3.8  Schematic  diagram  of  NOAA/AVHRR  da 


Table  3.5  Fonnat  of  HRPT  data  physical  i 


Byte  number 

Content  in  HRPT  data  file 

1-2 

2 bytes  of  scan  line  number 

3-8 

7-bit  year,  9-bit  day  of  year 
27-bit  millisecond  GMT  time  of  day 

9-12 

4 bytes  of  quality  indicators 

13-52 

4-byte  slope,  intercept  coefficient  for  each  channel 

53 

number  of  meaningful  zenith  angle  and  earth  location 
points  appended  to  scan 

54-104 

5 1 bytes  of  solar  zenith  angles" 

105-308 

204  bytes  of  earth  locations" 

309-448 

140  bytes  of  telemetry  data 

449-14104 

13656  bytes  of  LAC/HRPT  image  data 

-14800  696  bytes  unused 


. physical  record  (14800  bytes)  r 


aprised  of  two  logical  records  (7400  bytes). 


! data  in  AVHRR  tape  are  packed t 


as  three  10-bit  samples  in  4 byte 
(8bit/byte)  groups,  right-justified,  with  the  first  two  bits  zero,  as  illustrated  below: 


The  packed  format  was  intended  to  arrange  the  image  data  with  as  much  spectral 
information  and  as  few  spaces  between  the  image  elements  as  possible.  However,  this 

to  extract.  The  NOAA  AVHRR  processing  software  written  by  me  was  able  to  unpack 
the  data  and  reformat  it  so  that  the  individual  10  bit  samples  were  stored  in  consecutive 
pairs  of  bytes  (i.e.  1 0 bit  data  to  a 1 6 bit  word)  via  bit-field  data  structure  and 
manipulation.  The  data  from  the  same  channels  were  also  rearranged  and  stored  in 
sequence  so  that  the  original  five-channel  interleaved  format  becomes  a consecutive 
channel  for  each  scan-line  (B1L  format).  The  physical  data  header  which  contains 

read  by  the  unpacking  procedure  and  stored  in  a separate  header  file  for  further  analysis. 

3.3.2.2  Geometric  rectification  and  geo-refercncine 

Satellite  images  arc  subjected  to  different  geometric  deformations  due  to  the 
Earth's  rotation,  oblatencss,  curvature,  speed,  attitude,  and  altitude  variation  of  the 
satellite.  The  scan  skew  and  the  projection  of  spherical  surface  on  a flat  image  also 


contributes  to  geometric  errors.  Those  deformations,  if  not  properly  accounted  for,  will 
prevent  meaningful  comparison  among  data  from  ground-truthing  and  images  acquired 
from  different  sources,  and  will  also  hinder  further  applications  of  the  geographical 
information  system. 

Geometric  rectification  of  the  AVHRR  images  were  performed  with  the  following 
four  major  steps:  ( 1 ) establishing  a set  of  ground  control  points  (GCP).  GCPs  are  points 

map.  Depending  on  the  image  resolution,  the  GCPs  may  be  a lake,  a highway 
intersection,  or  a permanent  land  feature.  For  NOAA  AVHRR  images  of  Florida, 
dependable  ground  control  points  are  the  distinct  shoreline  features  and  permanent  water 

adjustment  of  both  x and  y distance  on  a plane.  The  geodetic  coordinates 
(latitude/longitude)  is  not  a plane  projection,  and  when  applied  in  a small  area,  the 
discrepancy  may  not  be  noticeable.  However,  it  is  important  to  use  a plane  projection  in 
large  coverage.  The  HRPT  processing  program  is  capable  of  converting  the  geodetic 
coordinates  to  UTM  coordinates;  (3)  establishing  mapping  equation.  Two  first  order 
polynomial  equations  were  used  to  map  the  GCP  coordinates  (UTM  easting  and  northing) 
to  the  corresponding  AVHRR  image  coordinates  (line/column)  in  the  program.  The 
image  distortion  caused  by  rotation,  skew,  translation  and  scale  change  was  eliminated  by 
the  first  order  transformation.  The  mopping  equation  is  as  follows: 


(3.4) 
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v * + b,x  + bjy  (3.5) 

where  u is  the  image  line,  v ihe  image  element,  x the  UTM  easting,  y the  UTM  northing, 
a's  and  b's  are  the  mapping  coefficients  to  be  determined.  This  transformation  is  widdy 
used  and  is  sufficient  to  fit  satellite  images  on  UTM  maps  when  the  terrain  docs  not 
undulate  significantly.  Least  squares  technique  was  used  in  the  program  to  estimate  the 
transformation  coefficients.  The  two  sum  of  square  errors 

E[u-<a0  + a1x  + aJy)]!  (3.6) 

E[v-(b0  + b,x  + b!y))!  (3.7) 

solved  simultaneously;  and  (4)  resampling.  Resampling  is  a process  of  forming  a new, 
rectified  image  by  choosing  values  ftom  the  original  image.  The  nearest-neighbor 
resampling  scheme  was  used  in  the  program. 

from  satellite  observations  (Tarpley,  1988;  Saunders.  1988).  For  either  surface 
temperature  or  vegetation  index  studies,  cloud  contaminated  pixels  should  be  excluded 
from  the  analysis.  Clouds  generally  have  a higher  reflectance  in  short-wave  channels  and 
a lower  thermal  emission  than  the  earth's  surface  in  thermal  channels.  Therefore  a 
certain  digital  number  threshold  can  be  set  by  manually  observing  pixels  containing 


Rundquist.  1994).  This  procedure  is  most  effective  in  the  thermal  spectrum  as  cirrus 
clouds  arc  often  invisible  in  the  reflective  spectral  range  (Price.  1990).  Clouds  in  the 
night-time  images  are  even  more  difficult  to  detect,  not  only  due  to  the  unavailable  visible 

day-time  images,  band  3 of  the  AVHRR  (3.7  pm)  night  time  images  appear  broadly 
similar  to  the  corresponding  1 1 pm  images.  However,  Hunt  (1973)  reported  that  the 
emissivity  of  low  cloud  or  fog  at  3.7  pm  wavelength  is  significantly  less  than  at  12  pm 
producing  a significant  difference  in  measured  brightness  temperature  between  the  two 
bands  which  is  normally  much  greater  than  the  difference  caused  by  atmospheric 
absorption.  Since  both  day  and  night  time  AVHRR  images  are  to  be  used,  the  cloud 
cover  condition  was  determined  in  two  steps  using  both  temperature  sensed  from  the 
thermal  infrared  bands  3 and  5,  and  from  visible  band  1.  The  first  step  uses  the 
temperature  difference  detected  in  3,7  and  12  pm  (band  3 and  band  5 of  AVHRR  data). 

If  the  brightness  temperature  difference  (TirT,  J is  greater  than  5°C,  then  the  pixel  was 
flagged  as  cloud  contaminated.  The  5°C  was  suggested  by  Saunders  (1986)  and  proved 
appropriate  in  Florida  (Shih  and  Tan,  1997).  The  second  step  was  to  examine  whether 
the  HRPT  image  was  discernible  using  the  visible  band  I to  supplement  the  cloud 
identification.  The  cloud  cover  areas  were  then  delineated  and  eliminated  in  this  study. 

3.3.3  AVHRR  Data  Calibration 

meaningful  surface  characteristics.  AVHRR  visible  data  values  (Channels  1 and  2)  may 
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be  converted  lo  albedos  and  AVHRR  thermal  data  values  (Channels  3 and  4,  and  5)  may 


AVHRR  visible  channels  (cha 
conversion  between  the  digital  numbe 
albedo  measured  by  the  sensor  channe 


el  1 and  channel  2)  arc  designed  to  provide 
and  surface  reflectivity  (albedo).  The  percent 
is  computed  as  a function  of  the  input  data  value 


A|=  SjC  + 1,  (3.8) 

where  A,  is  the  percent  albedo  measured  by  channel  i,  C is  the  input  data  value,  and  S, 
and  I,  are  the  scaled  slope  and  intercept  values,  respectively.  The  scaled  slopes  and 
intercepts  were  calculated  front  the  calibration  coefficients  extracted  from  HRPT  data  file 
(Kidwcll  1991). 


The  thermal  channels  (channel  3, 4,  and  5)  are  designed  and  calibrated  before 
using  the  inverse  of  Plank's  radiation  equation: 


(3.9) 
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where  T = radiant  temperature  (K)  for  the  energy  value  E, 

C,=  constant  (1.191 0695*  1 0J  milliwatts/m;-steradian-cm"‘), 

C ,=  constant  (1.438833  ctn-K), 

E = energy  measured  by  the  sensor  (channel  I),  and 
E,  = S,C*  I, 

'Hie  temperature  obtained  by  this  procedure  is  radiant  temperature,  which  needs 
further  correction  for  atmospheric  attenuation  and  surface  cmissivity.  The  atmospheric 
correction  of  brightness  temperature  was  performed  by  the  water-body  calibration 
technique.  Hourly  water  surface  temperature  at  four  stations  in  Lake  Okeechobee  were 
obtained  from  the  South  Florida  Water  Management  District.  The  water  surface 
temperatures  were  used  since  they  provided  large  homogeneous  surface  temperature  with 

The  atmospheric-corrccted  radiant  temperature  is  the  actual  surface  temperature 
only  when  the  cmissivity  is  equal  to  1.0.  For  most  surfaces,  where  the  emissivity  is  not 
equal  to  1.0,  a further  adjustment  needs  to  be  undertaken.  Emissivity  adjustments  were 
performed  using  the  Stcfon-Boltzmann  relationship  as  follow: 


Table  3.6 


of  satellite  surfai 


temperature. 


Station  ID  CO  Station  Description 


Lat  Long. 


L002+T0  OKE 
L005+TO  GLA 
L006+T0  PAL 
LZ40+TO  PAL 


Lake  Okeechobee  Tower  North  (#2) 
Lake  Okeechobee  Tower  West  (#5) 
Lake  Okeechobee  Tower  South  (#6) 
LZ40  Weather  Station  on  Lake  Okcc, 


270507  804715 
265734  805838 
264920  804659 
265405  804721 


found  in  literature  (Salisbury  and  D'Aria,  1992;  Lillsand  and  Kiefer,  1987;  Rees.  1990; 
Griggs,  1968).  The  emissivity  values  of  Florida  land  use/covcr  used  in  surface 
temperature  adjustment  were  listed  in  Table  3.7.  These  values  are  generally  higher  than 
values  obtained  from  laboratory  measurements  due  to  the  large  area  coverage  and  they  are 


Emissiv 


; for  the  land  use/covcr  classes  in  Florida. 


21 

22 


Coastal  strand 


Shrub  and  brushland 


0.94 

0.97 

0.96 

0.92 


0.97 

0.97 

0.97 

0.99 

0.99 


0.99 

0.98 

0.97 

0.97 


Land  use/covcr  types  were  adapted  from  the  classes  in  Florida  Land  Cover  (Florida 
Game  and  Freshwater  Fish  Commission,  1992) 


J.4  Regional  ET  Estimation  Mode 


ion  wilh  Surface  Temperature 

:n  latent  heat  flux  IE  and  surface  temperature  can  be  expressed 


LE  - RJ.I  -o)  . i(Rl-oT?)  - G ♦ P‘^r°~r«)  (3.11) 

R,  = incoming  short-wave  radiation  flux  (W  m';); 

e = surface  emissivity  coefficient; 

Rl  = long-wave  atmospheric  radiation  flux  (W  n>'!) 
a = Boltzmann  constant  (5.67  x 10-8  Wm‘!  KJ); 
p - density  of  the  moist  air  (kg  m J); 

Cp  = specific  heat  of  moist  air(J  kg'1  K1); 

T.  -air  temperature  (K); 

T,  - surface  temperature  (K);  and 

G = ground  heat  flux  into  the  soil  (W  m'!). 

It  con  be  seen  from  Equation  3.11  that  the  latent  heat  flux  depends  on  a number  of 
and  surface  parameteis.  Surface  temperature,  T„  can  be  remotely  sensed 
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by  thermal  infrared  satellite  observations.  The  diffusion  resistance,  r„,  depends  on  wind 
velocity,  v,  and  roughness  length  of  the  evaporating  surface,  z„.  Asa  first  approximation, 
/fj,  Rl  and  v can  be  taken  as  constants  over  a regional  area.  The  surface  parameters 
albedo  o,  emissivity  e and  roughness  length  can  be  determined  in  the  field  or  derived 
Horn  visible-image  data. 

This  energy  balance  equation  contains  many  parameters,  several  of  which  arc 
insensitive  to  ET  estimation.  A simplified  form  was  first  proposed  by  Jackson  ct  al. 
(1977)  and  further  adopted  by  Scguin  and  Itier  (1983),  Nicucnhuis  et  al.  (1985)  and 
Seguin  et  al.  (1989).  Using  this  approach  the  daily  ET  value,  ETj,  is  given  as  a (unction 
of  the  instantaneous  value  of  the  difference  between  the  surface  temperature  and  the  air 
temperature  (both  measured  at  noon),  and  of  the  daily  value  for  the  net  radiation,  K„(: 

- BTj  = A +B{T,-TJ  (3.12) 

where  R^ETj  are  in  energy  unit  (Wm  J); 

T„  Ta  in  temperature  unit  (°C);  and 
A (W  m !),  and  B (W  m‘;  °C')  arc  constants. 

Theoretical  justification  for  Equation  3.12  can  be  found  by  reference  to  the 
surface  energy  balance  integrated  over  a one-day  period.  The  G value  as  given  in 
Equation  3. 1 1 can  be  neglected  over  a 24-hour  period.  The  constant  B is  an  average  of 
bulk  conductance  for  the  daily  integrated  sensible-heat  flux.  The  surface  temperature  vs. 
air  temperature  variation  0.-TJ  is  defined  as  SA TV  for  this  study.  According  to  the 


literature,  the  constants  A and  B vary  with  location.  Jackson  et  al.  (1977)  obtained  A = 0 
and  = 0.54  for  the  Phoenix,  AZ  area,  while  Seguin  and  Iticr  ( 1 983)  obtained  A=  I and 
B = 0.25  for  the  area  of  La  Crau,  France. 

Ground-based  meteorological  data  were  accessed  from  the  hydrological  database 
for  the  34  stations  of  the  SFWMD.  The  data  were  used  for  calibrating  Equation  3.12, 
including  air  temperature,  net  radiation  and  the  parameters  used  for  potential  ET 
calculations.  Air  temperature  and  net  radiation  values  at  the  ground  stations  were 
interpolated  to  form  temperature  and  radiation  map  layers  with  a grid  si2e  equal  to  the 
NOAA  satellite  ground  resolution  (I  km).  A kriging  method  with  a linear  variogram  was 
used  to  formulate  the  grid  using  the  SURFER  software  (Golden  Software,  inc..  Golden, 
CO)  for  both  air  temperature  and  net  radiation.  Both  air  temperature  and  radiation  maps, 

registered  into  the  IDRISI  raster  geographic  information  system  (Eastman,  1995)  for 
further  calculation.  Potential  ET  at  the  34  ground  stations  was  calculated  based  on  the 
modified  Penman  combination  equation  (Jensen  et  al.,  1990)  and  stored  in  the  GIS 
database  as  point  coverage.  It  should  be  noted  that  both  air  temperature  and  net  radiation 
have  been  interpolated  to  all  points  throughout  the  entire  study  area.  The  ET  then  could 
be  estimated  based  on  the  land-surface  characteristics,  using  satellite-derived  surface 

The  relationship  of  the  difference  between  net  radiation  and  potential  ET.  (R.-ET), 
and  the  (TrTJ  values  in  Equation  3.12  was  analyzed  using  a regression  technique 
involving  data  from  the  34  ground  stations.  The  established  relation  could  then  used  to 


calculate  ET  values  al  all  grid  points  using  surface  temperature,  air  temperature,  and  net 
radiation. 

3.4.2  Relationship  of  Vegetation  Index  and  Surface  Temperature 

Surface  temperature  is  primarily  dependent  on  the  thermal  inertia  of  the  soil  after 
meteorological  inputs  to  the  soil  surface  have  been  accounted  for  (Engman  and  Gurney, 
1991).  Thermal  inertia  is  a measure  of  the  thermal  resistance  of  a material  to  temperature 
changes.  Some  materials,  e.g.  dry  soil,  have  low  thermal  inertias  and  exhibit  wide  diurnal 
ranges  of  temperature.  Others,  e.g.  wet  soils,  have  high  thermal  inertias  and  then- 
observed  diurnal  temperature  ranges  are  relatively  low.  The  orbit  of  the  NOAA  satellite 
is  designed  to  overpass  at  14:00  and  02:00  hours  when  surface  temperatures  are  near  the 
daily  maximum  and  minimum,  respectively.  The  diurnal  difference  of  these  two 
temperatures  has  been  proven  useful  in  determining  thermal  inertia  and  soil  moisture 
(Van  de  Griend  et  al.,  1985).  Several  studies  (Jackson  et  al„  1977;  Jackson,  1982; 
Schmugge,  1988)  verified  that,  over  bare  soil,  variations  in  surface  temperature  tend  to  be 
highly  correlated  with  variations  in  surface  water  content. 

The  NDVI-DSTV  triangle  approach  (Figure  3.9)  is  based  on  the  hypothesis  that  a 

surface  temperature  variation.  The  top  of  the  triangle  corresponds  to  full  vegetation- 
covered  areas.  The  base  of  the  triangle  corresponds  to  areas  with  low  vegetation  index, 
either  bare  soil,  barren  or  urban  area.  In  between  the  triangle  top  and  base,  the  major 
portion  of  the  triangle  corresponds  to  land  with  various  fractional  vegetation  t 


r.  The 


spectivcly  to  high  i 


Figure  3.9  The  hypothetical  triangle  shape  that  would  result  from  the  relation 

between  vegetation  index  and  diurnal  surface  temperature  difference. 


The  two  base  points  also  correspond  to  the  water  content  of  the  surface.  This  is  due  to 
the  high  thermal  inertia  of  water  molecules  which  absorb  high  amounts  of  solar  energy 
with  little  temperature  change.  The  lower  right  comer  of  the  NDVI-DSTV  triangle,  with 
low  vegetation  index  and  high  surface  temperature,  is  the  location  of  diy  soil  land  cover 
while  the  lower  left  point  is  the  location  of  moist  soil. 

3.4.2. 1 NDVI-DSTV  Irianale  from  field  measurement 

The  NDVI-DSTV  triangle  concept  was  evaluated  using  the  field  measurements  of 
surface  reflectance  and  surface  temperature.  Several  ground  covets,  including  asphalt, 

water  with  70%  vegetation  cover,  and  water  with  0%  vegetation  cover  were  measured 
from  ground  level  to  evaluate  the  NDVI-DSTV  triangle  relationship. 

Field  surface  temperatures  were  collected  on  several  days  coincident  with  NOAA 
satellite  overpass  times  (0230h  and  I430h).  Surface  radiative  temperatures  of  selected 
land  covers  were  measured  with  an  Everest  Inlcrsciencc  (Model  210)  hand-held  infrared 
thermometer  set  for  a emissivity  of  0.97.  Surface  radiative  temperatures  were  taken  in  the 
four  directions  and  averages  were  computed.  Field  spectral  reflectances  were  measured 
with  GER-1500  hand-held  spectroradiomcter  following  the  procedures  stated  in  section 


3.4.22  NDVI-DSTV  triangle  from  satellite  data 

The  NDVI-DSTV  triangle  relationship  from  NOAA  AVHRR  was  verified  by 
plotting  the  scattergram  of  NDVI  and  DSTV  for  seven  major  land  cover  types  in  south 


rf  each  type  i 
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Florida.  The  land  cover  types  and  the  locations  o 
listed  in  Table  3.8.  Locations  of  the  land  cover  were  obtained  from  the  classified  Landsat 
TM  image  of  the  Florida  Land  Cover  (Florida  Game  and  Fresh  Water  Fish  Commission. 
1992). 

Table  3.8  The  seven  land  cover  types  and  locations  of  classes  used  in  the  scattcrgram 
Name  Land  cover  type  Location  of  class  used  in  scattcrgram 


Agriculture 

Pine 


Urban 


grass  land  and  agriculture 


hardwood  forest 

diy  prairie 
cypress  swamp 
barren  and  urban  land 


west  of  J.W.  Corbett  WMA 
cast  of  Naples.  Collier  county 
Everglades  National  Park 
cast  Dc  Soto  county 
Big  Cypress  National  Preserve 


3.4.3  ET  Estimation  from  NPVI-DSTV  Triangle 

surface.  Meteorological  parameters  can  be  obtained  from  ground  stations  and  their 
values,  generally  can  be  interpolated  to  sites  without  measurements.  On  the  other  hand, 
surface  parameters  are  difficult  to  measure  over  large  areas  and  can  not  be  interpolated  to 


: the  maj 


affect  ET  rate.  Presence  of  green  vegetation  is  a major  determinant  of  ET  from  the  land 
surface  due  to  enhanced  surface  roughness  increasing  turbulent  exchange  of  water  vapor 
and  plant  roots  extracting  water  from  the  soil  more  rapidly  than  the  water  can  diffuse  to 
the  soil  surface  (Smith  and  Choudhury,  1991).  The  vegetation  characteristics  which  are 
related  to  the  ET  rate  include  type  of  vegetation,  plant  density,  and  vigor.  A quantitative 
description  of  vegetation  and  soil  characteristics,  however,  is  difficult  to  acquire  at 
regional  scale.  Data  observed  by  satellite  remote  sensing  usually  contains  areas  with 
various  vegetation  density  and  soil  moisture  content  Each  pixel  of  the  satellite  image  can 
be  assumed  a combination  of  the  responses  from  these  sub-pixel  areas.  The  fraction  of 

relationships  with  NDVI  and  DSTV,  respectively. 

Considering  the  large  ground  cover  of  each  AVHRR  pixel  (~  1 km1),  it  is  assumed 
that  there  is  non-homogeneity  within  each  area  corresponding  to  a satellite  measurement 
In  Florida,  each  area  can  be  considered  as  a mixture  of  vegetation,  wet  soil,  and  dry  soil. 
The  surface  temperature  and  vegetation  index  sensed  by  satellite  are  also  considered  to  be 

dry  soil.  Letting  subscript  v correspond  to  vegetation,  w to  wet  soil,  and  d to  dry  soil,  the 

T,-*,T,  + *iT.  + w,Tt  (3.13) 

+ w,  NDVI , + v>,  NDVIj 


NDVI,  - iv,  NDVI, 


(3.14) 


where  w„  w„  and  w,  are  weighting  factors  corresponding  to  the  vegetation,  moist  soil  and 
dry  soil,  respectively. 

Similarly,  day-time  and  night-time  temperature  can  be  written  as 

+ (3.15) 

*1  - »r  T..  + Wj  T,  ^ + w,  Tj  ^ (3.16) 

Their  difference,  DSTV,  can  be  written  as 

DSTV,  = W,  DSTV,  + w,  DSTV.  + w,  DSTV, 

Combining  Equations  (3.14)  and  (3.17),  adding  the  third  equation,  w,  + w;  + w, 
denoting  T as  DSTV  and  V as  NDVI,  the  weighting  factors  can  be  solved  as 

. T,(V.-VJ)  + V.(Z-T.)+T_V,-ZV_ 

" TAK-V.)*VJ(T.-T,)  + T,V,-T.V. 

w _ nV*  ~ V.)+V,(T. -Tj)+  ZV,  - T.V, 

1 TAK-vj+v^r.  - T,)+T,V.-T.V, 

w _ T,(K -K)+v.(T.-r.)+ T.v.  - t_v_ 

' ~ TAK -y.)+vd(T.-r,)+  t,v.-t.v. 

The  calculation  of  the  weighting  factors  depends  on  the  vegetation  index  and 
diurnal  surface  temperature  variation  of  the  three  comer  types  and  this  procedure  can  be 
easily  written  into  programs  and  determined  efficiently. 


(3.18) 

(3.19) 
(3-20) 


(3.17) 
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DSTV  (°C) 


Figure  3.10  Relationship  of  vegetation  and  moisti 
comer  points  (vegetation,  wet  land,  a 


lure  coefficient  (VMC)  from  the  three 
rnd  dry  land)  of  the  NDVl-DSTV 
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of  dry  land.  The  VMC  for  the  vegetation,  wet  land  and  dry  land  are  listed  in  Table  3.9. 
Monthly  VMC  of  vegetation  are  determined  by  the  ratio  of  sugarcane  ET  and  potential 
ET  reported  by  Shihand  Gascho  (1980)  and  Jones  et  al.  (1984).  VMC  of  wet  land  was 
assumed  a value  of  0.95  indicating  the  ET  is  close  to  the  potential  ET  (ET„).  VMC  of  dry 
land  was  assumed  a value  of  0.05  indicating  the  ET  is  much  smaller  than  the  potential  ET 
due  to  high  percentage  of  buildings  and  pavement  in  urban  areas. 

Table  3.9  Vegetation  Moisture  Coefficient  (VMC)  for  agriculture,  wet  land  and  dry 


Jan  Feb  Mar  Apr  May  Jun  Jul  Aug  Sep  Oct  Nov  Dec 

Vegetation'  0.53  0.32  0.52  0.56  0.73  0.94  1.02  1.08  0.95  1.14  1.01  1.00 

Wet  Land’  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95 

Dry  Land’  0.05  0.05  0.05  0.05  0.05  0.05  0.05  0.05  0.05  0.05  0.05  0.05 

1 Crop  coefficients  of  sugarcane  from  Shih  and  Gascho  ( 1 980)  and  Jones  et  al.  ( 1 984). 

1 assumed  wet  land  surface  ET  close  to  ET„ 

1 assumed  dry  land  surface  ET  much  smaller  than  ET0 

Potential  ET  values  were  calculated  with  meteorological  parameters  from  the 
ground  network  and  interpolated  to  all  areas  in  south  Florida.  The  potential  ET  can  be 
interpolated  due  the  fact  that  they  arc  determined  from  meteorological  data  which  are 
considered  less  dependent  of  surface  condition.  Actual  ET  values  were  computed  as  the 
result  of  multiplying  potential  ET  by  the  VMC  derived  from  satellite  data. 


ET,  - VMC,  x ETft , 


(3.22) 


where  ET,  ■ actual  ET  for  pixel  i, 

VMC,  = vegetation  and  moisture  coefficient  for  pixel  i. 

ET0  , = potential  ET  calculated  by  traditional  method  and  interpolated  to  pixel 


3.4.4  Geographic  Information  System 

data,  the  geographic  information  system  (G1S)  was  used  to  manage  the  large  amount  of 
satellite  and  ground  data.  The  GIS  was  implemented  by  the  ERDAS  Imagine  software 
package  (ERDAS,  1995).  Satellite  derived  day  and  night  lime  surface  temperatures, 

index  were  stored  in  the  GIS.  The  weighting  factor  layers,  (w,,  tvj,  tv,),  vegetation  and 
moisture  coefficient  (VMC),  interpolated  potential  ET  and  actual  ET  were  calculated  in 
the  GIS  system.  Regional  ET  was  computed  by  summing  up  ET  at  the  pixels  within  the 


3.5  Model  Verification 


Verification  of  satellite  derived  regional  ET  is  difficult  due  to  lack  of  ET  records 
representing  areas  large  enough  to  compare  with  satellite  derived  values.  Traditional  ET 
estimation  methods,  regardless  how  complex,  still  represent  a point  estimation.  In  fact, 
the  remote  sensing  ET  estimation  method  in  this  study  still  utilized  traditional  ET 


i for  potential  ET  estimation.  Only  1 


complex,  surface-dependent  crop  coefficients  (vegetation  and  moisture  index)  which  are 
difficult  to  measure  by  traditional  methods,  were  determined  by  remote  sensing  methods. 
Based  on  these  limitations,  verification  of  the  model  was  twofold:  I)  compare  point  ET  at 
weatherstations  by  tradition  ET  estimation  method  and  by  remote  sensing  method;  2) 
compare  ET  at  wetland  lysimeter  sites.  The  iysimeters  were  located  within  a large  marsh 

3.5.1  Verification  with  Traditional  ET  Methods 

The  meteorological  data  from  26  stations,  including  pan  evaporation,  solar 
radiation,  air  temperature,  relative  humidity,  atmospheric  pressure,  and  wind  speed,  were 
obtained  from  the  South  Florida  Water  Management  District  (SFWMD).  Using  the 
meteorological  data  from  the  ground  stations,  the  ET  rates  estimated  by  remote  sensing 
method  were  compared  with  those  by  traditional  ET  estimation  methods.  Potential  ET  in 
the  remote  sensing  method  was  estimated  by  the  Penman  method.  ET  estimated  by  the 
Penman  method  was  used  to  compare  the  results  from  the  remote  sensing  method.  The 

in  the  Florida  land  cover  map  (Florida  Game  and  Fresh  Water  Fish  Commission,  1 992). 
The  land  cover  type  and  the  assigned  crop  coefficient  for  each  weather  station  are  listed  in 


Table  3.10.  Boundary  of  the  SFWMD  study  i 


of  the  ground 
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Table  3.10  La 


efficient  for  each  weather  station. 


UTM Crop 

Station  Name  Easting  Northing  CoefF.1  Land  Use/Cover  Type1 


BELLE  GL 
CFSW 
CLEW 
FT.  PI  2 
FT.  PIER 
HGS1 


KISS.FS 
LALFEX 
L 002 
LOOS 
L 006 
LEHIGH 
LISBON 
LZ40 

MARTIN  WP 

MBTS 

S140SPW 

SSA 

S61W 

S65C 

S78W 

STAT 

TAMITR40 

VER04W 

WPB.EEDD 


536,821  2,948498  0.70 

510,392  2,956,889  0.70 

508,253  2,958,582  0.70 

564,240  3,034,893  0.70 

548,022  3,026,976  0.70 

491,306  2,968,484  0.65 

546,305  2,789,694  0.90 

456,118  3,129,269  0.60 

429,595  3,108,315  0.80 

521,067  2,995,733  1.00 

502,262  2,981,772  1.00 

521.613  2,966,594  1.00 

435,285  2,942,863  0.70 

423.606  3,193,294  0.65 

520,933  2,975.365  1.00 

573,045  3,013,430  0.60 

558,244  2,791,400  0.95 

517,159  2,894,513  0.90 

563,010  2,951,360  0.60 

465.426  3,106,296  0,70 

488,630  3,030,665  0.95 

469,871  2,962,979  0.70 

571,716  2,942,413  0.90 

517,549  2,848,955  0.95 

553,275  3,056,163  0.70 

593,252  2,955,013  0.60 


Grass  and  Ariculture 
Grass  and  Agriculture 
Gass  and  Agriculture 

Grass  and  Agriculture 
Shrub  and  Brushland 
Mangrove  Swamp 

Hardwood  Swamp 

Open  Water 
Grass  and  Agriculture 

Open  Water 
Barren 

Coastal  Salt  Marsh 
Shrub  Swamp 

Grass  and  Agriculture 
Freshwater  Marsh  & Wet  Prairie 
Grass  and  Agriculture 
Shrub  Swamp 

Freshwater  Marsh  & Wet  Prairie 
Pincland 


1 land  use/cover  type  from  Florida  Land  Cover  (Florida  Game  and  Fresh  Water  Fish 
Commission,  1992) 
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Figure  3.1 1 Boundary  of  the  Soulh  Florida  Water  f 
location  of  the  ground  meteorological  t 
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west  of  Veto  Beach  and  15  miles  east  of  Yeehaw  junction  on  highway  60  (Figure  3.12). 
This  is  a fresh  waler  marsh  area  where  cattail  ( Typha  ilomingensis ) and  sawgrass 
( Caladium  jamaicense)  are  the  dominant  vegetation.  The  east  and  north  side  of  the  marsh 
is  an  agricultural  area  surround  by  levees  and  drained  to  grow  citrus. 

ET  estimated  by  the  remote  sensing  method  was  compared  with  the  lysimeter  ET 
at  two  marsh  environments.  Three  lysimetcis  were  installed  in  the  marsh  area.  To 

large  homogeneous  marsh  area.  The  lysimeters  were  installed  within  marsh  areas  with 
the  vegetation,  meteorological  condition  and  fetch  being  the  same  as  surrounding 
vegetation.  Lysimeter  1 (80°  4119"  W;  27°34'51"  N),  located  in  a cattail  marsh,  was 
planted  cattail  extracted  from  the  surrounding  marsh.  Lysimeter  2 (80°41'18"  W; 
27°35'12"  N),  located  in  a sawgrass  marsh,  was  planted  sawgrass  extracted  from  the 
surrounding  marsh.  Lysimeter  3 <80°41'14"  W;  27-3523"  N)  was  an  open  water  tank 

water  level  gauge  and  solar  panel  installed  in  the  Fort  Drum  marsh.  The  lysimeters  were 

minimize  oasis  effect  The  main  component  of  the  lysimeter  is  a 3.54  m diameter  and  0.9 
m deep  polyethylene  tank.  Water  levels  in  the  lysimeters  were  monitored  with  analog 

8.95  and  10.55  cm  by  automated  inflow  and  outflow  pump  and  flow  meters.  Data  was 
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stored  in  CR10  data  loggers  (Campbell  Scientific.  Inc.)  and  could  be  manually 
downloaded  or  transmitted  through  radio  transmission. 

ET  of  the  lysimeters  was  computed  by  the  water  balance  equation  expressed  as: 
ET,  = R + I-0  + A (3.20) 

where  ET,  - evapotranspiralion  (or  given  period  t (mm); 

R=  rainfall  (mm); 

I = inflow  by  water  pump  (mm); 

O = outflow  by  water  pump  (mm); 

A = D,  - D,,  change  in  water  level  in  lysimcter  for  period  t (mm). 

A weather  station  was  installed  between  lysimeler  2 and  3 to  monitor  the 
meteorological  parameters  in  die  marsh.  Solar  radiation,  net  radiation,  photosynthetic 
photon  flux  density,  air  temperature,  relative  humidity,  atmospheric  pressure,  wind  speed 


ri 

km 

■ 

CHAPTER 4 

RESULTS  AND  DISCUSSIONS 

4.1  SpwKai  Responses 

4.1.1  Sawgrass  and  Cattail  SiKcTOLCgamofOP 

Radiance  from  sawgrass  ( Caladium  jamaicense),  cattail  ( Typha  domingensis ) and 
the  Spectrlon  standard  white  panel  were  measured  by  GER 1500  spectroradiometer.at  the 
Fort  Dram  ET  experimental  sites  on  Oct.  19, 1996,  Dec.  23,  1996,  Mar.  28, 1997  and 
May  30, 1 997.  Spectral  reflectances  were  calculated  with  the  target  radiance  divided  by 
the  standard  white  radiance.  Reflectances  of  cattail  and  sawgrass  against  the  spectral 
wavelength  for  the  four  dates,  ranging  from  350  tint  to  1 100  nm,  arc  plotted  in  Figure  4.1 
to  Figure  4.4,  respectively.  Each  reflectance  curve  averaged  three  to  five  measurements 
from  different  angles  and  was  independently  calibrated.  As  shown  on  these  four  curves, 
cattail  averaged  higher  reflectance  in  the  red  to  infrared  spectrum  in  three  of  the  four 
measuring  dates  (Oct,  and  Mar.,  1996  and  May,  1997).  The  percent  of  reflectance  was 
directly  related  to  the  vegetation  growing  condition  as,  the  reflectance  percentage 
increased  from  new  transplant  in  Oct.  1996  to  mature  in  late  May  1997. 

4.1. 1.1  Normalized  spectral  comparison 

In  order  to  better  compare  the  difference  in  the  shape  of  sawgrass  and  cattail 
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_ Cattail  Sawgras 


— Cattail  Sawgrass 


350  nm  lo  1 100  nm  in  Dec.  23. 1996*  ^ * 
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Cattail  Sawgrass 


Figure  4.3  Reflectance  curves  of  cattail  and  sawgrass  in  the  spectral  wavelength  from 
350  nm  to  1 100  nm  in  Mar.  28, 1997. 
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reflectance  curves,  each  reflectance  was  normalized  with  respect  to  the  mean  value  of  the 

growing  condition  and,  as  a result,  the  spectral  variation  between  vegetation  species  could 
be  more  clearly  examined.  Figure  4.5  and  4.6  illustrate  the  normalized  sawgrnss  and 

unique  representation  of  leaf  structure,  the  four  normalized  curves  of  each  vegetation 
displayed  highly  similar  shapes.  The  four  dates  of  normalized  sawgtass  and  cattail  curves 
arc  plotted  in  Figure  4.7.  Cattail  and  sawgrnss,  both  being  wetland  sedge,  with  the 
exception  of  a few  spectral  regions,  displayed  similar  normalized  curves. 

In  order  to  locate  the  bonds  for  more  efliciently  discriminating  between  sawgrnss 
and  cattail,  the  spectrum  was  divided  into  five  narrow  bands:  blue  (425-475  nm);  green 
(525-575  nm);  red  (625-675  nm);  near  infrared  (725-800  nm);  and  first  infrared  (NIR) 
recess  (925-975  nm).  The  mean  spectral  difference  at  these  five  bands  was  computed  for 
the  original  reflectances  and  the  normalized  reflectances  in  Tabic  4.1  and  4.2, 
respectively.  As  evident  in  the  two  tables,  the  reflectance  difference  of  the  normalized 
curves  reduced  significantly.  The  most  consistent  deviation  occurred  at  the  red,  near 
infrared  and  the  first  NIR  recess  bands  indicating  that  these  three  spectral  regions,  625- 
675  nm,  725-800  nm,  and  925-975  nm  should  be  given  higher  priority  in  band  selection 
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Wavelength  (nm) 


Figure  4.5  Normalized  reflectance  curves  of  sawgrass  ( Caladium  jamaicense ) at  four 
dates.  Note  that  the  curves  arc  almost  identical  curves  except  in  few 
spectral  regions. 


Wavelength  (nm) 


Figure  4.6  Normalized  reflectance  curves  of  cattail  ( Typha  domingensis ) at  four 


Table  4.1 


Spectral  difference  of  cattail  and  sawgrass  at  6 wavelength  intervals  on 
four  dates.  Test  statistics  arc  in  parenthesis.  Differences  are  calculated  as 
sawgrass  reflectance  subtracted  by  cattail  reflectance. 


Wavelength 

(nm) 


10-19-96  12-23-96  03-28-97  05-30-97 


425-475 
525-575 
625-675 
725-800 
925-975 
350-1 100 


blue  -0.4362 

(0.0224) 

green  -1.9854 

(0.1018) 

red  -7.7271 

(0.3962) 

NIR  -6.0236 

(0.3088) 

first  NIR  recess  -7.4587 

(0.3824) 


2.2107  -0.2599  -2.6178 

(0.1050)  (0.0103)  (0.0866) 

5.7769  -5.4344  0.0148 

(0.2744)  (0,0083)  (0.1797) 

2.5379  -1.8091  -4.8227 

(0.1205)  (0.0716)  (0.1595) 

12.6335  -3.1510  -17.5021 

(0.6001)  (0.1247)  (0.5787) 

11.1177  -1.2871  -9.6142 

(0.5281)  (0.0509)  (0.3179) 


-5.2894  7.3169  -1.7452  -8.9560 

(0.2712)  (0.3475)  (0.0690)  (0.2961) 


ormalizcd  ? 


of  cattail  and  sawgrass  at  6 wavelength 
tistics  are  in  parenthesis.  Differences  are 
iss  reflectance  subtracted  by  normalized 


Date 

Wavelength  Interval  

(nm)  location 10-19-96  12-23-96  03-28-97  05-30-97 


425-475 


0.0195  0.0261  0.0028  -0.0122 

(0.0301)  (0.0394)  (0.0045)  (0.01%) 


525-575 


0.0148  0.0873  0.0144  -0.0251 

(0.0229)  (0.1318)  (0.0227)  (0.0404) 


-0.1811  -0.0209  -0.0253  -0.0115 

(0.2799)  (0.0317)  (0.0400)  (0.0186) 


725-800 


0,0631  0.0661  -0.0152  -0.0923 

(0.0977)  (0.0998)  (0.0241)  (0.1484) 


0.0501  -0.0427  0.0398  0.1020 

(0.0775)  (0.0646)  (0.0628)  (0.1641) 


350-1100 


0.0000  0.0000  0.0000  0.0000 

(0.0000)  (0.0000)  (0.0000)  (0.0000) 


' sawgrass  rcflcctanc 
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Vegetation  with  similar  spectral  response  are  difficult  to  separate  in  satellite 
image  classification.  With  spectral  measurements  from  field  spcctroradiometer  made  in 

vegetation  types.  In  Figure  4.7,  the  normalized  reflectance  curves  are  shown  for  different 
limes  of  the  year.  These  curves  indicates  that  winter  (Oct  19, 1996  and  Dec.  23. 1996)  is 


AVHRR 


The  relationship  between  spcctroradiometer  and  NOAA  AVHRR  were  evaluated 
through  the  red,  near  infrared  bands  and  their  derived  vegetation  indices.  The  results  are 


4.1. 2.1  Comparison  of  spectral  bands  and  vegetation  indices 

Spectral  response  of  red  and  near  infrared  bands,  as  well  as  generated  vegetation 
indices  ftom  spcctroradiometer  and  from  NOAA  AVHRR  on  the  four  dates  are  listed  in 
Table  4.3.  Due  to  the  available  data  archive  in  the  NOAA  Satellite  Active  Archive, 
AVHRR  data  from  the  nearest  available  dates  of  Oct  21, 1996  and  Dec.  24,  1996  were 
used  to  compare  against  the  field  measurements  of  Oct.  19, 1996  and  Dec.  23, 1996, 
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spectroradiometcr  measurements  were  significantly  higher  than  that  of  the  satellite 
measurements,  both  regarding  individual  band  reflectance  and  vegetation  indices.  This  is 

was  measured  in  a short  distance  with  the  subject  vegetation  occupying  most  of  the  view, 
while  satellite  measurements  were  taken  from  high  above  and  unavoidably  aggregated 
with  some  non-vegetated  background  (open  water,  exposed  soil  surface)  in  the 


4. 1.2. 2 Comparison  ef.correlwiQn 

The  correlation  of  the  spectral  bands  i 


id  vegetation  indices  between  field  and 
using  regression  analysis.  The  slope, 


between  field  and  satellite  measurements  arc  listed  in  Table  4.4.  As  shown,  vegetation 
indices,  in  general,  display  a higher  correlation  than  individual  spectral  bands  except  for 

correlation  to  vegetation  characteristics  than  do  the  individual  spectral  bands.  Among  the 
three  vegetation  indices,  NDVI  showed  the  highest  correlation  of 0.904. 


4.1.3  Florida  Vegetation  Reflectance  Curves 

Selected  reflectance  curves  of  Florida  wetland  vegetation  (cattail,  sawgrass, 

vegetation  (pine,  oak,  palmatto)  have  been  collected  by  the  GER  1 500  spectroradiometcr. 
Spectral  reflectance  was  calibrated  with  a Spcctralon  white  standard  plate.  The  curves  are 
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illustrated  in  Appendix  1 . Each  vegetation  response  curve  is  an  average  of  three  to  five 
measurements  from  different  angles.  The  upper  and  lower  90%  confidence  intervals  are 


Table  4.4  Correlation  analysis  between  spectroradiometer  and  NOAA  AVHRR 
spectral  bands  and  derived  vegetation  indices.  Data  was  calculated  from 


RVI 

NDV1 


0.58  - 0.68  pm 
0.73 -1.05  pm 

JR 

R 

IR-R 

IR+R 


0.78  7.80  0.496 
1.84  23.53  0.813 
0.61  1.71  0.791 
0.86  0.21  0.904 


SAVI 


IR-R 
IR  + R + L 


:<■  + £) 


0.87 


0.32  0.898 


The  collection  of  spectral  response  curves  of  Florida  vegetation  can  be  used  as  a 
spectral  signature  database.  This  database  can  aid  in  determining  the  appropriate  number 
and  range  of  spectral  bands  for  vegetation  discrimination. 


The  Landsat  TM  images  were 


3000  m,  respectively.  Figure  4.8  illus 
(c).  and  900  m (d)  spatial  resolutions. 


user  resolutions  by  averaging  their 
,5.6.8, 10, 20, 30, 40, 50.  turd  100 
>40, 300, 600, 900, 1200, 1500  and 
area  in  120  m (a),  300  m (b),  600  m 


Using  the  original  30  m pixel  size,  the  mean  and  standard  deviation  for  each  of 
the  land  use/cover  images  was  calculated  with  ERDAS  IMAGINE  (ERDAS,  1996)  image 
processing  software  and  the  results  are  listed  in  Table  4.5.  As  expected,  as  a result  of  the 
high  land  cover  variability  in  the  urban  area,  the  urban  area  had  the  highest  standard 

produced  the  lowest  standard  deviation. 


The  Mean  Local  Coefficient  of  Variance  for  the  near  infrared  band  of  the  Landsat 
TM  data  was  calculated  for  each  of  the  four  land  use/covcr  types  (urban,  agriculture, 
wetland,  and  mixture)  at  14  different  spatial  resolutions.  The  MLCV  and  its 
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Figure  4.8 


Landsat  TM  images  in  four  levels  of  spatial-resolution  degradation. 
Images  show  a residential  area  west  of  Fort  Lauderdale  in  1 20  m (a), 
300  m (b),  600  m (c),  and  900  m (d)  resolutions. 
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Table  4.5  Size,  local 

ion.  mean  and  standard  deviation  (S.D.)  for  the  selected 

Land  Cover  No.  of 

Type  Pixels 

Near  Infrared  Thermal 

(Km1)  Mean"  S.D."  Mean"  S.D." 

Urban  360,000 

324  70.64  23.82  136.07  4.28 

Agriculture  1,200,000 

1,080  44.99  21.22  136.17  4.50 

Wetland  1,200,000 

1,080  46.84  9.62  130.43  1.43 

Mixture  12,000,000 

10,800  52.46  20.34  132.34  4.67 

‘ Area  was  calculated  with 

original  pixel  size  30  m * 30  m. 

Mean  and  Standard  Devia 

ition  (S.D.)  values  arc  for  the  Digital  Number  (DN)  of  the 

: for  the  Digital  Number  (DN)  of  the 
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corresponding  resolution  are  listed  in  Table  4.6.  The  graphs  of  MLCV,  as  a function  of 
spatial  resolution,  are  plotted  in  Figure  4.9.  As  the  peaks  of  the  graphics  indicate,  the 
peak  of  urban,  agricultural,  wetland  and  mixture  lines  are  at  60  m,  300  m,  1200  m,  and 
300  m respectively. 

4.2.3  MLCV  of  Thermal  Images 

The  Mean  Local  Coefficient  of  Variance  for  the  thermal  band  of  the  Landsat  TM 
data  was  calculated  for  each  of  the  four  land  use/cover  types  at  14  different  spatial 
resolutions.  The  thermal  band  MLCV  and  its  corresponding  resolution  are  listed  in  Table 
4.7.  The  graphs  of  MLCV,  as  a function  of  spatial  resolution,  are  plotted  in  Figure  4.10. 
As  the  peaks  of  the  graphics  indicate,  the  peak  of  urban,  agricultural,  wetland  and  mixture 
lines  are  at  300  m,  450  m,  1 800  m,  and  600  m respectively. 

These  peaks  in  the  graphs  of  MLCV  against  spatial  resolution  identify  the  optimal 

land  use/cover  conditions.  For  example,  the  agricultural  fields  in  the  Everglade 
Agricultural  Area  measured  an  appropriate  resolution  of  300  m which  is  approximately 

than  the  size  of  the  sugarcane  fields,  the  neighboring  pixels  are  within  a single  stalk  or 
clump  of  sugarcane,  which  are  spectrally  homogeneous  and  results  a low  MLCV.  A 


-MLCV  i 


Tabic  4.6 


6 Mean  Local  Coefficient  of  Variance  (MLCV)  for  near  infrared  band  for 
urban,  agriculture,  wetland,  and  mixture  land  use/covcr.  Bold  number 
indicates  peak  value  location. 


Spatial 

30 


60 

90 

120 

130 


1300 

3000 


Urban 

2.0398 

2.4318 

2.3522 

2.1902 

2.1419 

2.0170 

1.7832 

1.6075 

1.2830 

1.0045 

0.7621 

0.5318 


Land  Use/Cover  Type 

Aerie. Wetland Mixture 

0.9549  0.3004  0.6428 

1.5968  0.3495  0.9175 

1.9933  0.3349  1.0069 

2.2697  03189  1.0568 

2.4903  03132  1.0898 

2.6262  0.3107  1.1111 

2.7820  03176  1.1321 

2.8419  0 3285  1.1357 

2.6995  0 3497  1.0955 

2.6079  0 3719  1.0645 

2.1887  0.3862  0.9949 

1.7675  0.4067  0.9650 

1.4726  0.1585  0.9754 

0.1174  0.9789 
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Spatial  Resolution  (m) 

-o-  urban  agric.  -o-  wetland  mixture 


Figure  4.9  Near  Infrared  band  Mean  Local  Coefficient  of  Variance  (MLC  V)  for 

urban,  agriculture,  wetland,  and  mixture  land  use/cover. 


Mean  Local  Coefficient  of  Variance  (MLCV)  for  thermal  band  for  urban, 
agriculture,  wetland,  and  mixture  land  usc/cover.  Bold  number  indicates 
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— o—  urban  agric.  -o-  wetland  mixture 


Figure  4.10  Thermal  band  Mean  Local  Coefficient  of  Variance  (MLCV)  for  urban, 
agriculture,  wetland,  and  mixture  land  usc/covcr. 
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agricultural  fields  aggregated  within  a single  pixel.  The  appropriate  resolution,  therefore, 
is  indicated  by  a peak  in  the  MLCV  graph,  which  occurs  when  the  resolution  is  close  to 

high. 

4.2.5  Comparison  of  Near  Infrared  and  Thermal  Ontiroal  Resolution 

Graphs  of  near  infrared  and  thermal  band  MLCV  as  a function  of  resolution  arc 
compared  for  the  four  land  cover  types  (urban,  agriculture,  wetland,  and  mixture)  in 

and  1800  m for  near  infrared  and  thermal  band,  respectively.  This  is  approximately  the 
same  resolution  of  a weather  satellite  such  as  the  NOAA/AVHRR  data  which  has  1100  m 
nadir  resolution.  In  other  words,  the  spatial  scale  of  the  NOAA/AVHRR  imagery  is 
suitable  for  wetland  studies  in  south  Florida.  The  graph  of  the  mixture  land  usc/covcr 
had  a peak  at  300  m and  600  m,  indicating  that  300  m and  600  m is  the  preferred  spatial 
resolutions  for  near  infrared  and  thermal  bands  if  only  one  resolution  can  be  selected  for 
the  entire  south  Florida. 
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Spatial  Resolution  (m) 


! (MLCV)  of  the  near  infrared  and 
ver.  Peaks  of  the  near  infrared  and 
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Figure  4.12 


Mean  Local  Coefficient  of  Variance  (MLCV)  of  the  near  infrared  and 
thermal  bands  for  agriculture  land  use/cover.  Peaks  of  the  near  infrared 
and  thermal  bonds  are  located  at  300  and  4S0  m,  respectively. 


Spatial  Resolution  (m) 


Figure  4.13 


Mean  Local  Coefficient  of  Variance  (MLCV)  of  the  near  infrared  and 
thermal  bands  for  wetland  land  usc/covcr.  Peaks  of  the  near  infrared  and 
thermal  bands  are  located  at  1 200  and  1 800  m.  respectively. 
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Mean  Local  Coefficient  of  Variance  (MLCV)  of  the  near  infrared  and 
thermal  bands  arc  located  at  300  and  600  m,  respectively. 
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4.3.1  Field  Temperature  apd  NDVI  MManfflatt 

Field  surface  lemperalures  were  collected  on  several  days  coincident  with  NOAA 
satellite  overpass  times  (0230h  and  1430h).  Surface  radiative  temperatures  of  selected 

(Model  210).  Surface  radiative  temperatures  were  taken  in  the  four  cardinal  directions 
(North,  East,  South,  and  West)  and  averages  were  computed.  Field  land  cover  spectrum 
were  measured  with  a GER-1500  field  spectro-radiometer  and  downloaded  to  computer 
via  RS232  port  for  further  analysis.  Surface  temperature  and  NDVI  derived  from  field 
measured  spectrum  reflectance  curves  for  several  selected  ground  covers  (asphalt, 

water  with  30%  vegetation  cover  and  water  with  0%  vegetation  cover)  are  listed  in  Table 
4.8.  Spectrum  reflectance  curves  for  vegetation,  soil  and  water  samples  for  the  selected 
ground  objects  are  plotted  in  Figure  4. 15. 

4.3.2  Image  Processing  of  NOAA  AVHRR 

NOAA  AVHRR  data  in  10-bit  packed  format  were  acquired  from  the  NOAA 
Satellite  Active  Achieve  via  the  internet  file  transfer  and  processed  by  my  AVHRR 
processing  program  and  ERDAS  Imagine  software  package.  Raw  images  and  ancillary 
data  were  extracted.  Figure  4. 16  illustrates  the  raw  Florida  image  from  the  NOAA 
satellite  with  the  Florida  state  boundary  in  the  background.  It  is  evident  that  the  image  is 
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Figure  4.15  Spectrum  reflectance  cu: 
bare  soil,  grass,  cabbage 


rves  for  selected  ground  objects  (asphalt,  concrete, 
iter  with  0%  vegetation  cover)  as  measured  by 
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distorted  due  to  scan  skew  and  earth  curvature.  Geometric  rectification  was  applied  to 
correct  the  geometric  deformations  of  the  satellite  images.  Figure  4.17  shows  the  gco- 
rcctificd  Florida  image  on  December  9,  1996. 

Normalized  difference  vegetation  index  (NDVI)  was  derived  from  the  infrared  and 
red  bands  of  NOAA  AVHRR  data.  NDVI  in  Florida  ranged  from  -0.05  to  0.6S  mainly 
depends  on  the  plant  density  and  vegetation  vigor.  Figure  4. 1 8 illustrates  the  NDVI  map 
of  Florida.  In  the  color  scheme  of  NDVI  map,  the  greener  the  color,  the  higher  the 
NDVI.  The  upper  St.  John's  river  basin.  Everglade  Agriculture  Area,  and  the  southwest 
coast  marsh  display  the  highest  NDVI.  Open  water  bodies,  urban  areas,  and  wetland  areas 

Surface  temperature  was  derived  and  calibrated  with  the  corresponding  hourly  lake 
surface  temperature  Diurnal  surface  temperature  was  computed  as  the  difference  of  the 
surface  temperature  at  day-time  pass  (~1400h)  and  at  night-time  pass  (~0200h).  Figure 
4.19  illustrates  the  satellite-derived  surface  temperature  on  Dec.  9,  1996.  As  shown,  there 
is  an  apparent  trend  of  temperature  decreasing  with  latitude;  the  further  north  the  location, 
the  lower  the  surface  temperature.  Urban,  metropolitan  areas  such  as  Miami,  and  West 
Palm  Beach  showed  high  surface  temperature  while  lakes  and  wetlands  displayed  lower 
overall  temperatures  in  the  daytime  satellite  scene. 


NOAA  AVHRR  image  of  December  9, 1 996  with 
ackground.  The  image  is  geometrically  distorted  di 


4.17  Geometrically  rectified  NOAA  AVHRR  image  of  December  9, 1996. 


Figure  4. 1 8 Normalized  difference  vegetation  index  derived  from  NOAA- 1 4 
AVHRR  image  for  December  9,  1996. 
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Figure  4. 1 9 Land  surface  temperature  derived  from  NOAA- 1 4 AVHRR  images  for 
December  9. 1996. 
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4.3.3  Potential  ET  by  Surface  Temperature  Method 

The  relationships  between  R„-ETani  SAW  were  calculated  using  a linear 
regression  model  and  data  from  the  34  weather  stations  on  four  dates.  The  regression 


Table  4.9  Coefficients  A and  B of  the  energy  balance  equation,  and  the  coefficient  of 

determination  (rJ)  from  the  regression  analysis  of  7V  vs  (R.-E7) 

Date A B r 

14  Dec.  1989  -0.40  0.26  0.64 

20  May  1990  -0.71  0.79  0.80 

12  Dec.  1992  -0.37  0.44  0.69 

11  Apr.  1993  -0.48  0.31  0.74 

Mean  -0.49  0.45 

Std  Dev.  ±0.133  ±0.207 

Coefficient  A,  the  intercept  of  the  /f„-£7’and  the  T,-Ta  regression  line,  ranged  from 
-0.40  to  -0.71  with  a mean  of -0.49  and  standard  deviation  ofO.133.  Coefficient  B,  the 
slope  of  the  regression  line,  ranged  from  0.26  to  0.79  with  a mean  of  0.45  and  standard 
deviation  of0.207.  The  results  are  in  good  agreement  with  the  results  concluded  by  other 
researchers.  For  instance,  the  average  B value  of  0.45  fell  within  the  range  of0.25  as 
reported  by  Scguin  and  Itier  (1983)  and  0.64  as  reported  by  Jackson  ct.  al.  (1977). 


Average  of  the  r1  values  is  0.72.  Figure  4.20  illustrates  the  relationship  between  SAW  and 


(R.-E1)  on  1 1 Apr.  1993,  i 


As  Table  4.9  shows,  the  wel  season  images  (20  May  1990  and  1 1 Apr.  1993) 
appear  to  have  higher  r1  values  (0.80  and  0.74)  as  compared  with  the  other  two  dty  seasor 
images  (14  Dec.  1989  and  12  Dec.  1992)  which  had  r2  values  of  0.64  and  0.69.  This 
might  be  due  to  the  surface  energy  flux  that  is  sensitive  to  the  dense  cover  of  vegetation 
grown  in  the  wet  season.  In  other  words,  using  satellite  derived  surface  temperature  to 
estimate  the  regional  ET  in  the  wet  season  may  be  slightly  better  than  that  used  in  the  dry 
season.  However,  further  studies  are  necessary  to  support  this  observation. 

ET  simulated  with  the  established  relationship  from  the  satellite-based  data  is 
plotted  in  Figure  4.21  against  ET  calculated  from  the  ground-based  data  for  the  34 
stations.  The  relation  of  the  satellite-based  ET,  AT,,,  and  the  ground-based  ET,  ET^  is 
expressed  as  follows: 

ETm  = 0,74  * 0.79  *£7^  r*± 0.76 

The  r1  value  of  0.76  implies  that  the  satellite-based  ET  estimate  is  an  indication  of 
the  ground-based  ET  estimation.  However,  the  slope  of  regression  was  0.79,  not  1.0. 

This  means  that  there  was  slight  discrepancy  between  two  approaches  in  ET  estimation. 

In  the  meantime,  it  should  be  noted  that  the  ground-based  ET  is  a point  estimation,  while 
the  satellite-derived  ET  is  a manifestation  of  an  I km!  area.  This  spatial  coverage 
difference  could  cause  a deviation  between  two  approaches. 


Ts-TaCC) 


Fig.  4.20  Relationship  of  the  difference  between  net  radiation  and 

evapotranspiration  ( R„-ET)  and  the  surface-to-air  temperature 
(TrTJ  on  1 1 Apr.  1993. 


Simulated  ET  from  the  linear  relationship  between  ASTV and  Rn-ET. , 


Fig.  4.21 
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4.3.4  NPV1  ant!  D.STVXriangje 

Relationship  of  surface  temperature  and  vegetation  index  was  evaluated  with  both 
the  field  measured  results  and  satellite  data.  Day-time  (14:00  to  15:00)  surface 
temperature  Td,  air  temperature  Ta,  night-time  (02:00  to  03:00)  surface  temperature  Tn. 

and  spectroradiometer.  Figures  4.22  to  4.24  show  the  NDVI  and  surface  temperature 
(TJ,  diurnal  temperature  difference  (Tj  - TJ,  and  surface-air  temperature  difference  (T,  - 
TJ  relationship  for  the  selected  ground  objects,  respectively.  It  can  be  seen  that  the 
points  with  the  highest  diumal  surface  temperature  difference  are  dry  surfaces  (asphalt, 
concrete,  dry  bare  soil);  the  points  with  the  lowest  diumal  surface  temperature  difference 
are  wet  surfaces  (water  covered  by  vegetation,  wet  bare  soil);  and  the  points  with  the 
highest  vegetation  index  arc  grass  and  vegetables.  These  points  generally  form  a triangle 
with  dry  surface,  wet  surface,  and  vegetation  at  three  comers. 


Satellite 


I NDVI 
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Figure  4.22  NDVI  and  diurnal  temperature  difference  (Td  - Tn)  relationship  for  the 
selected  ground  objects  measured  in  the  field. 


Figure  4.23  NDVI  and  surface-air  lempcralure  difference  (Ts  - Ta)  relationship  for  the 
selected  ground  objects  measured  in  the  field. 
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Ts(-C) 


Figure  4-24  NDVI  and  surface  temperature  relationship  for  the  selected  ground  objects 
measured  in  the  field. 
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Figure  4.2S  Satellite  derived  normalized  difference  vegetation  index  (NDVI)  versus 
diurnal  surface  temperature  deviation  (DSTV)  for  seven  land  use/cover 
types.  Agriculture,  urban,  and  marsh  are  die  three  comers  of  the  NDVI- 
DSTV  triangle. 


S I 1 ! t i I 


TheNDVI-DSTV  triangle  relationship  was  used  to  calculate  weighting  coefficients 
ofVMC  as  defined  in  Equation  3.15  to  3.17.  The  DSTV,  NDV1  and  weighting  foctois  for 
the  eight  selected  land  use/cover  on  Dec.  9,  1996  are  listed  in  Table  4.10.  Land  cover 
types,  Grass  and  Agriculture,  Barren  and  Urban  Land,  and  Fresh  Water  Marsh,  are 
assigned  as  the  three  comers  (agriculture,  dty  land,  and  wet  land)  of  NDVI-DSTV 
triangle,  respectively.  The  DSTV  and  NDVI  were  average  values  of  all  pixels  within  the 
selected  sites.  Weighting  factors.  w„  w„  Wj  indicate  the  resemblance  of  the  selected  land 
cover  type  to  agriculture,  wet  land,  and  dry  land,  respectively. 

As  shown  in  Table  4. 10,  the  mangrove  swamp  and  cypress  swamp  have  higher 
values  of  w,  weights  from  wet  land  (fresh  water  marsh).  The  VMC  can  be  determined 
based  on  the  NDVI-DSTV  triangle  relationship  for  locations  with  various  fraction  of 
vegetation  cover  and  soil  moisture  due  to  the  fact  that  the  points  at  the  comers  of  the 
NDVI-DSTV  triangle  arc  relatively  easy  to  measure  or  can  be  reasonably  assumed. 

The  potential  ET  values  at  ground  weather  station  were  computed  by  the  Penman 
method  and  interpolated  to  the  entire  study  area  using  the  Krigning  interpolation  algorithm 
with  linear  variogram.  The  interpolation  was  implemented  using  the  software  package 
SURFER  for  Windows  (Golden  Software  Co.,  Golden,  CO). 
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Table  4. 10  DSTV,  NDVI  and  weighting  factors  for  the  eight  selected  land  use/cover 

on  Dec.  9, 1996. 


Land  Cover  Type'  DSTV  NDVI 


Grassland  and  Agriculture 
Barren  and  Urban  Land 
Fresh  Water  Marsh 

Mangrove  Swamp 


Cypress  Swamp 


1.0000  0.0000  0.0000 

0.0000  0.0000  1. 0000 

0.0000  1.0000  0.0000 

0.5686  0.3063  0.1251 

0.3502  0.6643  0.0056 

0.4852  0.4534  0.0614 

0.7668  -0.0098  0.2429 

0.2120  0.5679  0.2201 


1 Land  cover  types  used  in  the  Florida  Land  Cover  Map  (Florida  Game  and  Fresh  Water 
Fish  Commission,  1992) 

A geographic  information  system  was  constructed  with  NOAA  AVHRR  raw 
bands,  satellite  derived  NDVI,  surface  temperature,  DSTV  and  Florida  land  use/cover 
map  Each  of  the  data  type  was  built  as  a data  layer  covering  the  study  area  in  the  GIS. 
The  three  weighting  factors,  VMC,  and  the  interpolated  potential  ET,  were  also  calculated 
and  integrated  into  the  GIS.  The  ET  rates  were  then  determined  as  the  potential  ET  layer 
multiplying  by  the  VMC  layer  for  each  of  the  grid  cell  in  the  GIS.  Figure  4.26  illustrates 
the  satellite-derived  and  ground-based  data  layers  and  the  calculation  of  ET  in  the 
geographic  information  system. 


Satellite-derived  and  ground-based  data  layers  and  the  calculation  of  ET  in 
the  geographic  information  system. 


Figure  4.26 


4.3.6  Verification 

Verification  of  ET  from  the  remote  sensing  method  was  made  with  traditional  ET 
methods  and  with  lysimeter  ET  methods. 

4.3.6. 1 Verification  with  traditional  ET  methods 

The  ET  estimated  by  the  remote  sensing  method  versus  ET  estimated  by 
traditional  (Penman)  method  at  1 1 weather  stations  on  the  five  days  is  shown  in  Figure 
4.27.  Their  relationship  can  be  shown  as 

ET**  = 1 .022  ET 0.712,  r1  = 0.647,  n - 55 

where  ETt!  is  the  remote  sensing  method  ET,  and  ET, is  the  Penman  method  ET. 

remote  sensing  ET  and  the  tradition  method  ET.  However,  afler  a close  inspection  of  the 
weighting  factors  and  the  VMC  of  the  stations,  an  excessive  low  VMC  value  was  found  at 
station  JBTS  (Joy  Bay  Weather  Station,  25°  13'  N,  80"  32'  W).  Station  JBTS  was  located 
at  the  southern  tip  of  the  Everglades  National  Park  adjacent  to  the  Florida  Bay.  Due  to 
the  low  infrared  reflectance  of  water,  JBTS  displayed  NDVI  values  below  zero  that 
produced  VMC  values  for  JBTS  much  lower  than  the  crop  coefficient  used  for  mangrove 
swamp  in  traditional  ET  calculation.  The  relationship  of  remote  sensing  ET  and  the 
traditional  method  ET  was  re-evaluated  excluding  the  five  JBTS  points,  and  the 


r3  = 0.785,  n = 50. 


3.588, 


158 

The  correlation  coefficient  was  significant  to  indicate  a positive  relationship 
between  the  ET  estimated  by  the  remote  sensing  method  and  by  the  traditional  method. 

The  NDVi-DSTV  method  can  not  be  used  with  open  water  surface  due  to  the 
negative  NDV1  which  results  in  a negative  VMC  in  the  calculation.  However,  this 

technique.  VMC  values  can  be  assigned  to  open  water  surfaces. 

Stations  located  at  the  boundary  of  a large  water  body  can  create  highly  unreliable 
results  due  to  large  pixel  scale  and  spectral  signature  mixing.  The  large  ground  coverage 

high  percentage  of  water  present  within  the  area  of  satellite  field  of  view. 

4.3.6  2 Verification  With  lysimeter  ET 

Remote  sensing  estimated  ET  were  compared  with  the  lysimetcr  water  balance  ET. 
recorded  ET  for  sawgrass  and  cattail,  respectively. 


ET  by  Penman  Melhod  (mm/day) 

• 08-17-96  * 09-24-96  *12-09-96  * 03-28-97  • 05-30-97 


27  ET  estimated  by  the  remote  sensing  method  versus  ET  estimated  by  the 
traditional  (Penman)  method  on  the  five  days.  ET  values  at  the  Jay  Bay 
weather  station  (JBTS)  on  the  five  dates  were  identified. 


^et  et0 
(mm/day)  (mm/day) 


DSTV  VMC 


R.  S,  Error 
ET 

(mm/day)  (%) 


08- 18-1996  5.84 

09- 24-1996  5.33 
03-28-1997  3.30 
05-30-1997  5.59 


1 cloud  contaminated. 


Table  4. 12  Lysimeter  recorded  ET  and  remote  sensing  estimated  ET  for  cattail  marsh 
in  four  days. 


Lysimeter  Penman  NDV1  DSTV  VMC  R S.  Error 
ET  ET0  ET 

(mm/day)  (mm/day)  (mm/day)  (%) 


08- 18-1996  N/A‘  5.98 

09- 24-1996  6.09  5.45 
03-28-1997  4.06  5.32 
05-30-1997  4.83  6.58 


3.08  -23.92 


As  shown  in  Table  4, 1 1 and  4. 12,  the  average  error  of  remote  sensing  and 
lysimeter  recorded  ET  was  8.23  and  14.92  percent  for  sawgrass  and  cattail,  respectively. 
The  largest  estimation  error  (-23.92%)  was  with  the  cattail  lysimeter  on  March  28,  1997. 
On  that  date,  the  satellite  image  pixel  at  the  cattail  lysimeter  location  was  not  identified  as 
cloud-contaminated  while  the  sawgrass  pixel  was.  However,  the  NDVI  of  the  cattail  was 
significantly  lower  than  the  average  NDVI  from  other  dates.  The  cause  of  the  low  NDVI 
value  might  be  due  to  the  existence  of  thin  cloud  which  was  not  detected  by  the  cloud 
screening  algorithm.  The  low  NDVI  resulted  in  a low  VMC  and  consequently  a low 
estimated  ET.  If  we  assumed  this  pixel  was  also  contaminated  by  cloud  cover  and 
excluded  from  the  error  calculation;  the  average  error  reduced  to  10.42  and  9. 1 percent 
for  cattail  lysimeter  and  for  the  two  lysimeters  combined,  respectively.. 

The  existance  of  clouds  was  a significant  problem  in  applying  remote  sensing 
methods  for  regional  ET  calculation.  In  the  NDVI-DSTV  algorithm,  both  day-  and  night- 
time satellite  images  must  be  cloud  free  in  order  to  correctly  derive  the  diurnal  surface 
temperature  variation.  However,  most  of  the  NOAA  AVHRR  satellite  images  contain 
some  clouds  over  an  extended  area.  This  situation  is  particularly  severe  in  the  summer  and 
fall  seasons  and,  in  most  night  satellite  images,  low  cirrus  clouds  and  fog  are  present  in  the 
images.  In  fact,  there  is  no  cloud-free  day-and-night  Florida  image  pair  in  the  more  than 
seventy  NOAA  AVHRR  images  studied  in  the  period  from  July  and  September,  1996, 

The  VMC  values  of  the  three  comers  of  the  NDVI-DSTV  triangle  arc  important  in 
determining  VMC  at  each  pixel.  The  wet  and  dry  land  VMC  were  assumed  to  be  0.95  and 
0.05,  respectively  as  the  high  and  low  actual  ET  to  potential  ET  ratios  for  the  entire  year 
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However,  VMC  should  be  determined  at  least  on  a seasonal  basis.  Average  monthly 
ratios  of  actual  to  potential  ET  in  a long  term  lysimcter  study  from  marsh  area  in  south 
Florida  could  provide  a more  reliable  monthly  VMC  for  the  ET  estimation  in  Florida. 

4.3.7  Regional  ET  in  Soulh  Florida 

Regional  ET  in  south  Florida  was  calculated  by  using  the  VMC  and  the  potential 
ET  layers  in  the  GIS.  The  ET  contour  map  is  plotted  in  Figure  4.28.  ET  contour  map  of 
interpolated  ET  on  May  30,  1996  is  also  plotted  in  Figure  4.29. 

As  shown  in  both  Figures  4.28  and  4.29,  the  highest  ET  rates  appears  over  the 
Lake  Okeechobee.  However,  the  interpolated  ET  in  traditional  method  displays  a higher 
value  along  the  SFWMD  boundary  due  to  the  extrapolation  from  further  distance  away. 
ET  contour  map  by  the  remote  sensing  method,  on  the  other  hand,  displays  a more  evenly 
distributed  and  detailed  contour  due  to  the  spatially  distributed  VMC  from  remote  sensing 
method  The  pattern  of  excessive  values  along  the  extrapolating  boundary  is  more 
noticeable  if  there  are  fewer  weather  stations  in  the  region. 


811  May  30,  IW.  ***  y 8 
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Figure  4.29  Contour  map  of  regional  ET  in  south  Florida  by  interpolation  of  traditional 
ET  method  on  May  30.  1997. 


CHAPTERS 

CONCLUSIONS  AND  RECOMMENDATIONS 


5,1  Conclusions 

The  hand-held  spectroradiomcter  is  useful  in  selecting  the  most  efficient  bands 

measured  in  different  months  can  be  used  to  determine  the  optimal  time  span  to  perform 
land  use/covcr  classification  of  satellite  images  for  areas  with  similar  vegetation,  in  this 
study,  the  most  efficient  bands  for  sawgrass  and  cattail  discrimination  were  found  at  red 
(625  to  675  nm).  near  infrared  (725  to  800  nm).  and  first  IR  recess  (925  to  975  nm).  As 
results  show  in  the  seasonal  reflectance  comparison,  fall-early  winter  was  the  best  season 
to  distinguish  cattail  and  sawgrass.  Hand-held  radiometer  measurements  and  satellite 
sensed  data  were  found  to  have  a higher  correlation  through  vegetation  indices  as  opposed 
to  the  raw  spectral  bands.  By  applying  the  ratio  and  difference,  the  vegetation  indices  can 

produce  a higher  correlation  between  hand-held  radiometer  and  satellite  sensor. 
Normalized  Difference  Vegetation  Index  (NDVI)  displayed  the  highest  correlation  among 
the  vegetation  indices  studied. 

Based  on  the  results  of  mean  local  coefficient  of  variance  (MLCV)  curves  in  the 
study,  the  optimum  spatial  resolution  for  classifying  urban,  agriculture,  and  wetlands  in 
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Florida  is  60, 300,  and  1200  m,  respectively  for  the  near  infrared  band.  The  optimum 
spatial  resolution  is  300, 4S0,  and  1800  m for  urban,  agriculture  and  wetland,  respectively, 
for  the  thermal  infrared  band.  For  south  Florida,  if  only  one  resolution  can  be  selected,  the 
optimum  spatial  resolution  is  300  m and  600  m for  near  infrared  and  thermal  band, 
respectively.  As  to  satellite  imagery,  these  are  the  most  appropriate  resolutions  to  identify 
the  land  use  activities  for  respective  regions.  A potential  ET  study  using  NOAA  AVHRR 
thermal  bands  at  a resolution  close  to  the  optimal  spatial  resolution  resulted  in  a positive 
relationship  between  satellite  and  ground  based  ET. 

Vegetation  index  and  surface  temperature  extracted  from  the  NOAA  AVHRR  data 
are  related  to  the  fractional  vegetation  cover  and  soil  moisture  content.  A triangular  shape 
appears  by  plotting  various  Florida  land  covers  in  the  NDVI  and  the  Diurnal  Surface 
Temperature  Variation  (DSTV)  scattergram.  The  three  comers  of  the  triangle  are 
agriculture,  wetland,  and  urban  for  south  Florida  land  use/cover  types.  By  using  the 
weighting  factors  derived  from  the  NDVI  and  DSTV  triangle  relationship,  the  vegetation 
and  moisture  coefficient  (VMC)  can  be  calculated  for  areas  with  various  vegetation  covers 
and  soil  moisture  contents. 

The  NDV1-DSTV  triangle  algorithm  and  the  subsequently  determined  VMC  were 
proven  capable  of  estimating  ET  on  a regional  scale.  Due  to  its  fewer  ground  parameter 
requirements,  the  method  developed  in  this  study  can  be  used  for  ET  estimation  over 

words,  by  obtaining  coefficients  from  three  homogeneous  land  cover  areas,  one  can  then 


three  cover  types  used  in  this  study  arc  not  available,  with  this  method,  values  from  three 

5.2  Rreommendaiipng  for  Further  Research 

Vegetation  species  mapping  is  possible  with  high  precision  using  spectral  feature 
analysis  of  data  from  hand-held  spectroradiometer.  This  information  can  be  further 
applied  to  vegetation  classification  of  multi-spectral  data  from  next  generation  airborne  or 
space-borne  data. 

Spectral  reflectance  of  crops,  natural  vegetation,  as  well  as  different  soil  types  and 
minerals  of  Florida,  can  be  measured  from  hand-held  spectroradiometer  resulting  in  a 
database  of  reflectance  signatures.  Such  a database  can  be  used  as  a knowledge  bank  to 
study  the  mixed  spectrum  in  satellite  images.  A spectral  unmixing  model  can  then  be 
applied  to  determine  the  components  and  percentage  of  vegetation  species  and  soil 
backgrounds  within  each  satellite  pixel. 

in  order  to  establish  reliable  ratios  of  actual  to  potential  ET,  long  term  ET  study  of 
wetland  vegetation  is  recommended.  These  parameters  should  be  developed  locally  and 
are  important  for  the  accurate  appheauon  of  remote  sensing  methods  for  the  determination 
of  regional  ET  estimation  in  Florida. 


APPENDIX  1 

SPECTRAL  REFLECTANCES  OF  NATURAL 
TARGETS  IN  FLORIDA 
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Figure  Al-4.  Spectral  reflectance  curves  for  eucalyptus  ( Eucalyptus camakhilensis) 
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Wavelength  (nm) 


Spectral  reflectance  curves  for  a single  eucalyptus  (Eucalyptus 
camaldulemis)  leave  Bold  face  line  is  the  mean  reflectance.  The  thin  lines 
are  upper  and  lower  90%  confidence  intervals. 
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Figure  AI-6. 
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Figure  A I -7. 


urves  for  sugarcane  ( Saccharum ).  Bold  face  line  is 
The  thin  lines  are  upper  and  lower  90%  confidence 


Wavelength  (nm) 


Spectral  reflectance  curves  for  saw  palmetto  (Serenoa  repens).  Bold  face 
line  is  the  mean  reflectance.  The  thin  lines  are  upper  and  lower  90% 
confidence  intervals. 
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Wavelength  (nm) 


Figure  Al-9.  Spectral  reflectance  curves  for  catlail.  Bold  face  line  is  the  mean 

reflectance.  The  thin  lines  are  upper  and  lower  90%  confidence  intervals. 
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Figure  Al-10.  Spectral  reflectance  cu 
secumdatum).  Bold  fa 
upper  and  lower  90%  < 
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Figure  Al-ll.  Spectral  reflectance  curves  for  Bibb  lettuce  Bold  face  line  is  the  i 


efleclance. 


Wavelength  (nm) 


Figure  Al-12.  Spectral  reflectance  curves  for  dry  muck  soil.  Bold  face  line  is  the  mean 
reflectance.  The  thin  lines  are  upper  and  lower  90%  confidence  intervals 


APPENDIX  2 

PROGRAM  CODES  OF  NOAA  AVHRR  PREPROCESSING 

This  is  a group  of  programs  written  in  C language  to  perform  the  NOAA  AVHRR 
data  extraction  and  image  processing.  The  names  of  the  subroutines  and  their  functions 
are  listed  below: 

1.  MAIN.CPP:  main  program  to  manage  all  subroutine  calls; 

2.  HRPT.CPP:  subroutine  to  extract  NOAA  AVHRR  raw  image  and 

header  auxiliary  data 

3.  FUP.CPP:  subroutine  to  rotate  a raw  image  file 

4.  SET.CPP.  subroutine  to  set  the  input  and  output  parameters 

5.  ALBEDO.CPP:  subroutine  to  calculate  narrow  band  albedo  in  NOAA 

AVHRR  channel  I and  2 

6.  TEMP.CPP:  subroutine  to  calibrate  radiant  temperatures  from  NOAA 

AVHRR  image  channels  3,  4 and  5 

7.  GEOREF.CPP:  subroutine  to  geo-rectify  NOAA  AVHRR  image 
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II 

//  NOAA/HRPT  DATA  AND  IMAGE  PROCESSING  PROGRAM 


//  INCLUDE  THE  FOLLOWING  FUNCTIONS: 

//  I . READ  PACKED  FILE 

//  2.  ROTATE  IMAGE 

II  3.  ALBEDO  CALIBRATION 

//  4.  TEMPERATURE  CALIBRATION 

//  5.  VEGETATION  INDEX  CALCULATION 

//  6.  FORMAT  CONVERSION 

II  7.  GEOMETRIC  CORRECTION 

II  8.  DOS  SHELL 

//  Program  Code:  C++ 

II  Author:  C.  H.  Tan 
II  Version:  1.03 
II  Date:  Jun  12, 1993 


/*  The  following  is  the  main  program  which  call  subroutines  to  implement  the  data 
extraction  and  image  processing  tasks. 

PROGRAM  MAIN.CPP 


II  define  included  head  files 
((include  <stdio.h> 

((include  <stdlib.h> 
((include  <etype.h> 


//  define  external  subroutines 
extern  void  Read_HRPT(void); 
extern  void  Albedo(void); 
extern  void  Vegetation Jndex(void); 
extern  void  Temperaturefvoid); 
extern  void  Flip(void); 

extern  void  Georefjvoid), 
extern  void  Set_Ftle(char  •); 
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void  Main_Help(void); 
void  Select_Command(void); 

struct  Header  { 

ini  ii,  II,  ie.  Ie,  nc,  bt; 


FILE  •infI,*inS,*infl,*oufl,*ouG,*ouft; 

char  infl_namc(20],inf2_namc[20],inft_name[20]; 

charoufl_name(20],ouf2_name[20],ouft_naine[20]; 


void  main(void)  { 
clrscrO; 

printli“\n\l Welcome  lo  NOAA/HRPT  Data  and  Image  Processing  ProgramW); 
prinlfl;"\n\t\t\tVeiaion  1.04.  Aug.  1993W); 

MainHelpO; 

Select_Command(); 


void  Selecl_Command(void)  { 
char  *p.cmd[20); 

printfT\nNOAA->-); 

p=gets(cmd), 

forO"0;i<20;i++)  | cmd(i]=lolower(cmd[i]);  ( 

i((cmd[0]='h'&&cmd[  I ]— e'&&cmd(2]=T&&cmd[3]='p') 
Main_HeIpO; 

else  iRcrndfO]— d'&&cmd[  I ]=’o'&&cmd[2J— s')  ( 
if(*(p+3)=^0’)  system("cominand.conr); 
else  system(p+4); ) 

else  i(|cmd[0)='='h'&&cmd[  I ]='r‘&&cmd[2]=’p'&£:cmd[3]’='t') 
Rcad_HRPTO; 

else  itTcmd[0]=T&&cmd[l)=T&&cmd[2]=T&&cind[3]— p') 
FlipO; 

else  if(cmd[0]— a’&&cmd[  1 ]=T&&cmd[2]='b,&&cmd[3)==,d') 
AIbedo(); 

else  if(cmd[0]— t'&&cmdf  I ]“'e'&&cmd[2]=’m'&&cmd[3]='p') 


else  iHcmd[0]=V&&cmd[  I ]='e'&&cmd(2]='g,&&cmd(3]=,i') 
Vegetation JndexO; 

else  iHcmd[0]='g'&&cmd|  I ]=,e'&&cmd[2]='o'&&cmd[3]“'r') 
GeoreH); 

else  itlcmd[0]--,e’&&cmd[  1 ]=T&&cmd[2]=='a'&&cmd[3]=='s') 

ELAS  FormatO; 

else  iHcrad[0]— q'||cmd(OJ“’e'&&cmd|  I )«-'x'||cmd[0)=W) 
exit(0); 

else  printlf"Unknown  command  ....  %s",cmd); 


Iain_Help(void)  ( 

printf("\ii\tHRPT  : read  raw  HRPT  data"); 
pnntH"\n\lFLIP  : flip  image  file"); 
printH"\n\tALBD  : albedo  calculation"); 
printfl"\n\tTEMP  : temperature  calculation'); 
printfl"\n\tVEGI  : vegetation  index  calculation"); 
printl("\n\tELAS  ; convert  tile  to  ELAS  format'1); 
printl("\n\tGEOR  : germetric  correction"); 
printf("\n\tDOS  <cmd>  : dos  command  <cmd>"); 
printft”\n\tHELP  : this  messaged"); 
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SUBROUTINE  TO  EXTRACT  NOAA  HRPT  RAW  IMAGE  AND  HEADER 
AUXILIARY  DATA 


II  PROGRAM  HRPT.CPP 


"define  swap(a,b)  (unsigned  char  c;c-a.a=b,b=c,| 


void  Read_HRPT(void); 

void  HRPT_Header(void); 

extern  void  Set_File(char  •); 

extern  void  Set_Parameter(char  *); 

extern  struct  Header  (int  il,  II,  ie,  le,  nc,  bt;)  H; 

extern  FILE  •infl,,in(2,*inft,*oufl,,ouG,*ouft; 


void  Read_HRPT(void)  { 
int  ij,k,l, line, element; 
long  int  skip; 

int  value[10242],ch(5][2048); 

struct  Hd  ( 

unsigned  d3: 10; 
unsigned  d2: 10; 
unsigned  dl:  10; 
unsigned  null:  2; 

)h; 

union  DN  ( 

long  inti; 

unsigned  char  ch[4); 
struct  Hd  h; 

) dn[34l4J; 

charcmd[20),,p,*endptr,path[]="NOAA-HRPT‘; 
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for(N);i<2(y++)  { cmd[i]=tolower(cmd[i]); ) 
nd[0]=' s'&&cmd[  I ]— f)  Set  Filefpath), 
iHcmd[0]— s'&4tcmd[  1 J='p')  Set_Parameter(path); 
if|cmd[0]='s'&&cmd[l]='k'&&cmd[21=='i'&&cnid[3]='p')( 
skip=strtol(p+5,&endptr,  1 0); 
print!!"SKIP  = %ld\n",skip); } 

else  ificmd[0]=='h'&&cmd[l)— e'&&cmd[2]=T&&cmd(3]='p')  j 
printf)  "\n\tSF  : set  input/output  files"); 

set  parameters"); 

<i> : skip  i byte  in  input  file"); 
printH"\n\tRUN  : run'); 
print!! "\n\tEX  : eat"); 

print^"\n\tHELP  : this  message"); 


line=H.n-H.il+l; 

clemcnt=H.Ie-H.ie*l; 
fseek(infl  ,skip,SEEK_SET); 
fseek(infl,14800L*(H.il-l),SEEK  CUR); 
for(l=0;l<line;l++)  ( 
printl!"\033[4D%3.0!%%”,l*100./line); 
HRPT_Header();  //readHRPThe 

fread(dn,4,3414,infl);  //  34 1 4 4 byte  di 

fseek(infl  ,696L,SEEK_CUR);  //  696  empty  spa 

fbr(i=0 J=0;i<34 1 4;i++>  { 

swap(dn[i].ch[3],dn[i].ch[0]); 
swap(dn[i).ch[2j,dn[i].clii  I ]); 
value[j]=dnp].h.dl; 
value[i+l)=dn[i).h.d2; 
value(j+2]=dn[i].li.d3; 
j+-3; } 

for(H>;i<5;i++)  { 115  channels 


//2048  samples/channel 
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k+-5; ) } 
for(iiO;i<5;i++)  { 

for(j=0;j<element;j++)  ch[i][j]=ch[i][j+H.ie-l); 
lwrite(ch[i],2,element,oufl); ) 


) fclose(oun);  fdose(ouf2), 

else  if(cmd[0]=,e'&&cmd[  1 ]“V||cmd[0]=W)  return; 
else  printifUnknow  command  ...  %s",cmd); 


) //program 


void  HRPTHeaderfvoid)  { 

im  i.P.k: 

long  int  slopc[5],imcept[5]; 
unsigned  char  hcader[448]; 
struct  HRPT  HEAD  { 

int  year, 
int  julian; 
float  time; 
float  S[5]; 
float  I[5]; 

int  solar[Slj; 
float  lon[SI]; 
float  lat[5l]; 


fread/header,  1 ,448,  inf] ); 

head.line=header|0J*256+headcr[l]; 
hcad.ycar=header(2]»  I ; 

head.julian=(headcrI2J%2),256+header[3); 

head.time=<header[41%8),2S6.,256.*256.+header|5]*256.*256. 

+header[6]-256.+header(7]; 

fbrfi=0,p=(y<5;i++)  { 

slope[i]  = header[p+ 1 2]"256L*256L,256L+header[p+- 1 3 J*256L*256L 
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+header[p+14),256L+header[p+15]; 
intceptp]  = hcader[p+161"256L*256L,256L+headerfp+17]"256L,256L 
♦headerf  p+ 1 8]  •256L+header[p+ 19]; 
head.  S[i]=(floal)slope[i)/l  07374 1824.; 
head.I[iHfloal)imccptp)/4194304.; 
p4=8;  ) 

head.meaning=header[S2], 
for(i=0;i<bead.meaning;i-H-)  { 

head.solarp]=headerp+53]/2; 

head.lalpHheaderp,4+104]*256.+hcaderIi*4+105]yi28.; 

head.lonp]=(65536.-(headerpM+106]*256.+headcrp*4-H07J))/128.; 


fwrite(&head,sizeoHhead),  1 ,ouf2); 
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SUBROUTINE  TO  ROTATE  A RAW  IMAGE  FILE 

/*  The  raw  image  need  to  be  rotated  1 80  degree  in  the  cases  of  ascending  satellite 

ground  track. 


//  PROGRAM  FLIP.CPP 

//define  included  head  files 
(/include  <stdio.h> 

//include  <stdlibh> 

//include  <ctype.h> 


void  Flip(void); 

extern  void  Set_File(char  •); 

extern  void  Set_Parameter(char  •); 

extern  struct  Header  |int  0,  II,  ie,  le.  nc,  bt;)  H; 

extern  FILE  *infl,*inf2,*inft,*oufl,*ouf2,*ouft; 

void  Flip(void)  { 


int  ch2[5000],flip2[5000]; 
unsigned  char  chi  [SOOO],flipI[SOOO]; 
long  int  line,elenicnt,skip,byte; 
char  cmd[20],path[KNOAA-FLIP"; 


printtf"\n%s->",path|; 

gets(cmd); 

for(i=0;i<20;it+)  { cntd[i]=tolower(cmd(i)); ) 

i«cmd[OJ=’s'&&cmd[l]='P)Set_File(path); 

else  if|cmd[0|— 5'&&cmd[  1 ]=' p'7  Set_Parameter(path); 

else  if(cmd[0]='h'&&ctnd[l)='e'&&cmd[2)-=T&&cmd|3]='p')  ( 


else  if[cmd[0]=V&&cmd[  I ]='u,&&cmd[2]=,n')  ( 
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printfCrunning ....  "); 
line=H.ll-H.i!+l; 
dement=(H.Ie-H.ie+l)*H.nc; 
byte=H.bt; 
forfi=0;i<line;i++)  ( 

primH'\033[4D%3.0B4%",i*IOO./lin«); 
skip=(line- 1 L-(long)i)*demem'byte; 
fseek(infl  ,skip.SEEK_SET); 
if(byte=2L)  ( 

fread(ch2,byte.clcmcnl,infl ); 
for(j=Oj<eleiiientj++) 

flip2[j]=ch2[element-j]; 
(wrilc<nip2,byle,element,oufl ); ) 
ifibylc=IL)  ( 

fread(ch  1 ,byle,elemcnl,infl ); 
for(j=Oj<demcnij'H-) 

flipl[j]=chl[elemenl-j]; 
fwrite(flip  I ,byte,element,oufl ); ) 
} fdose(oufl),  fdosefmfl): 


else  iHcmd[0]— e’&&cmd[  I ]=Y||cmd[0]“W)  rMum; 
else  prinlffUnknow  command  ...  %s",cmd); 


SUBROUTINE  TO  SET  THE  INPUT  AND  OUTPUT  PARAMETERS 


Input  and  output  parameters  including  file  names,  line  and  column  numbers,  < 
channels,  output  byte  per  pixel. 


II  PROGRAM  SET.CPP 

//  define  included  head  files 
((include  <stdio.h> 

((include  <stdlib.h> 

((include  <string.h> 

((include  <ctype.h> 

void  Set _Fde(char  pathO); 

void  Sct_Parameter(char  path[]); 

extern  strnct  Header  {int  il,  II,  ie,  le,  nc,  bt;>  H; 

extern  FILE  •inn,,inf2,*inft,,oufl.,ouf2,,ouft; 

extern  char  infl_name[20],inG_name|20I.inIi_namc[20]; 

extern  char  oufl_name[20],ouf2_nameI20|,ou!t_namc[20]; 

void  Set_File(char  pathQ)  ( 
inti,f; 

char  cmd[20],*p; 

foK;;)  f 

printfi"\n%s-FILE->",path); 

p=gets<cmd); 

for(i=0;i<20;i++)  { cmd[i)=tolower(cmd[i]); ) 
iHcmd[0]=’i'&&cmd[  I ]— n'&&cmd[2]='f&&cmd[3]—  1 ')  ( 
ifl{infl=fopen(p+5,"ib"))=NULL) 
printHInCannot  open  file !"); 

fbr(i=0;i<I0;i++)  infl  name[i]=toupper(cmdp+5J);)) 
iftcmd[0]— f&&cmd[l)— n'&&cmd[2]='f&&cmd[3]— 2')  ( 
iHCtnG=fopen(p+S,-rb"))=NULL) 
printir'nCannot  open  file !“); 
else  (for(i=0;i<IO;i++)  inC_name[iHoupper(cmdIi+5]);}) 
iHcmd[0]=V&&cmd[  1 )='n,&&cmd[2J='f&&cmd[3]=T)  ( 
ifi(inft=fopen(p+5,"rt"))=NULL) 
printfi'InCannot  open  file !"); 
else  (for(i=0;i<IO;i++)  infi_name[i]=toupper<cmd[i+5]);)) 
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ilTcmd[0]=‘o'&&cmd[l]—u'&&cmd[2]=' f&&cmd[3]==T)  { 
oufl=fopen(p+5."wb"); 

foi<i=0;i<10;i++)  oufl  _name[i]=toupper(cmd[i+51);  I 
if(cmd[0]=,o,&&cmd[  1 ]='u'&&cmd[2]=,f&&cmd[3]='2•)  { 
ou(2=fopen(p+5,"wb"); 

for(i=0;i<10;i++)ouf2_nanie[i]=toupper(cmd[i+5]);) 
ificmd(0]='o'&&cmd[  1 ]=,u,&&cmd[2]='f&&cmd[3]=,t')  ( 
ouft=fopen(p+5,-wf); 

for(i=0;i<10;i++)ouft_namc[ij=loupper(crad[i«]);) 
ificmd[0]=’c'&&cmd[l]=='f)  { 

iH!strcmpi(p+3,-infl"))  infl=W; 
ifT!slrcmpi(p+3,-inl2"))  inE=W; 
i«!strempi(p+3,"inft"))  inft-W; 
if(!strcmpi(p+3,"oufl "))  oufl=W; 
il(!strcmpi(p+3,"ouf2"))  ouf2=W. 
ifi!strcmpi(p+3.-ouft-))ouft=W;) 
if(cmd[0]==T&&cmd[  I )=?)  ( 
f=0; 

iHinfl  1-NULL)  {prinllfVtINFl  = %s\n\infl_name);ft+;} 
if{inf2!=NULL)  (printlC\tINF2  = %s\n",inl2_name),f++,) 
ifiinftl=NULL)  { prim(("\lINFT  = %s\n',inft_name);ff+; ) 
iftoufl  !=NULL)  (printf("\lOUFI  = %s\n-',oufl_name);f+->-,| 
iRoufi!=NULL)  jprimlT\lOUF2  = %s\n",ouC_namc);fn-;) 
il(oufl!=NULL)  jprinttCVtOUFT  = %s\n",ouft_name);ft-n) 
if(fi=0)  primfTMNo  opened  file  !\n"); ) 
iRcmd[0]='h'&&cmd[  1 ]=’e'&&cmd[2]=T&&cmd[3]=='p')  ( 
printfl>\t\tLF  : list  opened  file  name"); 
printfC\n\t\tINFl  <fii> : open  first  input  file  <&>"); 
printfi"\n\t\tINF2  <fn>  : open  second  input  file  <fi»"); 
printH"\n\t\tlNFt  <fit> : open  text  input  file  <&>■); 
printR"\nVt\tOUFI  <fh> : open  first  output  file  <fh>*); 
printfi"\n\t\tOUF2  <fn>  : open  second  output  file  <&>”); 
printfC\n\t\tOUFT  <(h> : open  text  output  file  <fh>“); 
printfi"\n\t\tCF  <fs>  : close  file  stream  <&>"); 
printfi"\n\t\tEX  : exit"); 
printfi"\n\t\tHELP  : this  messaged"). ) 
ificmd[0)=W||(cmd[0]=’e'&&cmd[  1 ]='x'))  return; 
else  printffUnknown  command  %s",cmd); 


void  Set_Parameter(char  path[]) 
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char  cmd[40),*p,*endptr; 


M;;) 


prinlH"\n%s-PARM->",path); 

p=gets(cmd); 

for(i=0;i<40;i++)  { cmd[i]=tolower(cmd[i]), ) 
fbr(i=0;i<40;i+'t-)  { 

if(cmd[i]=“,i'&&cmd[i+l)=T)  H.il=strtol(pti*3,&endptr.  10); 
itf  cmd[i]=T&&cmd[i+ 1 ]=T)  Hll=strtoI(p+i-*-3,&endptr,  10); 
iHcrndji)— i'&&cmd[i+ 1 ]— e')  H.ie=stitol(p+i+3,&endptr,  1 0); 
iHcrad[i)=Tdt&cmd[i+lj— e')  H.le=stnol(p+i+3,&endpIr,  10); 
iHcmd[ij='n'&&cmd[i+ 1 )='c')  Hnc=strtol(p+i+3,&endptr.  1 0); 
Hcmd[i]=V&&cmd[i+lj='l')  H.bt=slrtol(p+i+3,&endptr,  10), ) 

il(cmd(0]='h'&&cmd[  1 ]— e'&&cmd(2]=T&&cmd[3]='p')  ( 


prinl(!"\n\l\!BT  : byte  per  pixel"); 
printfl"\n\t\lEX : exitin'); ) 


i((cmd[0]— r&&cmd[ I J—p1)  { 

print(("\n\t\tIL  = %6d\lLL  = %6d",H.il,H.ll); 
printH"\n\t\tlE  = %6d\tLE  = %6d",H.ie.H.le); 
print(T\n\t\tNC  = %6d\tBT  = %6d\n",H.nc,H.bt); ) 
il((cmd[0]— e'&&cntd[  1 ]='x')||cmd[0]=W)  return; 
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SUBROUTINE  TO  CALCULATE  NARROW  BAND  ALBEDO  IN  NOAA  HRPT 
CHANNEL!  AND 2 


II  PROGRAM  ALBEDO.CPP 

//  define  included  head  files 
((include  <stdio.h> 

((include  <ctype.h> 


void  Albedo(void); 

extern  void  Set_Fde(cliar  •); 

extern  void  Set_Parameter(char  •); 

extern  struct  Header  {int  il,  II,  ie,  le,  nc,  bt;J  H; 

extern  FILE  •infl.’inC.'inft.'oufl  ,'ouI2.-ouft; 

extern  char  ouG_name[20]; 


void  Albedo(void)  { 
float  a!bedo[2][2048], 
int  ch[2][2048); 
unsigned  char  albout[2][2048]; 
int  ij,l,line,element; 

char  cmd[20],*p,path[]=‘NOAA-ALBD'; 
struct  HRPT  HEAD  ( 


float  S(5]! 
float  I[5]; 

int  meaning; 
intsolar[51]; 
float  lon[5I]; 
float  lat[SI); 


printH"\n%s->“,path); 

p=gets(cmd); 

for(ii=0;i<20;rH-)  { cmd[i|=tolower(cmd[i]), } 
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SSSSi,? 


for(l=0;l<line;l++)  { 

prints '\033(4D%3.0l'/o%M-IO0./linc); 


— m’&&cmd[3]='p')  { 


^^o[0|.4,da™m,°uC); 

iH*(p+5)=T)  ( ' 

fprinIftouft,"\lPixel  Chi  Ch2  ... " 

"line  #%4d\n\H.i|-H); 
for(j=Oj<eleniemJ++) 
fprintHoufl,‘\l%4d  %7.3f  %7.31\n", 
j+H  ie.albedo[0)[j),  albedo!  1)01); ) 
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else  { 

primf(-\iPixel  Chi  Ch2  ...  lineWMdW, 
H.il+I); 

forQ=Oa<elemeiitj-H-) 
primfCIWMd  %7.3f  %7.31\n', 

j+H.ie,aIbedo[0)[j],albedo[ !][)]); } 


else  iHcmd[0]— e'&&cmd[l]=\'||cmd[0]='W)  return; 
else  printH'Unknow  command  ...  %s",cmd); 

}//for  loop 
(//program 
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SUBROUTINE  TO  CALIBRATE  RADIANT  TEMPERATURES  FROM  NOAA  HRPT 
IMAGE  CHANNELS  3, 4 AND  5 


//  PROGRAM  TEMP.CPP 


//  define  included  head  files 
((include  <stdio.h> 

((include  <stdlib.h> 

((include  <5tring.h> 

((include  <c!ypeh> 

((include  <conio.h> 

((include  <math.h> 

void  Temperature! void); 

extern  void  Set_Filc(char  •); 

extern  void  Sct_Parameter(char  •); 

extern  struct  Header  {int  i!  II,  ie,  le,  nc,  bt,(  H; 

extern  FILE  •infl,,inf2,,infi,'oufl,"ou(2.’ou(l, 

extern  char  ouf2_name[20]; 


void  Temperaturc(void)  { 

float  Temp[3][2048]; 
unsigned  char  Temp_8[3][2048]; 
int  ch[3][2048]; 

char  cmd[20],’p,flag=’0,,path[]="NOAA-TMPT", 
float  cl=I . I9I0659E-5.C2-1 .438833; 
float  e,tl,wave[3); 
struct  HRPTHEAD  { 

float  time; 
float  S[5); 
float  I[S]; 

int  solar[S  lj; 
float  lon[SI]; 
float  lat[S  1 ]; 


forfc)  ( 


printf["\n%s->",path); 

p-gels(cmd); 

for[i=0;i<20;i-H-)  ( cmd[i)-tolower(cmd[i]); ) 

if(cmd[0]— s'&&cmd[  1 ]=-f)  Set_File(path); 
else  if(cmd[0]— s'&&cmd[  I ]— p')  Sel_Parametcr(palh); 
else  il[cmd[0]— h'&&cmd[  I ]=='c'&&cmd[2]==T&&cmd[3] 
printf[“\n\tSP  : scl  parameters'), 

prinlli"\n\tSF  : set  inpul/oulput  files"); 

printf["\n\tRUN  : mn">; 
print!) "\n\lDUMP  <n>  : n=0  lo  screen,  n=l  to  file"); 


NF2=head  file, " 

"•OUFl=output  (1  byte)\n"); 
printH"  0UF2=output  (2  bytes),  OUFT=text  dump" 

" (•  required  filc)\n");J 

else  if)cmd[0]=V&&cmd[  I ]='u'&&cmd[2)='n')  ( 
prinlfi'running ....  ■); 

wave[0)=267l.4; 

wave[  I ]=927.8;  //  assume  temp  275-320  in  the  beginning 

wave[2]=838.4; 

forO=0;l<H.II-H.il+l;l++)  ( 

printf("\O33[4D%3  0f%%",|-IO0./(H.II-H.il+l)); 
fseek(infl,4L*(H.le-H,ie+l),SEEK_CUR); 
for(i=0;i<3;i++)  { 

fread(ch[i],iH.le-H.ie+l,infl); 
fread(&head,sizeo ([struct  HRPT  HEAD),  I,inl2), 
for(j=0J<H.le-H.ie+ 1 j++)  ( 


if(Temp|i][j]>225.&&Temp[i)[jJ<=27S  &&flagl=T)  ( 
wave|0]=2668.I5; 
wave[l]=927.34; 
wave[2]=838.08; 


e=head.S[i+2]*ch[i]D]+l>ead.I[i+2); 

tl=l.+cl*wave[i]“wave[i]*wave(ij/e; 


Temp[i)lj]=c2"wave[iJ/log(tl)-273.  i 5; 


flag=T; 
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goto  count jigain; 

) 

iflTemp[i][j]>180.&&Temp[i][j]<=225.&&flag!=‘2')  { 
wave[0]=2663.5; 
wave[l  ]=926.8; 
wave(2]=837.75; 
flag=’2'; 

goto  count_again; 

) 

ifCTcmp[i]D]>41-)  Temp[i]Q]=41.; 
iflTemp[i][j)<-10.)  Temp[i][j]— 10.; 
Temp_8[i][j]=(unsigned  char)((Temp[i](j]+10.)*S.); 

fwrile(Tenip_8[il,  I ,Hle-Hie+ 1 ,oufl ); 
fwrile(Temp[i],4,H.le-H.ie+ 1 ,ou!2); ) 

) fc!ose(oufl);  fclose(ou(2); 

) 

else  iflcmd[0]='d'&&cmd[l]='u'&&cmd[21— 'm'&&cmd[3]=,p')  1 
oul2=fopen(ouf2_namc,"rb"); 
for(I=0;l<H.II-H.il+ 1 ;l++)  ( 

(read(Temp(0J,4,H.le-H.ie*l,oul2); 
(read(Temp[l),4,H.le-H.ie+!,ouG); 
fread(Templ2].4.H.le-H.ie+l.ouG); 
in*(p+s)=-io  { 

fprintfioufl,"\tPixel  Ch3  Ch4  Ch5" 

" ..  line«%4d\n",H.il+l); 
for(j=0  j<H.le-H.ie+ 1 j++) 
fprintf(ouft,"\t%4d  %8  3f  %8,3f  %8.3f\n" 
j+H.ie,Temp[0][j],Tetnp[l][j],Temp[2][j]l; ) 

else{ 

printH"\tPixel  Ch3  Ch4  CHS  ...  ” 

"line  #%4d\n",H.il+l); 
fot(i=Oj<H.le-H.ie+l  j++) 
printf("\t%4d  %8.3f%8.3f  %8.3ftn\ 

j+H.ie,Temp[0][j],Temp[l]Q],Temp[2)[j]); ) 


else  if(cmd[01='e'&&cmd[l]=V||cmd[0]=W)  return; 
else printlTUnknow command  ...  %s”,cmd); 


SUBROUTINE  TO  GEO-RECTIFY  NOAA  AVHRR  IMAGE 


subroutine  to  geo-rectify  a image  with  the  Earth  Locations  extracted  from  the 
HRPT  header.  Also  include  subroutines  to  convert  from  Lat/Lon  coordinate  to 
UTM,  affine  transformation,  and  resampling 


PROGRAM  GEOREF.CI 


void  Geo  reR  void); 
double  UTM  N,UTM_E; 
double  X[6); 
int  I jjeo.e _geo; 
long  res_x,res_y; 

void  LL_to_UTM(double,double); 

void  Control_Point(void); 

void  Affine(void). 

void  Resample(void); 

void  invert(double  A[6][6] ), 

unsigned  char  Read_DN(int,int); 

extern  void  Set_Filc(char  *); 

extern  void  SetParameleRchar  •); 

extern  struct  Header  (int  il,  II,  ie,  le,  nc,  bt;)  H; 

extern  FILE  •infl,*inl2,*inft,*oufl,*ouf2,“oufl; 

void  GcoreRvoid)  { 

char  cmd[20),  *p,  •endptr,path[]=,NOAA-GEOR-; 


fort;)  { 

printR"\n%s->",path); 


p=gets(cmd); 

for(i=0,i<20;i++)  { cmd[i]=iolower(cmd[i]), ) 
if{cmd[0)— s'&&cmd[IJ='P)  Set_fiie(path); 
else  ificmd[0]='s'&&cmd[  1 ]— p1)  Set  Parameter(path); 
else  if(cmd[0]— g'&&cmd[  I |=='c'&&cmd[2)“'p'&&cmd[3J=T)| 
Control  PointO; ) 

else  if?cmd[0]='iJ&&cmd[  I )=’e'&&cmd[2|“'s,&&Qmd[3]='m')| 
sample=cmd[5); 
iRsample— !•) 

printlCResampling  method  = Nearest  neighborin'); 
else  il(sample='2') 

printlCResampling  method  = Bilinear  interpolation\n"); 
elseiHsample=T) 

printlCResampling  method  = Cubic  convolution^"); 
printICUnknown  resampling  method  l\n");  ) 


else  if(cmd[0)=y&&cmd[  1 ]— s')! 
res_y=strtol(p+3,&endptr,  10); 

else  iHcmd|0]='h,&&cmdll]■='e,&&cmd[2)=T&&cmd|3]=,p,)  ( 


printH"\n\t  1=3  cubic  convolution"); 

printlC\n\tRUN  : run"); 
printfl;"\n\tEX  : exit"); 

printfl>\tHELP  : this  message"); 
printH"\n\nHint:  *INFl=image  file.  •INFT=control  point  " 
"file"); 

printH"\n\t"OUFl=output  image  file  (*  required  file)\n"); 
printH"\tfor  GCPF  modual  sel  INF2=header  file,  OUFT=control" 
" point  fileVn"); 


else  iflcmd[0]='x'&&cmd[  I ]=  s ) ( 
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else  ittcmd[0]='r,&&cmd[  1 )='u'&&cmd[2]='n')  { 
Affinef); 

Rcsample(); 

) 

else  iftcmd(0]='e'&&cnid[l]=Y||cmd[0]='\0')  return 
else  prirnffUnknow  command  ...  %s",cmd); 

j //for  loop 
) //program 


void  ControlPoinl(void)  ( 
intijj; 

struct  HRPT  HEAD  ( 


float  time; 
float  S[5]; 
float  1[5J; 

int  solar[S  1 ]; 
float  lon[SI]; 
float  lat[S  1 ]; 


FILE  *utmfile; 

int  fcc,lcc,raw_re.raw_le1f; 


l_geo=0; 
e _geo=0; 

printf('\nFirst  element  in  raw  image  ” "); 
scanff%d',&raw_fe), 

scanff%d_.&rawje); 

printfJ"\nUpper-left  comer  UTM  easting  cord.  = "); 
5canff"%f',&min_e); 

printfi*\nUppcr-left  comer  UTM  northing  cord.  = ”); 
scanft"%F,&max_n); 


fec=(raw_fe-25  J/40+ 1 ; //first  ele 


lec-(raw_le-25)/40;  //last  element  control  point 
fbt<i=0;i<H.II;i+=40)  { 

fread(&head,sizcoRstruct  HRPT_HEAD),l,inE); 
fscck(inf2,39L"sizeoHstruct  HRPT_HEAD),SEEK_CUR); 
fbr(j=fcoj<=leca++)  I 

LL_to_UTM((double)head.lat[i],(double)head.lon|j]); 
iH*40+25-raw  fe+l; 

//  min_e=min_e<UTM_E  ? min_e  : UTM_E; 

//  max  n=max_n>UTM_N  ? max_n : UTM_N; 

fprintRutmfile,"%3d  °/=3d  %8.3f  %8.3f  %13.3f%13.3fW7 
head.lat[j],hcad-lon[i),UTM  E.UTM  N); 

} 


rewind(utmfile); 

do{ 

P=fscanHutmfilc,'%d%d%B4P/ol%f'.&i,&j,&t,&t.&e,&n) 

new_l=(max_n-n)/res_y; 

new_e-<e-min_eVres_x; 

1 _geo=l _geo>newJ  ? l_geo:newJ; 
c_geo=c_geo>new_e  ? e_geo:new_e; 
if(fl=EOF) 

fprintf(ouft,"%3d.O  %3d.O  %6.2f  %6.2i\n".ij,new_l,new 
( while(fl=EOF); 

printfT\tGEO-CORRECTED  FILE:\n"); 
printH"\t\t\tLine  = %4d\n",l_geo), 
printff\t\t\tElement  = %4d\n",e_geo); 
fclosefouft);  fc!ose(utmfile); 


void  LL_to  UTM(double  lat,  double  Ion)  ( 
intzone=17; 

double  a=6378206.4,b=6356583.8,k0=0.9996,pi=3.  141592653589793; 


nl,n2,n3,n4,n5,pl,p2,p3,p4,pS,p6,e2,ep2,ep4; 

ss  1 ,ss2,ss3,ss4,ss5,ss6,sp  I ,sp2,epl  ,cp2,cp3,cp4,cp5,tp  1 ,tp2,tp4; 

Ap,Bp,Cp,Dp,Ep,S,A6,B5, 1,11, 111,1  V,V; 


rpd=pi/180.; 

nl=(a-bV(a+b); 

n2=nl*nl; 


ss  1 =sin<  1 ./3600. *rpd);  II a nO") 
ss2=5sl*ssl; 
ss3=ss2*ssl; 
SS4=5S3*55l; 


sp  1 “sin(phi*rpd);  II  sin(phi) 

sp2“spl*spl; 
cpl=cos(phi,rpd);  //  cos(phi) 
cp2=cpl*cpl; 
cp3=cp2*cpl; 
cp4-cp3'cpl, 
cp5=cp4*cpl; 

lpl=spl/cpl,  //tanfphi) 

lp2=tpl*tpl; 
tp4=tp2*tp2; 
e2=(a,a-b*b)/(a-a), 
ep2=c2/(l.-e2); 
ep4=ep2'ep2; 
nu=a/sqit(l  -e2*sp2); 

Ap=a*(l.-nl+5./4.*(n2-n3)+81./64.“(n4-n5)); 

Bp=3./2.*a*(nl-n2+7./8  *(n3-n4)+55./64.,n5); 
Cp=15V16.,a,(n2-n3+3./4.*(n4-n5)); 

Dp-35./48.*a,(n3-n4+l  l./16.*n5); 

Ep=3 1 5/5 1.*a*(n4-n5); 

S=Ap’phi,60.,.000290888208666-Bp*sin(2.*phi*rpd)+Cp*sin(4.,phi,rpd) 

-Dp’sinfb.’phi'rpdj+Ep’sinfS.'phi’rpd); 


I = S*kO; 

II  =nu*spl,cpl,ss2*kO/2.*l.E8; 

in  = 5s4,nu,spl’cp3*(5.-tp2+9.,ep2,cp2+4.,ep4,cp4)*k0/24.,I.E16; 

IV  = nu*cpl*ssl*kO*l.E4; 

V - s$3*nu*cp3,(l.-ip2+ep2*cp2)*kO/6.*l.EI2; 

cenlon-(30-/one)"6.t3.i 

pl“fabs(cenlon-!on)*.36; 

p2=pl*pl; 

p3=p2*pl; 

p4=p3*pl; 
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p5=p4*pl; 

p6=p5*pl; 

A6  = p6*ss6*nu*5pl,cp5*(6L-58,*tp2+tp4+270.*ep2*cp2-330.*ep2*sp2) 
•k0/720.,l.E24; 

BS  = p5*Bs5*nu*cp5*(5.-18,*tp2+tp4+l4.*ep2*cp2-58.*ep2*sp2) 
*k0/120.*l.E20; 

UTMN  = 1+H*p2+Ill*p4+A6; 

it(ccnIon<lon) 

UTMJE  = 500000.-(IV,pl+V*p3+B5); 

UTM_E  = SOOOOO.+(IV*pl+V*p3+B5); 


void  Affine(void)  { 
int  fread.row,col.k,poinl=0.i; 

double  random_x,random_y, control  X,control_Y,S0,Sum  error=0.; 
double  AI[2][6],ATA[6J[6),ATL[6]; 
double  LI[2],VI[2],S[6]; 


for(iow=0;row<6;row++)  { 
ATL[row]=X[row]=0. , 

for(col=0;col<6;col++)  ( ATA[rowJ[col)=0.; ) 

for(row=0;row<2;row++)  { 

LI[row]=VI[row]-0., 

for(col=0;col<6;col++)  ( AI[row][col]=0.; ) 


/•  form  Al.  LI,  ATA.  ATL  metrics  •/ 


ftead=f5canftinft,"%IIV.I(%II%ir,&control_X,&controi_Y, 
& random  x,&random_y); 

ifffread=EOF)  break. 


AI[0][0]=A1[1][3J=1.; 
AJ{0][1  ]=AI[  1 ][4]=randon 
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AI[0)[2)=AI[  1 ]£S]=iandom_y; 
AI[0][3]=AI[0][4]=AI[0][5]=0.; 
AI[1][0]=A1[I][1J=AI[1][2]=0.; 
LI[0]=control_X; 

LI[l]=control_Y; 

for(row=0,row<6,row++)  { 
for(col=0;col<6;col++)  ( 
for(k=0;k<=I;k-t-t-)  { 

ATA(row][col]+=AI[k][row]*Al[k][colJ; 


for(row“0;rov.<6,row-f-i ) ( 

fbr(k=0;k<=l;k++)  ( 

ATL(i©w]+=AI[k][row]*LI[k]; 


) while(fread!=EOF); 

/•  invert  ATA,  multiply  inverted  ATA  by  ATL  to  get  X V 
invert(ATA); 

for  (row=0;row<6;row++)  ( 

X[row]=0.; 

for  (col=0;col<6;co!++)  ( 

X[row)+=ATA[row][col],ATL[coI]; 


rewind(inft): 

do( 

fread=fscanR)nft,"%ll%ll%ll%ir,&controi_X,&conlrol_Y, 

&random_x,&random_y). 


AI[0][0)-AI[I][3]=I; 

Al[0][l  ]=AI[  I ][4]=random_x; 
AI[0][2]=AI[  1 ][5)=random_y; 
AI[0][3]=AI[OJ[4]=AI[OJ[5J=O.; 
AI[1][0]=AI[I][I]=AI[I][2]=0; 
U[0]«controi_X; 


LI[l]=control_Y; 


for  (k=0;k<=1  ;k++)  ( 

VI[k]=0.; 

for  (col=0;col<6;coH-+)  ( 

VI[k]4=AI[k][col]*X[col); ) 

VI[kj-=LI[k); 

) 

Sum_erTor+=(VI[0]*VI[0J+VI(  1 ]•' VI[  1 ]); 
point++; 

) whilc(fread!=EOF); 

//  calculate  standard  error  of  unit  weight 

S0=sqrt(Sum_error/(double)(2’point-6)); 

//  calculate  standard  errors  of  transformation  parameters 

for  (row=0;row<6;row++)  ( S[row]=SO'sqn(ATA[row][row]); ) 

printfl:"\n\nTRANSFORMATION  PARAMETERS  AND  THEIR  STANDARD 
ERRORS:\n\n"); 

printH"\tMapping  Parameter  Standard  ErroiW); 

printti"  \n“); 

for  (i“0;i<3;i++)  { 

printR'\ta%d  =%13.7f\tSa%d  =%10.7f\n".i.X[i],i,S[i]);( 
for  (i=0;i<3;i++)  ( 

printH"\tb%d  =%13.7f\tSb%d  =%10.7ftn\i,X[i+3],i,S[i+3]);) 
printH"\n\nSTANDARD  ERROR  OF  UNIT  WEIGHT  = %10.4f\S0); 


void  invert(doub!e  A[6)[6)) 


for  (k=0;k<6;k++)  { 
for  (j=0 j<6 j++)  { 
iffl  l=k)  ( A[k][jJ/=A[k][k], ) 


A[k][k]=l./A(k](k]; 


for  (i=0;i<6;i++)  { 
if(i  l=k){ 

for  (j=0 j<6 j++)  { 

if  0 1=  k)  ( A[i][j]  -=  A[i][kJ’A[k](jJ; ) 

Ap][k]  •-  -AIk][k]; 

} 


void  ResampleO  { 
ini  i j,i*,iy; 

unsigned  char  dn[2000]; 


printH"\nHOW  MANY  LINES  IN  GEO-CORRECTED  IMAGE  ? "); 
scanR"%d",&l _geo); 

printR-HOW  MANY  ELEMENTS  IN  GEO-CORRECTED  IMAGE  ? "); 
scant)  "?'od’,&e_geo); 
printl!"\nRunning ....  "); 
forfi“0;i<I _geo;i++)  { 

prinllT\033(4D%3OP/o%',i*  100./l_gco), 
forO-Oj<e_gcoj++) ( 


x=X[0]+X[l  ]*i+X[2]*j; 
y”X[3)+X[4]*i+X[5]*j; 

■yt. 

ix=(x-ix)>0.5  ? ix+1  ; ix  ; 
iy=(y-iy>>0.S  ? iy+I  : iy ; 

if(ix>H.ll||ix<H.il||iy>H.le||iy<H.ie)dn[jJ=0; 

else  { 

iftsample— ?•)  { 

printir\nBilinear . . . to  be  linished\n"); 

dsei«sample='3'){ 
prinlf!“\nCubic ..  to  be  iinished\n"); 


dn[j)=Read_DN(ix,iy); 


fwrite<dn,  1 ,c_g  co.oufl  ); 


unsigned  char  Read_DN(int 


nty){ 


unsigned  char  dn; 

ele=H.le-H.ie+l; 
skipe=x*eie*H.bt+y*H.bt; 
fseck(infl  ,skip,SEEK_SET); 

frcad(&dn,  H bl,  I ,infl ); 


SUBROUTINE  FOR  IMAGE  FORMAT  CONVERSION 


/*  The  following  is  a subroutine  to  write  NOAA  HRPT  image  to  ELAS 

PROGRAM  CONV.CPP 


//define  included  head  files 
//include  <stdio.h> 

//include  <ctype.h> 

//define  swap(a,b)  (unsigned  char  c;c=a;a=b;b=c;| 

void  ELAS_Format(void); 

extern  void  Sct  File(char  •); 

extern  void  Set_Parameter(char  "); 

extern  struct  Header  (int  il,  II,  ie,  le,  nc,  bt;}  H; 

extern  FILE  •infl,"inf2,"infl,"oufl,*ouf2,"ouft; 


void  ELAS_Format(void)  { 
int  i,l,line,element; 
union  read{ 

unsigned  char  bt[2]; 

} in2[  10240]; 

unsigned  char  inl[10240]; 

unsigned  char  elas  head[1024]; 

char  cmdpOI.'p.pathlKNOAA-ELAS"; 

FILE  •chl,*ch2,*ch3; 


printH'\n\t  

printlf\n\t  n=2  2 bytes  input,  2 byte  output,"); 
printf("\n\t  n=[]  I byte  input,  1 byte  output"); 


printfC\n\tEX  : exit"); 
printft"\n\lHELP  : this  message"); 

printR"\n\n\nHint:  "INF]  =file  to  be  converted,  *OUFl-output" 
" file"); 

printit"  (*  required  file)\n");) 
else  ifTcmd[01='r,&&cmd[  I ]='u,&&cmd[2]='n')  { 
printH'running ....  "); 

for(i=0;i<1024;i++)  (elas_head[i]=0;) 
elas  head[2)-4;  //  Header  record  length  1024 

elas  head[6)=(H.lc-H.iet-l),H.bt/256; //NBPR 4, 5.6, 7 
elas_head[7]=(H.lc-H.ie+ 1 )*H.bt%256; 
elas_head[  10]=H.il/256; 
elas_head[l  l)=Hil%256, 
elas_head[14)=H  11/256; 
elas  head[15]=H,ll%256; 
elas_head[l8)-H.ie/256; 
elas_head[  1 9]=H.ie%256; 
elas  head[22]=H.le/256; 
elas_head[23|=H.le%256; 
elas_hcad[27J=Hnc; 
elas  head[30]=16; 
elas_head(3lj=225; 


element=H.le-H  ie+T ; 
fwrite(elas_head,  1 . 1 024,oufl ); 
for0=O;l<line;l++)  ( 
printl?"\033[4D%3.0f>/o%",l"100./line), 

i*T*(p+4)='2')  { 

fread(in2,2,H.nc"eIemcnt,infl ); 
for(i=0,i<Hnc'clcmcm;it-+) 

swap(in2[i).bt[l),in2[i].bt[2J); 
fwrite(in2,2,H.nc,element,oun); } 
else  if(*(p+4)=T)  ( 

fread(in2,2,H.nc*element,infl); 

for(i=0,i<H.nc*element,i+*) 

in  1 [i)=(unsigned  char)(in2[i].org/4); 
fivrite(inl,I,H.nc*element,oufl); ) 


/INC  24,25,26,27 
//HEAD  IDENTIFIER 
//(432I),  28-31 


il.oufl);  | 


else  { // 1 byte  in,  1 byte  out 
fread(inl,l,H.nc*elcmcnt, 
fwrite(inl,  I,H.nc*elemcnt 


else  if(cmd[0]=='e'&&cmd[  1 ]=y||cmd[0]=W)  return; 
else  printf("Unknow  command  ...  %s",cmd); 


)//for  loop 
[//program 


APPENDIX  3 


ALGORITHM  AND  SOURCE  CODE  FOR 
TRADITIONAL  ET  CALCULATION 


Modified  FAQ  Pe 


The  original  F AO-24  Penman  combinalion  moihod  (Doorenbos  and  Pniilt,  1977) 
was  revised  by  several  authors  in  the  late  90's.  The  modified  Penman  combination 
equation  has  the  form: 


A ET=-^-(Rn  - G)  + -1—  6.43  1 (ej  - e.  ) 


where  A - latent  heat  of  vaporization  (MJ  kg'') 

A = slope  of  vapor  pressure  curve  (kPa°C‘); 
Y - psychrometric  constant  (kPa°C). 

Rn  = net  radiation  (MJ  m';  d"'); 

G = ground  heat  flux  (MJ  m'?  d'1); 

Uj  = wind  speed  at  2-m  height  (m  s '); 
a*.  b„  = wind  function  coefficients. 


i follow: 


A=2-501-2361xl037\  MJkg  ' 


The  saturation  vapor  pressure  wi 


imatedby: 
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215 


e° 


'16.787-116.9' 
T + 2373 


kPa 


Slope  pf  Ihe  Salutation  Vapor  Preairc  Curve 

The  slope  of  ihe  saturation  vapor  pressure  ct 
the  equation  for  saturation  vapor  pressure: 

A_(7'  + 2373)1  ’k 

where  e°  = saturation  vapor  pressure  (kPa); 

T = air  temperature  (°C). 


The  psychrometric  constant,  y,  represents  a t 
gained  from  air  and  the  sensible  heat  transformed  inti 


where  P = atmospheric  pressure  (kPa); 

cp  = specific  heat  of  moist  air  at  constant  pressure 
X = latent  heat  of  vaporization  (kj  kg'1) 


e,  1.013  kj  kg  ' °C; 
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presented  by  Wright  (1 982) 


where  u.  - wind  speed  at  2 meter  height  (m  s'1); 
D “ day  of  the  year. 


Priestley-Taylor  Method 

The  Priestley-Taylor  method  is  an  empirical  radiation-based  equation  with  the 
following  general  form  (Jensen  et  a!..  1990) 

lET-a-£-(R,-G) 

where  A = Latent  heat  of  vaporization  (MJ  kg'1) 
o - Priestley-Taylor  constant; 

A = slope  of  vapor  pressure  curve  (kPa°C); 

Y = psychrometric  constant  (kPa°C); 


radiation  (MJm’d'1); 
und  heat  flux  (MJ  nr2  d"1). 
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The  Penman-Monteith  combination  method  included  a surface  and  aerodynamic 
resistance  function  to  the  Penman  equation. 


4¥H¥1 


where  p = atmospheric  density  (kg  m'5); 

0,  = specific  heat  of  moist  air  (Id  kg'1  °C'); 

r,  = aerodynamic  resistance  (s  m'1); 

rc  = canopy  resistance  (s  m'1); 

z»  = height  of  the  wind  speed  measurement  (m); 

z,,  = height  of  the  humidity  and  temperature  measurement  (m); 

= roughness  length  for  momentum  transfer  (m); 
z„.  = roughness  length  for  vapor  transfer  (m); 


= wind  speed  at  height  z.  (m  s '); 
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Some  oflhe  parameters  can  be  computed  as  follow: 

7W  = 0.123  h*.  he  is  mean  height  of  crop  canopy  in  mm; 

d -<2«)hc; 


The  FAO  Blaney  Criddle  method  is  based  on  the  general  linear 
between  measured  reference  ET  and  the  Blaney-Criddle/factor  for  m 
that  were  investigated  (Doorenbos  and  Pruitt,  1977), 

ET=a*bf 

/-/>  (0.46  T+ 8.13) 

a ■ 0.0043  RHm-n/N-  1.41 

b - 0.82  - 0.41 « 10'!  1.07  (wW)  + 0.066  Ud 


where  p = mean  daily  percent  of  annual  daytime  hours; 


n/N  = actual  to  possible  sunshine  hours; 
Rh^,  = minimum  relative  humidity; 
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ET  PROGRAM  INPUT  DATA  FORMAT 


Latitude,  Longitude, 


Solar  radiation.  Net  radiation 

Mean  airremperature  difference  (T„  - T,„) 
Maximum  relative  humidity. 

Mean  relative  humidity 
Atmospheric  pressure 

Wind  speed  at  2 meter  height 
Mean  day  time  wind  speed 
Mean  night  time  wind  speed 


a in  the  Blaney-Criddle  equation 

Green  Fetch 

Canopy  resistant 

Height  of  wind  measurement 

Height  of  humidity  and  temperature  measurement 

Crop  height 


[] 


decimal  degree 


MJ  mJ  d'1 


°C 

•C 


[] 


ETEslin 


//  Evapotranspiration  Calculalion  by 

//  I.Panman  Method 

//  2.  Blaney  Criddle  Method 

II  3.  FAO  Pan  Evaporation  Method 

//  ’ 4.  Penman_Monteith  Method 

//  Programer:  Chih-Hung  Tan 

//  Date:  Oct  29, 1997  (version  1.00) 


((include  <stdio.h> 

/(include  <stdlib.h> 

((define  PI  3.14159 

/(define  Cp  1.013/*  specific  heat  of  moist  air,  kJ/kg/C  •/ 

II  declare  global  fiinctions 

void  Read_Data(); 

void  Blaney_Criddle(void); 

void  Pcnman(void); 

void  FAO_Pan(void); 

void  Penman_Monteith(void); 


//  declare  global  variables 

float  Rs,Rn,RHmax,RHmin,RHmean,airt,U2,Ud,Un,P; 
intDOY; 

float  Lat=27,53,EIev=10.,pan=7.5; 

Read_Data(); 

PenmanO; 

Blaney_CriddleO; 

FAO_Pan0; 

Pcnman_MonteithO; 

getchO; 


Did  Read_Data(void)  ( 
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FILE  *data_file; 
char  name[30]; 

printi(”\n\nWeather  station  data  file  = "); 
i(((data_file=fopen(gets(name),,'it"))=NULL)  { 
prinlff'Data  file  not  existlW);  exit(l);} 
fscanfidata_file,"%d%f%P/oP/ot%t%P/oP/oP/«f%f', 

&DOY,&Rs,&Rn,&airt,&RHmax,&RHmin,&RHmean,&P,&U2,&Ud,&Un); 

printfCDOY  Rs  Rn  airt  RHmax  RHmean  RHmin  P U2  Ud 

’printfi" W); 

printfi"%3d%6.2f%6.2f%6.2f%6.2f%6.2f%6.2f%6.2f%6.2f%6.2f%6.2f\n", 

D0Y,Rs,Rn,airt,RHmax,RHmm,RHmean,P,U2,Ud,Un); 


float  LE,e0,Delta,gamma,Wf,VPD,PenET,P_T_ET,G-0.21; 

float  alpha=l. 18;  //locally  verified  by  Albetew  1995 

LE=2.501-0.002361*airt; 

eO  = 0.61 1 *exp((17.27  * airt)/(airt+237.3)); 

VPD=eO*(  1 ,-(RHmax+RHmin)/200.); 

Delta  = (4098*e0y((airt+237.3)*(airt+237.3)); 
gamma  = (Cp  * P)/(0.622*LE)*  IE-3; 

aw=0.4+1.4*exp(-((DOY-173.)/58.)*((DOY-173.)/58.)); 

bw=0.605+0,345*exp(-((DOY-243.)/80.)*((DOY-243.)/80.)); 

Wf=awfbw‘U2. 

i«RHmax>=90.)  RHmax=90.; 
i((RHmin<=30.)  RHmin=30.; 
e=0.892-0.0781*Ud+0.00219*Ud*Rs 
+0.000402*RHmax*Rs+0.000 1 96,Ud/Un*Ud*RHmax*Rs 
+0.00000236*Ud*Ud*RHmax*Rs 
-0.0000086*Ud*Ud/Un/Un*Ud*RHmax 
-0,0000000292*Ud/Un*Ud*Ud*RHmax*RHmax*Rs 
-0.0000161*RHmax*Rs*Rs; 

PenET=((Delta/(Delta+gamma))*(Rn-G)+(gamma'(Delta+gamma))*6.43*WPVPD)/LE; 

P_T_ET=alpha*(DeIta/(Delta+gamma))*(Rn-G)/LE; 

//  printfi"e0=%6.3i\n".e0); 

//  printf("VPD=«/o6.3ftn",VPD); 

//  printl("Delta=%64(Vr, Delta); 


printfl"gamma=%6.4f\n",gamnia); 

prinlH"aw=%6.3f\n",aw); 

prinlH"bw=%6.3(Vn",bw); 

prinlf(”W(=%6.3f\n*,W0; 

primffD/D+r=%6.4(\n",(Della/(Della*gamma))); 

prin!fTc=%6.3l\n",c), 

prmtR"Penman  ETO  = %8  3f  mm/day\iT,PenET); 
printfTP_Tlor  ETO  = %8.3f  mm/day\n\P_T_ET); 


IBIaney  Criddle(void)  { 
floal  a,b,p,n,N,B_C_ET; 
float  LE,della,ws,lolalN=0.; 
float  A,B,RsO; 


for(i=0;i<365;i++)  ( 

delta=0.4093  *sin(2  •PI*(284.+i+l )/365.); 
ws=acos(-tan(Lal’PI/ 1 80  )-tan(delta||, 
N=ws“(24./PI); 
totalN+=N; 


A=3I.54-0.273*Lat+0.00078*Elev; 

BMf.S+^dS'LaH-.OOOdl  *Elev; 

RsO=A+B*cos(2.*Pl*DOY/365,-2.93); 

n-((RsfRsO)-.35)*N/.6l; 

prinlf1"A=%6.2f\n'.B|, 

prin!f|"B=%62f\n",A), 

prinlfl:"Rs0=%6.2f\n",R50); 

pnnlH"n=%6.2f\n",n); 

delfa=0.4093’sin(2.*PI*(284.+DOY)/365.); 

ws=acos(-tan(Lal*PI/180.)’tan(ddta)); 

N-ws*(24./PI); 

p=N*100./totaIN; 

if(RHmin>100.)  RHmin=IOO.; 
if(RHmin<0.)  RHmin=0.; 
if|n>N)  n=N; 
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if(Ud>10.)  Ud-10.; 
a=0.0043  *RHmin-(n/N)-  1 .41; 
b=0.82-0.004 1 •RHmin+ 1 07>n/N+0.066,Ud 


-0.006*RHmin*n/N-0.0006*RHmin*Ud; 
B_C_ET=a+b*p,(0.46*airt+8. 13); 

printf("p=%6.2f\n*,p); 

printf("N=%6.2f\n",N); 
prinlJi"a=%7  3i\n",a), 
prinH"b=%7.3ftn",b); 

printfr"B_Cdle  ETO  = %8.3fmm/day\n",B_C_ET); 


FAO_Pan(void)  ( 
float  kp,  ETpan ; 
float  Fetch=1000.; 

kp=0.475-0.00024*U2+.005 1 6’RHmcatH-OO  1 18’Fetch 
-.00001 6*RHmean*RHmcan-.000001 01  'Fetch'Fetch 
-000000008"RHmean"RHmean'U2-.00000001"RHmean"RHmean*Fetch; 
ETpan-kp'pan, 
printfTkp=%7  2l\n',kp); 

printfl"FAOpan  ETO  = %8.3f  mm/day\n",ETpan); 


Pcnman_Monteith(void)  j 
float  eb,VPD,gamma,Delta,LE,G=0.21; 
float  ra,rc,gammap,K,P_M_ET; 
float  d,zOm.zOv,hc,zw,zp, 

rc=25.;  //canopy  resistance  (s/m)  from  Albetew  1995 
zw=2000.;  //wind  speed  measured  at  heigh  zw  (mm) 
zp=1200.;  //temperature  & humidity  measured  at  heigh  zp  (mm) 
hc=700.;  //mean  crop  height  for  cattail  & sawgrass  (mm) 

LE=2.501-0.002361*airt; 

e0  = 0.611  "exp((l 7.27  • airt)/(airt+23 7.3)); 

VPD=e0*(  1 ,-(RHmax+RHmin)/200.); 

Delta  - (4098.*e0)/((airt+237.3)*(airt+237.3)); 

gamma  = (Cp  • P)/(0.622,LE)’  1 E-3; 

z0m=0.123*hc; 

z0v=0.1*z0m; 

d=0.66667*hc; 

raKlogttzw-dj/zOmJ’logftzp-dVzOvMA  1 *.4 1 *U2); 


K=1710.-6.85*ain; 
ganunap=gamma*( ! ,+rc/ra); 

P_M_ET=((Delta/(Della+giunmap)),(Rn-G)+<gamnia/(Dclta+gammap)),K/ra*VPD)/LE-, 
//"  prinII("K=  %7.2l\n',K); 

II  prinlH'ra=  %7.3ftn",ra); 

//  printH'gammap=%6.3t\n",gammap); 

primfTPMnlh  ETO  = %8.3f  mm/day\n",P_M_ET); 
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